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ABSTRACT 


This  report  describes  the  research  performed  over  the  period  of  1  January  1999  to  31  December 
1999  on  a  MURI  under  Office  of  Naval  Research  contract  N00014-96-1-1 173  on  the  topic  “Acoustic 
Transduction:  Materials  and  Devices”  brings  together  work  from  the  Materials  Research  Laboratory 
(MRL),  the  Applied  Research  Laboratory  (ARL),  and  the  Center  for  Acoustics  and  Vibrations  (CAV)  at 
The  Pennsylvania  State  University.  As  has  become  customary  over  many  years  research  on  the  program  is 
described  in  detail  in  the  87  technical  appendices  to  this  report  and  only  a  brief  narrative  description 
connecting  this  research  is  given  in  the  text. 

Perhaps  the  most  outstanding  accomplishment  of  the  year  is  a  “spin  on”  from  our  earlier  single 
crystal  studies  now  involving  Brookhaven  National  Laboratory  and  Professor  Gonzalo  s  group  in  Madrid, 
Spain.  Using  exceptionally  homogeneous  polycrystal  lead  zirconate  titanate  samples  prepared  in  MRL, 
precise  synchotron  x-ray  analysis  has  confirmed  a  new  monoclinic  phase  at  lower  temperature  in 
composition  close  to  the  important  morphotropic  phase  boundary.  This  work  demands  a  re-thinking  of 
both  intrinsic  and  extrinsic  contributions  to  response  in  this  most  important  practical  transducer  material 
family.  Domain  Engineering/Domain  Averaging  in  lead  zinc  niobate:lead  titanate  (PZN.PT)  in  lead 
magnesium  niobate:lead  titante  (PMN:PT)  and  inn  barium  titanate  (BaTi03)  continues  to  offer  single 
crystal  systems  with  outstanding  transducer  and  actuator  properties  and  new  insights  into  the  field 
induced  strain  mechanisms  in  all  perovskite  type  piezoelectrics.  Excellent  progress  with  the  new  high 
strain  irradiated  P(VDF:TrFE)  relaxor  ferroelectric  copolymer  system  has  helped  catalyze  a  new  DARPA 
initiative  in  this  area  and  a  re-awakening  of  interest  in  the  whole  area  of  electrostrictive  polymer  systems. 

A  primary  objective  of  this  MURI  grouping  was  to  help  shorten  the  time  constant  for  new 
materials  and  device  concepts  to  be  applied  in  practical  Navy  Systems.  We  believe  this  has  now  been 
realized  in  joint  work  on  the  composite  cymbal  type  flextensional  arrays  for  large  area  projectors,  and  in 
the  progress  made  towards  a  micro- tonpilz  array  system. 

Original  work  on  new  step  and  repeat  piezoelectric  high  strain  systems  continues  to  make  good 
progress  now  using  commercial  motion  rectifiers  to  produce  both  linear  and  rotary  systems  with  high 
torque  capability.  New  composite  designs  are  pushing  toward  1  mm  diameter  motors  in  the  size  regimen 
where  there  are  real  difficulties  for  conventional  electromagnetic  designs. 

A  new  area  of  activity  this  year  is  in  piezoelectric  transformers  where  a  circular  symmetry  design 
in  conjunction  with  controlled  inhomogeneous  poling  is  shown  to  offer  capabilities  which  are  of  real 
interest  for  energy  recovery  actuator  power  systems.  Basic  studies  have  evolved  a  new  environmental 
SEM  technique  for  high  resolution  domain  wall  studies  without  changing  problems.  Work  is  continuing 
on  reliable  measurements  of  electrostrictive  constants  in  simple  solids  confirming  by  both  direct  and 
converse  methods  and  permitting  the  first  generalization  of  trends  in  these  fundamentally  important 
coupling  constants. 
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Design  and  Performance  of  a 
High  Force  Piezoelectric  Inchworm  Motor 

Jeremy  Frank,  Gary  H.  Koopmann,  Weiching  Chen,  and  George  A.  Lesieutre 
Center  for  Acoustics  and  Vibration,  Penn  State  University,  University  Park,  PA  16802 


ABSTRACT 

A  linear  inchworm  motor  was  developed  for  applications  in  adaptive,  conformable  structures  for  flow  control.  The  device  is 
compact  (82  x  57  x  13  mm),  and  capable  of  unlimited  displacement  and  high  force  actuation  (150  N).  The  static  holding 
force  is  350  N.  Four  piezoceramic  stack  elements  (two  for  clamping  and  two  for  extension)  are  integrated  into  the  actuator, 
which  is  cut  from  a  single  block  of  titanium  alloy.  Actuation  is  in  the  form  of  a  steel  shaft  pushed  through  a  precision 
tolerance  hole  in  the  device.  Unlimited  displacements  are  achieved  by  repetitively  advancing  and  clamping  the  steel  shaft. 
Although  each  step  is  only  on  the  order  of  10  microns,  a  step  rate  of  100  Hz  results  in  a  speed  of  1  mm/s.  Since  the  input 
voltage  can  readily  control  the  step  size,  positioning  on  the  sub-micron  level  is  possible. 

Keywords:  Inchworm  actuator,  piezoceramic  actuator,  piezoelectric  actuator 

1.  INTRODUCTION 

The  possible  applications  of  a  high  force,  high  displacement  jinear  actuator  are  numerous,  and  the  development  of  such  a 
device  has  received  considerable  attention  in  recent  years,  as  evidenced  by  the  patents  disclosed.1*16'  Two  types  of  linear 
piezoelectric  actuators,  ‘quasi-static  inchworm’  and  ‘ultrasonic  travelling  wave’  have  been  most  common.17**'  With  the  goal 
of  improving  the  maximum  dynamic  actuation  force,  the  quasi-static  inchworm  actuator  was  developed  to  actuate  adaptive, 
conformable  structures  in  flow  control  systems  (Figure  1). 


Figure  1  Application  of  the  actuator  in  a  flow  control  structure 
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2.  DESIGN  OF  THE  INCHWORM  ACTUATOR 


In  previous  designs,  the  dynamic  actuation  force  has  been  limited  by  machining  clearance  issues  and  passive  clamping 
schemes.  With  the  current  design,  similar  to  a  design  investigated  by  Daimler-Benz,23  we  achieve  higher  forces  by  using  the 
piezoceramic  stacks  to  directly  clamp  a  shaft,  and  by  minimizing  the  tolerance  between  the  shaft  and  the  structure.  As  seen 
in  Figure  2,  the  design  concept  consists  of  four  integrated  piezoceramic  stacks,  two  for  clamping  and  two  for  extension.  The 
structure  of  the  actuator  consists  of  an  outer  frame,  and  an  inner  block  which  is  able  to  move  vertically  on  structural 
monolithic  springs.  In  step  1,  the  lower  clamping  stack  is  energized,  clamping  the  steel  shaft  to  the  inner  block.  Next,  in  step 
2  the  pushing  stacks  are  energized,  lifting  the  inner  block  and  steel  shaft  a  distance  of  10  to  15  microns.  Finally,  the  upper 
clamping  stack  is  energized,  clamping  the  shaft  to  the  outer  frame.  At  the  same  time,  the  lower  clamping  stack  and  the 
pushing  stacks  are  de-energized,  allowing  the  spring  flexures  to  return  the  inner  block  to  its  original  position,  ready  for  the 
next  step.  When  this  process  is  repeated  at  high  frequency,  motion  is  achieved. 


Steel  shaft  (pusher) 


Figure  2  Conceptual  schematic  of  the  inch  worm  actuator 

Using  21  mm  long  PZT  stacks,  poled  in  the  d33  direction,  the  blocked  force  is  2000  N  with  a  corresponding  free  defection  of 
approximately  20  microns  (see  Figure  3).  Because  of  the  small  deflection,  any  clearance  between  the  clamp  and  the  shaft 
will  result  in  a  significant  loss  of  potential  clamping  force.  For  this  reason,  the  hole  through  which  the  steel  shaft  passes  is 
reamed,  ground  and  honed  to  as  precise  a  tolerance  as  possible.  The  shaft  is  of  hardened,  precision  ground  stainless  steel, 
and  is  chosen  specifically  to  match  the  machined  hole.  This  procedure  reduces  the  tolerance  between  the  shaft  and  the  hole 
to  about  5  microns,  resulting  in  a  minimized  loss  of  clamping  force. 


Figure  3  Force/displacement  characteristics  of  the  piezoceramic  stacks 

A  picture  of  the  prototype  actuator  is  shown  in  Figure  4.  The  structure  is  cut  using  an  EDM  process  from  a  single  block  of 
titanium  alloy,  resulting  in  a  sleek  design  with  no  assembly  necessary.  Of  particular  interest  are  the  monolithic  flexures,  also 


cut  directly  out  of  the  material,  that  allow  relative  motion  of  the  moving  parts.  The  inner  block  of  the  titanium  frame,  for 
example,  is  connected  to  the  outer  frame  only  by  four  monolithic  flexures.  These  flexures  allow  vertical  motion  but  not 
horizontal  motion  of  the  inner  block,  and  provide  the  necessary  pre-compression  force  to  the  piezoceramic  stacks.  This  force 
prevents  the  stacks  from  receiving  a  tensile  load,  which  is  known  to  cause  failure.  The  two  clamping  stacks  are  also  pre¬ 
compressed  into  the  frame  with  monolithic  flexures.  The  flexures  support  the  upper  and  lower  clamp  (seen  in  Figure  4), 
which  provide  the  clamping  force  on  the  shaft  when  the  stacks  are  energized.  The  clamp  flexures  also  support  the  vertical 
shear  load,  which  must  not  be  supported  by  the  piezoceramics. 


Figure  4  Photograph  of  the  inchworm  actuator 

3.  TESTING  THE  PROTOTYPE  ACTUATOR 

The  prototype  actuator  was  tested  to  characterize  its  performance.  The  variables  of  the  signals  driving  the  stacks  included 
waveform  shape,  level,  duty  cycle  and  relative  phasing.  With  a  need  for  flexibility  in  mind,  LABVIEW  software  was  used  to 
generate  inchworm  drive  signals  through  a  PC  digital  to  analog  output  board.  A  typical  set  of  drive  signals,  seen  in  Figure  5, 
shows  that  the  pusher  stacks  are  driven  with  a  sine  wave  while  the  clamping  stacks  are  driven  with  overlapping  square  waves. 
Since  the  clamping  mechanism  is  only  active  when  the  piezoceramics  are  energized,  the  overlap  in  the  clamp  signals  is 
necessary  to  prevent  backslip  when  actuating  an  external  force.  The  use  of  square  waves  as  drive  signals  limits  the  upper 
range  of  available  drive  frequencies  by  inducing  higher  order  longitudinal  resonances  in  the  piezoceramic  stacks.22 
However,  only  a  square  wave  can  provide  the  signal  overlap  that  is  essential  to  ensure  clamping  of  the  shaft  at  ail  times. 


Figure  5  Typical  drive  signals  generated  with  LABVIEW 

A  special-purpose  calibration  fixture  is  used  for  characterization  of  the  inchworm,  as  shown  in  Figure  6.  Its  main  feature  is  a 
stiff  beam  pivoted  on  one  end  in  a  precision  bearing.  The  fixture  is  able  to  selectively  drive  loads  that  are  primarily  inertial  in 
nature,  primarily  stiffness  in  nature,  or  a  combination  of  both.  The  beam  can  be  driven  against  a  compression  spring  for  a 
stiffness  load,  or  a  mass  can  be  attached  to  the  beam  to  increase  the  inertial  load.  A  high-sensitivity  strain  gauge  load  cell 
and  an  LVDT  position  transducer  are  used  to  measure  performance.  Further,  an  analog  to  digital  input  board  is  used  to 
record  the  output  force  and  velocity  in  LABVIEW,  a  feature  of  the  experimental  setup  that  will  eventually  enable  closed-loop 
operation  of  the  inchworm. 


Figure  6  Schematic  of  the  inchworm  actuator  performance  characterization  experiment 

Three  initial  experiments  were  devised  to  characterize  the  actuator’s  speed,  holding  force  and  dynamic  force.  First,  the 
actuator  was  driven  against  no  load  and  the  velocity  was  measured  as  a  function  of  drive  frequency.  As  seen  in  Figure  7,  the 
behavior  was  fairly  linear  up  to  about  150  Hz,  increasing  to  1.2  mm/s.  Above  150  Hz  however,  the  square  clamping  signals 
began  to  cause  problems,  and  the  performance  degraded. 
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Figure  7  Free  velocity  of  the  actuator  as  a  function  of  drive  frequency 

Next,  the  static  holding  force  of  the  clamps  was  determined  by  applying  a  DC  voltage  to  one  clamping  stack  and  measuring 
the  load  at  which  the  shaft  slips  backward.  As  seen  in  Figure  8,  the  voltage  was  increased  up  to  100  V  and  a  corresponding 
increase  in  holding  force  was  observed.  The  two  clamps  were  tested  independently  and  displayed  similar  but  not  identical 
performance,  probably  because  of  variations  in  machining  tolerance.  This  indicates  that  even  tiny  imperfections  in 
machining  can  result  in  a  significant  loss  of  performance.  The  static  holding  force  of  the  top  and  bottom  clamps  are  shown  to 
be  275  N  and  375  N,  respectively. 


Bottom  Clamp 
-►-Top  Clamp 


Figure  8  Static  holding  force  of  the  clamps 

Next,  the  dynamic  drive  signals  were  used  to  actuate  the  shaft  against  a  spring  load  to  measure  the  maximum  dynamic  force, 
or  stall  load.  The  maximum  force  was  measured  as  a  function  of  the  input  voltage  to  both  the  clamping  and  pushing  stacks, 
which  was  increased  up  to  124  V.  As  shown  in  Figure  9,  the  stall  force  increased  linearly  with  the  voltage  input  to  a 
maximum  value  of  about  130  N.  The  power  amplifiers  were  unable  to  generate  square  waves  at  more  than  124  V,  but  for  the 
voltage  limitation  of  the  piezoceramics  (150  V)  the  maximum  dynamic  pushing  force  can  be  extrapolated  to  about  155  N. 
Although  this  force  is  high,  it  is  disappointing  in  that  it  is  more  than  100  N  lower  than  the  measured  static  holding  force.  The 
loss  in  dynamic  force  is  probably  a  result  of  unwanted  oscillations  of  the  piezoceramics,  caused  by  the  square  wave 
excitation. 
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Clamp,  Pusher  Stack  Voltage 


Figure  9  Maximum  dynamic  pushing  force  (stall  load)  as  a  function  of  input  voltage 

With  evidence  that  the  square  waves  were  causing  limitations  in  dynamic  holding  force,  we  tried  alternative  signal  shapes  in 
order  to  improve  performance.  First,  sine  waves  were  used  to  drive  both  the  pushing  and  the  clamping  stacks.  As  expected, 
the  motion  was  quieter  than  with  square  waves,  since  the  clamps  were  gradually  clamping  onto  the  shaft.  However,  the 
driving  force  was  limited  to  less  than  60  N  (see  Figure  10)  because  of  the  lack  of  overlap  in  the  clamp  signals.  For 
applications  where  high  speed  but  low  force  is  required,  sinusoidal  signals  could  be  used,  since  the  smoother  clamping 
mechanism  allows  the  actuator  to  be  driven  at  much  higher  frequencies.  In  the  next  part  of  the  experiment,  ramped  square 
waves  were  used  to  drive  the  clamps  in  order  to  remove  the  high  frequency  problems  caused  by  the  square  waves.  The 
ramped  square  signals  were  really  square  waves  with  sine  waves  spliced  in  between.  The  hope  was  that  the  ramped  waves 
could  provide  the  necessary  overlap  between  the  clamp  signals  but  eliminate  the  problems  caused  by  the  discontinuous 
square  wave.  As  seen  in  Figure  10,  the  ramped  waves  achieved  better  stall  loads  than  the  sinusoidal  signals,  but  did  not 
perform  as  well  as  square  waves.  From  this  experiment  it  was  determined  that  although  the  square  waves  reduce  the  dynamic 
force  capability  of  the  actuator,  they  are  needed  to  cause  the  necessary  overlap  between  the  clamping  signals. 


Clamp,  Pusher  Voltage  (V) 


Square  Waves 
Sine  Waves 
Ramped  Squares 


Figure  1 0  Dynamic  force  tests  for  alternative  clamping  signals 


4.  CONCLUSIONS 


A  linear  inchworm  motor  was  developed  for  structural  shape  control  applications.  Features  of  the  subject  device  include 
small  size  (82  x  57  x  13  mm),  unlimited  stroke,  speeds  of  about  1  mm/s,  dynamic  forces  approaching  150  N,  and  static 
holding  force  capability  greater  than  275  N.  The  device  is  machined  as  one  continuous  structure  with  monolithic  flexures, 
and  actuates  a  steel  shaft  through  its  center.  Large  displacements  are  achieved  by  repetitively  advancing  and  clamping  the 
shaft.  The  dynamic  actuation  force,  although  high,  is  limited  by  the  use  of  square  waves  in  the  driving  signals.  Other  types 
of  signals  were  tried,  but  did  not  improve  the  stall  load. 
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Abstract 

A  new  concept  in  linear  piezoelectric  actuators  is  developed  for  applications  in 
adaptive,  conformable  structures  for  flow  control.  Motivated  by  a  desire  for  high 
actuation  force  (>lkN)  and  simplified  drive  signals,  the  design  takes  advantage  of  self¬ 
locking  wedges  to  lock  the  clamping  elements.  The  concept  relies  heavily  on 
knowledge  and  manipulation  of  the  friction  coefficients  between  several  surfaces,  so 
the  choice  of  coatings  and  lubricants  are  a  major  part  of  the  investigation.  Since  the 
wedges  are  self-locking  in  one  direction,  the  actuation  force  is  limited  only  by  the  size 
of  the  piezoceramic  and  the  strength  of  the  actuator  structure.  The  device  contains  a 
single  piezoceramic  stack  (8x8x42  mm,  PZT  5H),  so  the  drive  signals  and  amplifiers 
are  drastically  simplified  from  previous  designs.  A  prototype  of  the  concept  is 
developed  and  experimentally  tested.  At  a  drive  frequency  of  200  Hz,  the  free  velocity 
is  8  mm /s  with  a  travel  of  25  mm.  An  actuation  force  of  250  N  is  achieved  with  the 
prototype.  The  wedge  concept  also  reduces  the  amount  of  precision  necessary  in 
machining  and  assembling  the  device. 


1.  Introduction 

The  possible  applications  of  a  high  force,  high  displacement  linear  actuator  are 
numerous,  and  the  development  of  such  a  device  has  received  considerable  attention  in 
recent  years,  as  evidenced  by  the  patents  disclosed.1 1  The  clamping  mechanism  in 
most  of  the  designs  consists  of  a  piezoceramic  element  clamping  directly  (more  or  less) 
onto  a  moving  part.  Since  the  displacements  of  the  piezoceramics  are  usually  on  the 
order  of  10-20  microns,  this  means  that  the  machining  accuracy  must  be  within  a  few 
microns  to  ensure  a  high  clamping  force  (Figure  la).  Such  tight  tolerances  can  be  very 
expensive,  and  any  inaccuracy  quickly  leads  to  a  loss  of  actuation  force.  The  new 
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Figure  1  a)  Direct  clamping  method  b)  Passive  wedge  clamping  method 

wedge  concept  employs  a  passive  mechanism  for  clamping,  relaxing  die  need  for  low 
surface  tolerances  and  eliminating  the  piezoelectric  clamping  stacks  (Figure  lb).  For 
certain  values  of  the  parameters  shown  the  wedge  will  lock,  preventing  the  upper  block 

from  sliding  down. 


2.  Design  of  the  Wedgeworm  Actuator 

The  motion  of  a  wedgeworm  actuator  is  quite  similar  to  that  of  previously  designed 
Shworm-type  linear  actuators17*23'  (Figure  2).  A  central  pusher  element  consisting  of 
two  halves  is  separated  by  the  piezoelectric  stack  and  held  together  by  flexure  cut 
from  the  structure.  The  flexures  are  used  to  apply  a  pre-compressive  force  to  the  stack^ 
necessary  to  keep  the  stack  from  going  into  tension.  As  seen  m  Figure  lb),  the  left 
wedge^are  constrained  to  move  horizontally  and  are  in  contact  with  the  pusher  svu-face 
such  that  the  surfaces  are  self-locking  and  the  pusher  cannot  slide  < down.  For  this 
schematic,  the  right  wedges  are  not  in  contact  with  the  pusher.  The  right  and  left 
wedges  are  rigidly  attached  and  pulled  left  or  right  such  that  only  one  set  or  the  other  is 
in  contact  wilh  the  pusher  at  any  time.  Here,  the  left  wedges  are  pulled  to  the  right  by  a 
spring  force,  keeping  them  in  constant  contact  with  the  pusher.  Motion  of  the  device  is 
a?  follows:  1)  The  piezoceramic  stack  is  de-energized  and  the  pusher  is  locked  in  place 
by  the  left  wedges.  2)  The  stack  is  energized  with  voltage.  The  bottom  of  the  pusheris 
locked  by  the  clamp  and  cannot  move  down,  but  the  top  is  able  to  extend  up.  As 
pusher  moves  up,  the  top  left  wedge  moves  to  the  right  so  that  contact  is  mamtamed 
3)  The  stack  is  de-energized.  Now,  the  top  is  locked  in  the  extended riP°SltT’b 
flexures  pull  the  bottom  of  the  pusher  up,  completmg  one  step.  This  time,  the  bottom 
left  wedge  moves  to  the  right  to  maintain  contact  with  the  pusher.  The  process  is 
repeated  at  a  high  driving  frequency  (-200  Hz)  to  give  a  high  actuation  rate  and  smooth 

border  to  use  the  wedge  concept  as  an  actuator  clamping  mechanism,  it  is  necessary  to 
understand  the  physics  that  make  the  wedges  self-locking.  The  equilibrium  equations 
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Pusher 


Piezoceramic 


1)  Locked  Position 


2)  Stack  Energized 


3)  Stack  De-Energized 


Figure  2  Operational  mechanism  of  the  wedgeworm  concept 

for  a  static  system  consisting  of  two  wedges  can  be  simplified  to  a  single  equation 
relating  two  applied  forces  on  the  wedges,  C  and  F  (Figure  3).  If  the  ratio  is  negative 
the  wedges  are  self-locking.  If  C/F  is  positive,  then  a  downward  force  F  must  be 
supported  by  a  horizontal  force  C,  so  the  wedges  are  not  self-locking.  In  Figure  3  the 
self-locking  equation  is  a  function  of  the  coefficient  of  friction  at  three  contacting 
surfaces  and  the  wedge  angle.  As  a  possibility  for  unlocking  the  wedge  for  two-way 
motion,  a  change  in  the  coefficient  of  friction  at  the  wedge  mterface  is  assumed.  Thi 
could  be  accomplished  by  extending  a  teflon  (low  friction)  pad  from  the  wedge  surface 
with  a  piezoceramic  stack.  For  the  actuator  to  have  the  maximum  amount  of  travel,  the 
wedge  angle  0  should  be  as  high  as  possible.  The  problem  is  subject  to  two  constraints, 
however.  First,  the  self-locking  equation  must  be  negative  (locked)  using  the  original 
friction  parameters.  With  the  teflon  pad  extended,  the  equation  should  be  positive  (not 
locked),  which  gives  the  second  constraint  on  the  parameters.  A  non-linear  simulated 
annealing  optimization  routine  is  used  to  maximize  the  wedge  angle  subject  to  the  self¬ 
locking  and  non-self-locking  constraints.  A  summary  of  the  optimization  problem  is 
shown  in  Figure  4.  The  objective  function  is  to  maximize 


C  (1-  AA)-sinfl-(M+/^)-C0Sfl 
F  ~  (/A  -H  /^ )  *  sin6>  +  (i  -  M/^2 )  • cos^ 

Where 

•  F  =  axial  load 

•  C  =  applied  force  on  the  wedge 

•  p.  =  coefficient  of  friction 

•  Q  =  wedge  angle 


Figure  3  Self-locking  wedge  equilibrium  equation 
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Note  that  in  the  equation,  this  means  that  pt  is  replaced  with  ps,  the  friction 
coefficient  for  teflon.  There  are  five  variables  in  the  problem:  the  wedge  angle  and 
four  friction  coefficients  (three  original  contacting  surfaces  and  one  for  teflon).  Shown 
in  Figure  4,  the  permissible  range  for  the  wedge  angle  is  from  0  to  90  degrees.  The 
chosen  range  for  the  friction  coefficients  is  0.002  to  0.78,  selected  as  maximum  and 
minimum  realistic  friction  values  from  handbooks.  This  range  has  since  been  shown  to 
be  difficult  to  realize  in  practice  for  machined  and/or  coated  surfaces,  particularly  the 
upper  limit  on  friction.  The  program  uses  a  non-linear  optimization  technique  to 
maximize  the  wedge  angle  by  varying  the  four  friction  coefficient  values.  The  solution 
always  indicated  that  the  coefficients  pi  and  p3  should  be  as  high  as  possible,  while  the 
teflon  (slipping)  coefficient  ps  should  be  as  low  as  possible  (Figure  5a).  Interestingly, 
P2>  the  vertical  friction  surface,  also  went  to  the  lower  bound,  indicating  that  the 
vertical  sliding  surface  should  be  lubricated.  With  a  reasonably  high  safety  factor,  the 
maximized  wedge  angle  was  just  over  65  degrees  (Figure  5b),  so  this  was  chosen  as  the 
actuator  wedge  angle.  To  ensure  a  smooth  horizontal  sliding  surface  for  the  wedges, 
hardened  steel  rods  were  integrated  into  the  wedgeworm  frame  (Figure  6).  The  device 
is  12  cm  in  length,  but  could  be  made  smaller  with  a  structurally  optimized  structure. 


max(0)  min(-0) 


_  ( 1  -  pl-p3)-sin0 -  (pi  +  p3>cos0  <  ^ 
”  (pi  -  p2)-sin0^(l-  pl-p2)-cos0 

(1  -  pS-p3)-sinQ  -  (pS  +  p3)-cos0  > 

X  (pS  -  p2)-sin0  i-  (1  -  pS-p2)*cos0 


O.O<0<9O.O 
0.002  <  jj.3  <  0.78 
0.002  <pS<0.78 
0.002  <  pi  <0.78 
0.002  <  p2  <  0.78 

Figure  4  Optimization  problem  statement 


3.  Testing  the  Prototype  Actuator 

The  prototype  wedgeworm  actuator  was  tested  to  characterize  its  performance.  Since 
the  wedgeworm  clamping  mechanism  is  passive,  the  drive  signal  is  a  single  frequency 
sinusoid  that  could  be  generated  with  a  simple  circuit  This  is  one  of  the  this  is  that  the 
only  values  to  vary  in  experimental  characterization  experiments  are 


Co«fici»nt  of  Friction 
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Figure  5  Optimization  results  showing  a) 


friction  coefficients,  and  b)  wedge  angle 
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Figure  6  The  prototype  wedgeworm  actuator 
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the  signal  voltage  and  frequency.  To  characterize  the  wedgeworm  prototype  the  free  velocity  (dnvrng 
against  no  load)  is  first  measured  as  a  function  of  the  drive  signal  frequency.  The  drive 
voltage  (peak-to-peak,  positively  biased)  is  held  constant  at  144  V.  This  correspon  o 
an  electric  field  of  11.5  kV/cm.  As  seen  in  Figure  7,  the  actuator  velocity  increases 
linearly  with  drive  frequency.  This  is  expected  since  the  step  size  should  remain 
constant  with  frequency,  while  more  steps  per  second  are  taken  as  the  (hive  frequency 
is  increased.  At  a  drive  frequency  of  200  Hz,  the  measured  free  velocity  is  8  1 _mm/s. 
It  seems  that  higher  speeds  would  be  possible  if  the  frequency  were  increased  fimdier, 
but  the  amplifier  power  limit  was  reached  at  200  Hz.  A  velocity  of  8.1  mm/s  at  200  Hz 
corresponds  to  a  step  size  of  40.5  microns,  which  approaches  the  theoretical  free 
displacement  of  the  stack,  42  microns.  Thus,  it  seems  that  the  backlash  under  no  lo 
is  almost  zero  for  the  wedgeworm  prototype.  Next  the  maxnnum  actuation  force,  or 
stall  load,  was  measured  as  a  function  of  drive  voltage.  The  actuator  was  driven  with  a 
100  Hz  signal  with  a  positively  biased  peak-to-peak  voltage  varying  from  65  V  to  150. 
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Drive  frequency  (Hz) 

Figure  7  Free  velocity  of  the  wedgeworm  vs.  drive  frequency 


V  (the  maximum  drive  voltage  for  die  stacks).  The  actuator  pushed  against  a 
compression  spring,  and  was  allowed  to  push  until  a  stall  condition  occurred  Using  the 
original  wedge  surfaces  (stainless  steel  on  stainless  steel  with  p.  =  0.3),  the  maximum 
actuation  force  was  240  N  for  a  150  V  drive  signal.  In  an  attempt  to  increase  the  stall 
load,  a  thin,  hard  coating  of  aluminum  oxide  was  applied  to  the  wedge  surfaces.  The 
intent  was  to  increase  the  coefficient  of  friction  at  the  wedge  interface,  thereby 
increasing  the  actuation  force.  As  seen  in  Figure  8,  the  maximum  actuation  force 
increased  with  drive  voltage  to  246  N  at  150  V,  so  the  actuation  force  was  not 
increased  significantly  with  the  aluminum  oxide  coating.  The  relationship  between 
stall  load  and  drive  voltage  is  not  quite  linear,  and  seems  in  fact  to  be  approaching  a 
limit  around  275  N.  This  indicates  that  as  the  actuation  force  becomes  very  high  there 
may  be  some  slip  in  the  wedge  clamps,  limiting  the  stall  load. 


Figure  8  Maximum  actuation  force  vs.  drive  signal  voltage 
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4.  Conclusions 


A  linear  wedgewonn  actuator  was 

'gXSZSSSSSSXXX&BA* 

Ld  a  free  velocity  of  8  mm/s  at  a  drive  frequency  of  200  Hz. 
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Abstract 

A  new  rotary  motor  that  uses  a  piezoelectric  ceramic  drive  element  is  described.  The  main 
principle  underlying  its  operation  is  rectification  of  a  small  cyclic  displacement  of  the  active 
piezoelectric  element  using  a  mechanical  diode.  The  specific  device  described,  about  the  size  of 
an  audio  cassette  tape,  consists  of  an  active  stack,  a  mechanism  by  which  to  convert  linear 
motion  of  the  stack  to  rotary  motion  of  a  shaft,  and  two  over-running  roller  clutches.  To  obtain 
rotary  motion,  the  stack  is  driven  with  an  oscillatory  signal.  On  the  driving  half  of  each  cycle, 
the  forward  motion  of  the  working  end  of  the  stack  is  converted  to  rotation  of  the  shaft  when  the 
stack  torque  exceeds  that  of  the  load.  On  the  recovery  (over-running)  half  of  each  cycle,  a 
second,  fixed,  roller  clutch  prevents  the  load  from  back-driving  the  shaft.  Experiments  to  date 
have  demonstrated  steady  torques  of  about  0.5  N-m  and  speeds  of  about  450  RPM. 

Introduction 

Considerable  research  has  addressed  the  use  of  smart  (or  active)  materials  such  as  piezoelectric 
ceramics  for  solid-state  induced-strain  actuation  [Crawley  and  de  Luis,  1987].  The  actuation 
energy  density  of  an  active  material,  approximately  proportional  to  the  product  of  blocked  stress 
and  free  strain,  is  a  suitable  figure-of-merit  for  comparing  the  relative  performance  of  alternate 
drive  materials  [Giurgiutiu  and  Rogers,  1996].  The  specific  energy  density,  the  volumetric 
energy  density  divided  by  the  material  density,  may  be  a  more  appropriate  metric  for 
applications  in  which  weight  is  more  critical  than  volume.  Corresponding  effective  material 
power  densities  may  be  defined  by  multiplying  by  a  drive  frequency. 

When  used  in  a  device,  some  kind  of  generalized  lever  ("motion  amplification")  system  is  often 
required  to  maximize  energy  transfer  ("stiffness-match")  to  a  deformable  elastic  load  [Giurgiutiu 
and  Rogers,  1997].  Additional  parasitic  volume  and  weight  is  associated  with  this  lever  system, 
making  the  energy  density  of  a  device  smaller  than  that  of  the  active  material  used  in  it. 

In  many  realistic  applications  having  volume  or  weight  constraints,  the  actuation  energy  density 
available  from  solid-state  induced-strain  materials  or  devices  is  inadequate  [Lesieutre  and 
Koopmann,  1998].  In  these  cases,  it  is  sensible  to  consider  the  use  of  smart  materials  as  the 
active  drive  elements  in  non-solid  state,  large-stroke  devices  such  as  motors.  Examples  of  such 
generalized  motors  include  inchworms  [Chen  et  al.  1997;  Duong  and  Garcia,  1996;  Frank  et  al. 
1999;  Teter  et  al.,  1998;  Zhang,  1999],  ultrasonic  wave  motors  [Glenn  and  Hagood,  1997; 
Uchino,  1998],  and  pumps  [Mauck  and  Lynch,  1999],  In  all  of  these  devices,  the  active  elements 
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are  driven  in  an  oscillatory  manner  (sometimes  at  high  frequency,  sometimes  at  resonance),  and 
some  means  of  rectification  is  provided  to  develop  one-sided  motion.  Often,  this  rectification 
method  involves  some  kind  of  generalized  diode:  pumps  use  check  valves,  while  inchworms  and 
wave  motors  use  friction-based  clamps.  Power  density  and  specific  power  density  are  more 
appropriate  figures-of-merit  for  such  motors  than  energy  density,  as  their  performance  is  usually 
rate-,  not  stroke-,  limited. 

In  a  motor-type  device,  some  kind  of  generalized  gearing  system  may  also  be  used  to  maximize 
power  transfer  ("impedance-match")  to  a  load.  Additional  parasitic  volume  and  weight  is 
associated  with  such  a  gearing  system,  as  well  as  with  the  method  used  to  rectify  oscillatory 
motion,  making  the  power  density  of  a  device  smaller  than  that  of  the  active  material  used  in  it. 

Finally,  the  development  of  smart  materials-based  motors  invites  comparisons  to  competing 
approaches  such  as  hydraulics  and  electric  motors.  Consulting  a  catalog  of  commercial  motors 
under  100  hp  (75  kW)  leads  to  a  rough  power  density  of  100  W/kg  for  typical  electric  motors. 
Electric  motors  and  piezoelectric  material-based  motors  have  an  advantage  over  hydraulic 
systems  in  that  power  may  be  transferred  over  long  distances  with  relatively  light  wires. 
Piezoelectric  material-based  motors  have  advantages  over  typical  electric  motors  in  that  they 
offer  better  potential  to  conform  with  geometric  requirements  associated  with  tightly-integrated 
adaptive  structures,  and  in  the  potential  for  reduced  electromagnetic  field  generation. 

Background:  Piezoelectric  Inchworm  Motors 

Figure  1  shows  a  canonical  piezoelectric  linear  stepping  motor  ("inchworm"),  comprising  two 
clamping  elements  (considered  to  be  in  fixed  locations)  along  with  a  central  "pusher"  element 
capable  of  length  changes  and  of  sliding  relative  to  the  clamps. 


Figure  1 .  A  piezoelectric  inchworm  motor 
One  cycle  of  motion  proceeds  as  follows: 

(1)  The  right  clamp  is  released; 

(2)  The  pusher  element,  held  in  place  on  the  left,  is  extended; 

(3)  The  right  clamp  is  engaged,  holding  the  pusher  element  in  place  on  the  right; 

(4)  The  left  clamp  is  released; 

(5)  The  pusher  element,  held  in  place  on  the  right,  is  contracted;  and 

(6)  The  left  clamp  is  engaged. 

With  typical  piezoelectric  stack  lengths  on  the  order  of  2  cm,  typical  (quasi-static)  maximum 
step  sizes  are  on  the  order  of  20  pm.  Large  strokes  may  be  obtained  by  taking  repetitive  steps. 
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If  the  steps  are  taken  at  high  rates,  substantial  stroke  rates  are  possible.  For  example,  a  step  rate 
of  1000  Hz  combined  with  a  step  size  of  20  pm  yields  a  rate  of  2  cm/sec.  Note  that  fixed- 
frequency  operation  offers  some  potential  for  optimizing  and  miniaturizing  the  necessary  drive 
electronics. 


In  order  to  compete  effectively  with  alternative  technology,  the  output  power  density  of 
piezoelectric  motors  must  be  increased.  The  keys  to  high  power  density  in  such  motors  are  high- 
frequency  operation,  along  with  rectification  of  small  solid-state  oscillatory  motion. 

Output  power  is  the  product  of  force  and  stroke  rate.  In  piezoelectric  inchworm  stepping  motors, 
the  maximum  output  force  is  the  lower  of  the  pushing  or  holding  force  capability.  The  pushing 
force  is  proportional  to  the  cross-sectional  area  of  the  driving  piezoelectric  element.  The  holding 
force  capability  of  the  clamping  mechanism  has  typically  been  the  limiting  factor  for  such 
devices,  and  has  motivated  considerable  continuing  research. 


The  maximum  stroke  rate,  however,  being  the  product  of  the  maximum  step  size  (proportional  to 
length)  and  maximum  quasi-static  step  rate  (limited  by  the  resonance  frequency,  inversely 
proportional  to  length),  is  roughly  independent  of  length  scale.  Ignoring  parasitic  mass,  the 
maximum  quasi-static  stroke  rate  is  given  roughly  by 


maximum  velocity 
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where  c  is  the  speed  of  a  longitudinal  wave  in  the  active  material,  and  emax  is  the  maximum  free 
strain.  For  typical  piezoelectric  ceramic  materials,  this  quantity  is  on  the  order  of  1.2  m/s.  Due 
to  parasitic  mass  and  other  non-idealities,  such  rates  are  not  approached  in  practice  by  quasi¬ 
static  stepping  devices. 

Thus,  increased  power  density  is  obtained  by  using  short  piezoelectric  elements  for  maximum 
drive  frequency.  An  approximate  expression  for  specific  power  density,  ignoring  parasitic 
effects,  is  given  by 
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In  addition,  in  real  devices,  the  holding  element  is  critical.  Ideally,  the  clamping  system  behaves 
as  a  mechanical  diode,  permitting  selective  motion  in  one  direction.  Real  holding  mechanisms 
exhibit  backlash  and/or  low  stiffness  over  some  small  (seating)  range  of  motion.  Because  the 
step  size  must  be  somewhat  larger  than  the  backlash  in  the  clamping  system,  a  piezoelectric  drive 
element  cannot  be  so  small  that  its  free  displacement  is  smaller  than  the  backlash. 


Maximum  power  density  is  therefore  more  realistically  obtained  by  using  a  piezoelectric  drive 
element  that  is  as  short  as  practical,  as  limited  by  the  backlash,  and  running  it  at  as  high  a 
frequency  as  possible.  Figure  2  illustrates  this  concept,  with  numerical  results  based  on  the 
following  assumptions:  0.05  deg  (3  arcmin)  backlash;  a  stack  blocked  force  of  3500  N;  a 
moment  arm  of  1.5  cm;  quasi-static  driving  at  half  the  stack  resonance  frequency  (about  5  kHz 
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Figure  2.  Specific  Power  Density  is  limited  by  backlash  for  short  stacks 
and  by  low  frequency  operation  for  long  stacks 


Note  that  very  high  frequency  operation  associated  with  very  small  devices  could  also  result  in 
substantial  heating,  the  effects  of  which  have  not  been  considered.  Furthermore,  in  applications 
requiring  the  coordinated  efforts  of  many  small  devices  to  deliver  substantial  total  power  to  a 
macroscopic  load,  there  is  likely  to  be  some  additional  parasitic  mass  and  volume  involved. 

Performance  of  Some  Piezoelectric  Stepping  Motors 

Continuing  research  at  the  Penn  State  Center  for  Acoustics  and  Vibration  (CAV)  has  pursued  the 
development  of  high  force  piezoelectric  stepping  motors  [Frank,  Chen,  Lesieutre,  and 
Koopmann,  1999].  These  motors  have  evolved  considerably  over  time,  primarily  in  the 
approach  used  to  perform  the  clamping  function.  Holding  mechanisms  initially  relied  entirely  on 
Coulomb  friction  ("inchworm"),  then  on  a  combination  of  friction  and  mechanical  advantage 
("wedgeworm"),  and  currently  on  mechanical  interference  ("roller-wedgeworm"). 

This  section  summarizes  the  performance  of  several  of  these  previous  piezoelectric  stepping 
motors. 

Figure  3  shows  the  "H3C  inchworm."  It  functions  very  much  like  the  canonical  inchworm 
shown  in  Figure  1  and  described  in  the  preceding  section.  It  uses  a  direct  frictional  clamp  that  is 
initially  pre-loaded,  so  the  device  is  locked  with  power  off.  At  a  peak  drive  level  of  120  Vpp, 
and  a  step  rate  of  about  500  Hz,  it  delivers  a  maximum  45  N  of  dynamic  force,  with  a  maximum 
(no  load)  velocity  of  about  6  mm/sec.  The  static  holding  force  is  substantially  higher. 

Three  stacks  are  required,  each  with  appropriate  control  signals.  Typically,  the  clamps  are 
driven  with  near-square  waves  to  maximize  clamping  force,  while  the  pusher  is  driven  with  a 
harmonic  signal.  Performance  of  the  device  is  somewhat  sensitive  to  the  phase  of  the  pusher 
signal  relative  to  that  of  the  clamp  signals,  and  the  optimum  phase  angle  is  load-dependent.  The 
design  of  the  device  limits  its  stroke  to  about  10  mm. 
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Figure  3.  The  H3C  Inch  worm 


Figure  4  shows  the  "H3DB"  inchworm.  It  differs  from  the  H3C  in  several  ways.  It  consists  of  a 
single  block  of  material  machined  so  that  an  inner  platform  is  suspended  from  an  outer  base  by  a 
set  of  flexures. 


Figure  4.  The  H3DB  Inchworm 


The  device  operates  to  move  a  shaft  through  the  block,  with  very  large  total  stroke.  Single  stacks 
on  the  floating  platform  and  on  the  outer  base  activate  direct  frictional  clamps  onto  the  shaft. 
There  is  a  small  tight-tolerance  clearance  between  the  clamp  face  and  the  shaft.  Two  stacks 
move  the  floating  platform  relative  to  the  outer  base. 
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At  a  near-peak  drive  level  of  120  Vpp,  and  a  step  rate  of  about  100  Hz,  the  H3DB  delivers  a 
maximum  130  N  of  dynamic  force,  with  a  maximum  (no  load)  velocity  of  about  1.5  mm/sec, 
along  with  very  accurate  position  control.  Asymmetry  in  the  pusher  drive  system  is  probably  the 
cause  of  this  relatively  low  speed.  Static  holding  force  is  on  the  order  of  450  N.  Tight  tolerance 
machining  makes  this  device  relatively  expensive,  and  the  effects  of  wear  on  long-term 
performance  is  an  issue. 

The  next  device  in  the  succession  of  piezoelectric  stepping  motors  was  built  to  explore  an 
alternative  clamping  mechanism.  As  shown  in  Figure  5,  it  aims  to  exploit  the  concept  of  a  "self¬ 
locking  taper."  For  some  combinations  of  wedge  geometry  (0)  and  interface  coefficients  of 
friction  (p),  no  force  (C)  is  required  at  the  clamp  to  resist  the  load  (F).  Devices  built  using  this 
principle  were  termed  "wedgeworms." 


Figure  5.  A  Self-Locking  Taper 

Figure  6  shows  the  "Wedgeworm  I."  The  wedge  clamping  mechanism  is  intended  to  function 
essentially  as  a  mechanical  diode,  that  is,  it  cannot  be  back-driven.  Because  of  uncertainties  in 
wedge  angles  and  coefficients  of  friction,  the  wedging  mechanism  did  not  work  perfectly. 

At  a  drive  level  of  150  Vpp,  and  a  step  rate  of  about  200  Hz,  it  delivers  a  maximum  250  N  of 
dynamic  force,  with  a  maximum  (no  load)  velocity  of  about  10  mm/sec.  The  static  holding  force 
is  substantially  higher.  Because  of  the  wedge  geometry,  the  stroke  is  limited  to  about  25  mm. 
One  of  the  main  advantages  of  this  design  is  that  only  a  single  stack  is  required  for  forward 
motion.  It  can  be  effectively  driven  using  a  single  fixed-frequency  harmonic  signal.  By 
switching  the  active  wedges  from  one  side  to  the  other,  reversible  motion  may  be  obtained. 
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Figure  6.  The  Wedgeworm  I 

Figure  7  shows  a  concept  for  the  "Wedgeworm  II."  The  primary  difference  from  the 
Wedgeworm  I  is  that  the  wedging  clamps  now  move  along  with  the  central  pusher  assembly, 
instead  of  moving  perpendicular  to  the  side  wall.  These  wedges  are  held  in  place  by  preloaded 
springs.  The  pusher  of  the  Wedgeworm  II  essentially  moves  along  a  rigid  channel  that  reacts  the 
wedging  loads. 


Figure  7.  The  Wedgeworm  II 

The  next  stage  in  device  evolution  involved  recognizing  that  the  wedge,  which  still  relies  to 
some  extent  on  friction  to  operate,  could  be  replaced  by  a  rolling  element.  In  this  situation, 
mechanical  interference,  not  friction,  provides  the  primary  clamping  force.  Figure  8  illustrates 
this  concept. 
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Figure  8.  A  Linear  Roller-Wedge 

Figure  9  shows  a  device  built  to  exploit  this  concept.  This  wedging  mechanism  provides  an 
effective  mechanical  diode,  especially  in  an  axisymmetric  configuration  in  which  the 
interference  wedging  loads  can  be  carried  as  hoop  stresses  in  an  outer  retaining  ring. 


Figure  9.  A  Linear  Roller-Wedgeworm 


A  Rotary  Roller-Wedgeworm 

A  recent  evolutionary  design  is  a  rotary  roller-wedgeworm.  This  device  is  built  around  a 
commercial  over-running  roller  clutch,  a  rotary  version  of  the  linear  roller-wedge  shown  in 
Figure  8.  Figure  10  shows  a  schematic  of  a  roller  clutch.  Initial  experimental  measurements  of 
roller  clutch  performance  indicate  a  typical  angular  backlash  of  about  0.3  deg. 

Figure  1 1  shows  a  rotary  roller-wedgeworm  device.  It  consists  of  an  active  stack,  a  mechanism 
to  convert  linear  motion  of  the  stack  to  rotary  motion  of  a  shaft,  and  two  roller  clutches.  To 
obtain  rotary  motion,  the  stack  is  driven  with  an  oscillatory  signal.  On  the  driving  half  of  each 
cycle,  the  forward  motion  of  the  working  end  of  the  stack  is  converted  to  rotation  of  the  shaft 
when  the  stack  torque  exceeds  that  of  the  load.  On  the  recovery  (over-running)  half  of  each 
cycle,  a  second  roller  clutch  prevents  the  load  from  backdriving  the  shaft. 
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Figure  1 1 .  A  Rotary  Roller-Wedgeworm 


The  device  shown  in  Figure  11,  driven  between  0  and  265  V  (corresponding  to  an  electric  field 
of  about  10  kV/cm),  and  at  450  Hz,  developed  a  free  angular  velocity  of  300  RPM,  and  a 
maximum  torque  of  about  0.32  N-m.  Figure  12  shows  the  measured  dependence  of  angular 
velocity  on  step  rate.  Other  researchers  have  also  recently  begun  to  investigate  similar 
approaches  [Zhang,  1999]. 

Recent  experiments  on  other  similar  devices  have  measured  steady  torques  of  about  0.5  N-m  and 
speeds  of  about  400  RPM. 


9 


Figure  12.  Angular  Velocity  vs.  Step  Rate  for  a  Rotary  Roller-Wedgeworm  (Torque  =  0.05  N-m) 


Summary 

In  many  realistic  applications,  the  actuation  energy  density  available  from  solid-state  induced- 
strain  materials  is  inadequate.  In  such  cases,  it  is  sensible  to  consider  the  use  of  such  materials  as 
the  active  drive  elements  in  non-solid  state,  large-stroke  devices  such  as  motors.  Furthermore, 
the  development  of  such  devices  invites  comparisons  to  competing  approaches  such  as 
hydraulics  and  electric  motors. 

In  order  to  compete  with  alternative  technology,  the  output  power  density  of  piezoelectric  motors 
must  be  increased.  The  key  to  high  power  density  in  such  motors  is  rectification  of  small  solid 
state,  high-frequency  oscillatory  motion.  Maximum  power  density  is  obtained  by  using  a 
piezoelectric  drive  element  that  is  as  short  as  practical,  as  limited  by  the  backlash,  and  running  it 
at  as  high  a  frequency  as  possible. 

The  evolution  of  a  number  of  piezoelectric  stepping  motors  was  described.  These  motors  differ 
primarily  in  the  approach  used  to  perform  the  clamping  function.  Holding  mechanisms  initially 
relied  entirely  on  Coulomb  friction  ("inchworm"),  then  on  a  combination  of  friction  and 
mechanical  advantage  ("wedgeworm"),  and  currently  on  mechanical  interference  ("roller- 
wedgeworm"). 

A  new  rotary  motor  ("Rotary  Roller-Wedgeworm")  was  described,  and  consists  of  an  active 
stack,  a  mechanism  to  convert  linear  motion  of  the  stack  to  rotary  motion  of  a  shaft,  and  two 
roller  clutches.  To  obtain  rotary  motion,  the  stack  is  driven  with  an  oscillatory  signal.  On  the 
driving  half  of  each  cycle,  the  forward  motion  of  the  working  end  of  the  stack  is  converted  to 
rotation  of  the  shaft  when  the  stack  torque  exceeds  that  of  the  load.  On  the  recovery  half  of  each 
cycle,  a  second  roller  clutch  acts  effectively  as  a  mechanical  diode,  and  prevents  the  load  from 
backdriving  the  shaft.  Experiments  to  date  have  demonstrated  steady  torques  of  about  0.5  N-m 
and  speeds  of  about  400  RPM. 
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The  letter  describes  a  piezoelectric  mocor  that  combines  the  merits  of  piezoelectric  materials,  such 
as  hi2h  power  density  generated  at  electromechanical  resonance,  and  a  precise  control  ot 
displacement.  The  motor  utilizes  a  direct  coupling  mechanism  between  the  stator  and  rotor,  where 
a  clutch  drives  the  rotor  via  locking  it.  The  direct  coupling  makes  it  possible  to  transmit  the  whole 
power  venerated  in  the  piezoelectric  element  to  the  rotor,  and  thus  achieve  the  high  efficiency  ot  e 
mocor  It  also  allows  the  combining  of  two  regimes  of  operation:  continuous  rotation  and  a  stepwise 
motion  within  a  360°  interval  with  a  high  resolution  of  angular  dispiacemenL  ©  1999  American 
Institute  of  Physics.  [S0003-695 1(99)03232-5] 


Among  actuator  materials,  piezoelectrics  such  as  lead 
zirconate  titanate  (PZT)  are  distinguished  by  the  possibility 
of  die  precise  control  of  produced  displacement  and  by  a 
hi sh  power  density  which  enables  operation  at  high  driving 
frecuencies.1  The  applications  that  make  use  of  these  fea¬ 
tures  include  for  example,  robotics.1  vibration  and  noise  can¬ 
cellation.  and  translational  stages  for  optical  systems  and 
scanning  tunneling  microscopes.—5  At  the  same  time,  the  ma¬ 
jor  demerit  of  PZT  is  that  the  piezoelectric  strain  is  very 
small.  Therefore,  in  order  to  convert  a  small  strain  into  a 
larae  displacement,  possible  solutions  include  either  a 
built-in  amplification  mechanism  such  as  a  ’*moonie  or  a 
tubuiar  torsional  actuator/  or  an  accumulation  of  the  dis¬ 
placement  over  many  periods  of  the  alternating  current  (ao 
dhvins  voltaae  such  as  with  ‘"inchworm  actuators  or  ul¬ 
trasonic  motors.6 

Ultrasonic  motors  also  take  advantage  of  the  high  power 
density  of  PZT,  and  operate  at  the  electromechanical  reso¬ 
nance  of  the  piezoelectric  element,  which  is  built  into  the 
stator.  The  displacement  is  accumulated  by  converting  the 
hiah  frequency  mechanical  vibrations  excited  at  the  reso¬ 
nance  into  the  motion  of  the  rotor  by  using  the  frictional 
contact  between  the  stator  and  the  rotor.6  However,  the  fric¬ 
tional  contact  limits  the  performance  of  the  motor,  (i)  the 
energy  dissipation  due  to  the  friction  reduces  its  efficiency, 
(ii)  the  torque  that  can  be  transmitted  from  the  piezoelectric 
element  is  limited  by  the  frictional  force,  and  (iii)  sliding  of 
the  rotor  due  to  its  inertia  does  not  enable  an  instantaneous 
start  and  stop  operation,  and  hence,  it  is  difficult  to  achieve 
precise  control  over  the  displacement  of  the  rotor. 

The  concept  of  the  motor  proposed  in  this  work,  it  so 
employ  the  high  power  density  of  PZT.  by  driving  the  motor 
at  resonance,  and  to  accumulate  the  produced  displacement 
bv  usina  a  direct  coupling  mechanism  between  the  stator  and 
rotor,  which  drives  the  rotor  via  locking  it.  In  this  wav.  all 
the  power  generated  in  PZT  at  the  resonance  can  be  trans- 
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mitred  to  the  rotor  because  now  there  is  no  energy  loss  in  the 
frictional  contract,  and  the  friction  force  does  not  limit  the 
torque.  Also,  the  locking  mechanism  allows  smooth  motion 
either  in  a  continuous  or  stepwise  fashion  within  a  360°  in¬ 
terval  with  precise  control  over  angular  positioning.  Tnus.  in 
this  motor  one  fully  employs  the  merits  of  piezoelectric  ma¬ 
terials. 

In  the  proposed  motor.  Fig.  1.  piezoelectric  ceramic 
tube.  Fig.  1(a).  generates  high  frequency  torsional  vibrations 


FIG.  1.  (a)  Piezoelectric  ceramic  segments  are  bonded  together  to  form  a 
tube  using  a  conductive  adhesive,  which  also  acts  as  electrodes  to  apply  the 
driving  electric  field.  E.  Since  the  direction  of  the  polarization.  P. .  alter¬ 
nates  between  the  adjacent  segments,  they  exhibit  a  coherent  shear  piezo¬ 
electric  strain  S?  =  d,5£.  and  that  results  in  a  torsional  deformation  of  the 
tube  on  an  angle  0.  lb)  Roller  clutch  transmits  the  displacement  to  the  rotor 
onlv  if  it  moves  in  the  counterclockwise  direction,  when  the  rollers  are 
wedged  between  the  tilted  slope  of  the  cam  and  the  surface  of  the  rotor.  In 
the  clockwise  direction,  the  rollers  disengage  and  do  not  transmit  the  dis¬ 
placement.  (c)  The  proposed  piezoelectric  motor  consists  of  piezoelectric 
cube  shown  in  (a).  “I**,  two  clutches,  shown  in  (b).  “2**  and  “3**.  housing 
for  clutches  ‘*4’*  and  **5’“.  metal  shim.  *‘7*\  and  the  rotor.  “6’L  Bolt  **S” 
ciemens  together  piezoelectric  tube  and  a  metai  cylinder  ‘*9".  which  is 
required  in  order  to  have  zero  angular  displacement  in  the  shim.  ’*7”.  at 
resonance. 
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in  the  stator  using  the  shear  piezoelectric  response,  as  was 
described  in  earlier  works.5’'  A  one-way  clutch.  Fig.  i(b)  and 
**2’*  in  Fig.  1(c).  is  connected  to  the  piezoelectric  tube.  ‘*1/* 
via  housing.  “4/*  so  that  it  can  move  together  with  the  tube, 
and  is  used  to  accumulate  the  angular  displacement  in  one 
direction.  When  the  tube  with  the  dutch  twists  counterclock¬ 
wise.  Figs.  1(b)  and  1(c).  the  clutch  locks  the  rotor  with  the 
rollers  [because  the  rollers  are  wedged  between  the  tilted 
slope  of  the  cam  and  the  surface  of  the  rotor,  Fig.  1(b)],  and 
transmits  the  displacement.  On  the  contrary,  if  the  tube  twists 
clockwise,  the  clutch  releases  the  rollers  and  they  do  not 
transmit  the  displacement  to  the  rotor  in  this  direction.  Thus, 
the  clutch  provides  an  accumulation  of  angular  displacement 
in  one  direction.  An  identical  second  clutch,  “3’'  in  Fig. 
1(c).  is  mounted  to  the  external  support,  and  is  used  to  pre¬ 
vent  motion  in  the  opposite  direction  [clockwise  in  Fig.  1(b)] 
when  the  external  torque  load  is  applied  to  the  rotor. 

To  build  the  prototype.  Fig.  1(c).  the  piezoelectric  ce¬ 
ramic  tube  1  was  tightly  clamped  between  two  metal  pans,  4 
and  *'9.'*  using  “8.”  The  dimensions  of  these  components 
were  evaluated  using  an  equivalent  circuit  analysis  of  piezo¬ 
electric  vibrations.8  They  must  be  such  that  at  resonant  fre¬ 
quency.  the  lengths  of  the  pan  above  a  metal  shim,  4‘7.”  and 
below  it  are  each  equal  to  \J4:  where  \  is  the  wavelength  of 
standing  shear  wave  excited  by  the  piezoelectric  tube.  At  this 
condition,  the  maximum  anguiar  displacement  occurs  at  the 
free  ends  of  the  stator: 

£ 

jS-aT-xf  dt5E)XQL.  (1) 

where  i.;  is  the  shear  piezoelectric  coefficient.  £  is  the  elec¬ 
tric  driving  field.  L  is  the  length  of  the  pan  above  the  shim. 
Roul  is  the  outer  radius  of  the  rube  and/or  stator,  and  QL  is 
the  mechanical  quality  factor  of  the  stator.3  At  the  same  time. 
/?=Q  at  the  shim,  and  therefore,  the  shim  can  be  used  for  the 
attachment  of  the  stator  to  the  external  suppon.  The  piezo¬ 
electric  tube  was  made  of  hard  PZT  ceramics  (APC-841 
from  American  Piezoceramics),  and  metal  parts  were  made 
of  hare  aluminum  in  order  to  have  a  higher  mechanical  qual¬ 
ity  factor  Ql.  The  approximate  dimensions  of  the  stator 
were  equal  to:  2  x  L  -  80  mm  (total  length)  and  Rout**$  mm. 
Tne  resonant  frequency, /r.  corresponding  to  these  dimen¬ 
sions  is  approximately  equal  to  10  kHz.  Both  clutches.  2  and 
3.  are  commercially  available  products  from  Torrington. 

Figure  2  demonstrates  proof  of  the  concept  of  the  pro¬ 
posed  motor.  When  an  ac  voltage  was  applied  to  the  piezo¬ 
electric  tube,  the  rotor  was  moving  with  the  revolution  speed 
ft-  which  increased  toward  saturation  with  increasing  voltage 
(open  circles).  This  behavior  of  ft  correlates  with  the  voltage 
dependence  of  the  amplitude  of  the  torsional  vibrations  0 
(closed  circles).  The  latter  was  measured  at  the  free  end  of 
the  stator  by  attaching  a  small  mirror  on  top  of  the  housing, 

A  and  by  measuring  the  distance  change  between  the  mirror 
and  the  optical  probe  of  a  MTI-2000  fotonic  sensor.3 

Two  more  features  should  be  noticed  in  Fig.  2.  The  first 
feature  is  chat  from  Eq.  ( I ).  one  would  expect  that  0,  and 
therefore.  Cl.  should  increase  linearly  with  the  voltage  (or 
electric  field).  The  deviation  from  the  linear  function,  seen  in 
Fig.  2.  is  due  to  the  strong  frequency  and  field  dependence  of 
the  piezoelectric  response  around  the  resonance.9  which  is 
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FIG.  2.  .Amplitude  of  torsional  vibrations  of  the  stator.  0,  and  the  revolution 
speed  of  the  rotor,  fl.  as  a  function  of  the  ac  driving  voltage  (the  frequency 
is  constant.  9.8  kHz). 

illustrated  in  Fig.  3.  When  the  amplitude  of  the  driving  volt¬ 
age  (electric  field)  is  increased,  the  position  of  the  resonance 
peak  in  0( f)  shifts  toward  lower  frequencies.  Fig.  3(a),  due 
to  the  softening  of  the  piezoelectric  material  under  the  high 
stress  developed  in  it  at  the  resonance.3,9  As  a  result,  this 
shift  in  the  resonance  peak  will  lead  to  the  saturation  in  the 
field  dependence  of  the  displacement  measured  at  a  fixed 
frequency  (for  example.  /=  9.9  kHz:  as  is  clear  from  a  com¬ 
parison  of  the  plots  in  Figs.  3(a)  and  3(b).  At  the  same  time, 
the  displacement  at  the  peak.  :’=/r.  still  follows  a  linear 
dependence  predicted  by  Eq.  f 1 i.  as  shown  with  open  circles 
in  Fig.  5< b). 

The  second  feature  is  that  the  rotation  had  some  thresh¬ 
old  and  m  the  studied  prototype  the  rotor  did  not  move  when 
the  driving  voltage  was  below  30  V rms:  Fig.  2.  This  thresh¬ 
old  is  due  to  the  backlash,  or  lost  displacement.  10.  which  is 
a  common  feature  of  clutches  available  on  the  market.  One 
can  evaluate  10  from  the  data  shown  in  Fig.  2.  Without 
backlash,  the  amplitude  of  the  torsional  vibrations  and  the 
revolution  speed  should  be  related  by:  ft •  360°  =  2X/?X/, 
where  /  is  the  frequency  of  the  ac  voltage  and  the  factor  2 
appears  because  the  clutch  drives  the  rotor  during  half  of  a 
period  of  the  ac  voltage,  when  the  displacement  changes 
from  -0  to  +£  With  the  lost  displacement  ft  x  360° 
=2 0—10,  and  therefore,  we  obtain:  4/2^0.01°. 

Even  though  the  backlash  is  fairly  small,  it  sets  some 
limits  on  the  performance  of  the  prototype  motor.  Still,  it  is 
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FIG.  3.  .Amplitude  of  torsional  vibrations  of  the  stator,  0,  around  the  reso¬ 
nance.  i  a/  Frequency  dependence  of  0  measured  at  different  levels  of  the  ac 
driving  field.  E.  showing  the  peaks  at  the  resonant  frequency.  fr .  (b)  Field 
dependence  of  0  measured  at  a  rixed  frequency.  9.9  kHz.  (plus  signs),  and  at 
the  resonance.  t\  (open  circles i. 
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FIG.  A.  The  angular  displacement  of  the  rotor.  d>,  measured  in  the  stepwise 
motion  as  a  function  of  the  numoer  of  the  cycles.  N,  in  the  burst,  and  the  nt 
to  me  linear  law.  The  magnified  part  of  the  plot  is  shown  in  the  inset,  in 
orccr  to  demonstrate  that  there  exists  a  minimum  value  of  ;V  required  to 
move  the  rotor,  and  the  minimum  achievable  displacement,  .  The  fre¬ 
quency  of  the  ac  voltage  is  9.6  kHz. 


surprising  chat  the  motor  could  operate  even  using  commer¬ 
ce  oily  available  roller  clutches  which  in  principle  are  not  de¬ 
signed  for  operation  with  such  small  angular  oscilladons.  of 
the  order  of  0.01°  (Fig.  2).  With  improved  clutch,  the  lost 
disoiacemenc  can  be  completely  eliminated,  as  it  is  done -for 
example  in  watch  mechanisms.  This  would  allow  better  con¬ 
trol  over  the  motion,  and  would  improve  the  torque  transmit¬ 
ted  :o  the  rotor,  which  is  aiso  affected  by  the  backlash. 

Hence,  the  experiment  summarized  in  Fig.  2  proves  the 
nrcDOsed  concept,  that  ls.  to  use  the  clutch  as  a  coupling 
mechanism  to  accumulate  the  angular  displacement  gener¬ 
ated  by  the  piezoelectric  nice  and  produce  unidirectional  mo¬ 
ron  of  the  rotor.  Moreover,  because  the  clutch  drives  the 
rotor  by  locking  it.  this  ailows  achieving  not  only  a  continu¬ 
ous  rotation,  but  also  a  controlled  stepwise  motion  of  the 
rotor  over  a  finite  angle. 

In  order  to  produce  the  stepwise  motion,  the  driving 
voltage  can  be  applied  in  the  form  of  the  tone  burst,  using  a 
given  number.  N,  of  ac  cycles  at  a  frequency  /.  The  angular 
displacement.  d>.  of  the  rotor  accumulated  over  /V  cycles  can 
be  written  as:  d>  =  [2/3- A3]x.V.  where  /?  is  the  amplitude 
of  the  vibrations  of  the  stator  corresponding  to  the  same  fre- 
cuencv  of  the  driving  voltage,  and  we  took  into  account  the 
backlash.  1(3.  Figure  4  shows  the  example  of  the  angular 
displacement  as  a  function  of  .V.  measured  using  the  ac  volt¬ 
age  with  frequency  9.6  kHz  and  amplitude  62  V  rms.  The  fit 
of  die  linear  relationship.  d*.V.  to  the  data  yields  the  slope 
ecual  to  0.058^0.002  [deg],  which  has  a  meaning  of  the 
disolacement  of  the  rotor  per  cycle  of  driving  voltage.  Using 
trie  obtained  value  of  the  slope,  c biN.  and  /?= 0.04°  mea¬ 
sured  at  9.6  kHz.  we  can  calculate  A/?**0.02°T  which  is  in  a 
reasonable  agreement  with  the  value  A/?** 0.01°  estimated 
from  the  data  shown  in  Fig.  2. 


The  resolution  limit  for  the  stepwise  motion  can  be  de¬ 
termined  by  reducing  the  number  of  cycles,  N ,  and  measur¬ 
ing  the  smallest  achievable  displacement  of  the  rotor,  <bm in. 
At  62  V  rms  and  9.6  kHz.  we  found  that  this  prototype  motor 
could  be  easily  controlled  down  to  ^mjn“0.5o,  as  shown  in 
the  inset  to  Fig.  4.  By  changing  the  driving  voltage  to  25 
V  rms.  we  could  achieve  even  better  control,  0min— 0.2°, 
because  of  the  increase  in  the  number  of  cycles  required  to 
move  the  rotor  {N-50.  compared  to  N=20  in  Fig.  4). 
Hence,  this  experiment  demonstrated  that  even,  with  the 
available  clutch,  the  motor  can  produce  a  stepwise  motion 
with  a  precise  control  over  angular  positioning. 

In  previous  works.:*'  it  was  shown  that  at  low  frequen¬ 
cies  the  torque  developed  by  the  piezoelectric  tube  is  de¬ 
scribed  by 
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-^ou  ,-<) 
-■*  -U^  our 


Xd!5£. 


(2) 


where  Rm  and  Roul  are  inner  and  outer  radii  of  the  tube, 
respectively,  and  s is  the  shear  elastic  compliance  of  PZT. 
The  analysis  of  the  resonance  vibrations  using  the  equivalent 
circuit  approach3  predicts  that  the  same  formula  should  be 
aiso  valid  around  the  resonance,  which  makes  it  possible  to 
evaluate  the  torque  generated  by  the  stator  from  the  dimen¬ 
sions  of  the  tube,  and  matenais  properties,  s ^  and  al5 .  One 
can  expect  that  after  elimination  of  the  backlash,  the  dutch 
will  transmit  the  torque  to  the  rotor  completely  because  the 
motor  utilizes  the  direct  coupling  mechanism,  via  locking  the 
rotor.  Thus,  the  performance  of  the  motor  will  be  fully  de¬ 
termined  by  the  piezoelectric  element,  since  now  the  fric¬ 
tional  force0  does  not  limit  the  torque. 

In  conclusion,  a  design  of  a  piezoelectric  motor  was  pro¬ 
posed.  which  takes  advantage  of  the  torsional  motion  gener¬ 
ated  at  the  resonance  in  the  stator,  and  of  the  direct  coupling 
between  the  rotor  and  the  stator  to  accumulate  the  displace¬ 
ment  and  produce  a  large  torque.  It  was  shown  that  the  direct 
coupling  allows  for  smooth  morion  either  in  the  continuous 
or  in  the  stepwise  fashion  within  a  560°  interval  with  a  pre¬ 
cise  control  over  angular  positioning. 
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This  paper  reviews  recent  developments  of  compact  ultrasonic  motors  using  piezoelectric 
resonant  vibrations.  Following  the  historical  background,  ultrasonic  motors  using  the  stand¬ 
ing  and  traveling  waves  are  introduced.  Driving  principles  and  motor  characteristics  are 
explained  in  comparison  with  the  conventional  electromagnetic  motors. 


Keywords:  ultrasonic  motor;  piezoelectric  actuator;  standing  wave  motor 


INTRODUCTION 

In  office  equipment  such  as  printers  and  floppy  disk  drives,  market  research 
indicates  that  tiny  motors  smaller  than  1  cm3  would  be  in  large  demand  over  the 
next  ten  years.  However,  using  the  conventional  electromagnetic  motor 
structure,  it  is  rather  difficult  to  produce  a  motor  with  sufficient  energy 
efficiency.  Piezoelectric  ultrasonic  motors,  whose  efficiency  is  insensitive  to 
size,  are  superior  in  the  mm-size  motor  area. 

In  general,  piezoelectric  and  electrostrictive  actuators  are  classified  into 
two  categories,  based  on  the  type  of  driving  voltage  applied  to  the  device  and 
the  nature  of  the  strain  induced  by  the  voltage;  (1)  rigid  displacement  devices  for 
which  the  strain  is  induced  unidirectionally  along  an  applied  dc  field  (servo 
displacement  transducers  and  pulse  drive  motors),  and  (2)  resonating 
displacement  devices  for  which  the  alternating  strain  is  excited  by  an  ac  field  at 
the  mechanical  resonance  frequency  (ultrasonic  motors).  The  AC  resonant 
displacement  is  not  directly  proportional  to  the  applied  voltage,  but  is,  instead, 
dependent  on  adjustment  of  the  drive  frequency.  Although  the  positioning 
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accuracy  is  not  as  high  as  that  of  the  rigid  displacement  devices,  very  high 
speed  motion  due  to  the  high  frequency  is  an  attractive  feature  of  the  ultrasonic 
motors. 

The  materials  requirements  for  these  classes  of  devices  are  somewhat 
different,  and  certain  compounds  will  be  better  suited  for  particular  applications. 
The  servo-displacement  transducer  suffers  most  from  strain  hysteresis  and, 
therefore,  a  PMN  electrostrictor  is  preferred  for  this  application.  The  pulse- 
drive  motor  requires  a  low  permittivity  material  aiming  at  quick  response  with  a 
limited  power  supply  rather  than  a  small  hysteresis,  so  that  soft  PZT 
piezoelectrics  are  preferred  to  the  high-permittivity  PMN  for  this  application.  On 
the  contrary,  the  ultrasonic  motor  requires  a  very  hard  piezoelectric  with  a  high 
mechanical  quality  factor  Qm,  in  order  to  minimize  heat  generation  and 
maximize  displacement.  Note  that  the  resonance  displacement  is  equal  to 
a-dEL,  where  d  is  a  piezoelectric  constant,  E,  applied  electric  field,  L,  sample 
length  and  a  is  an  amplification  factor  proportional  to  the  mechanical  Q. 

This  paper  deals  with  ultrasonic  motors  using  resonant  vibrations, 
putting  a  particular  focus  on  miniaturized  motors.  Following  the  historical 
background,  various  ultrasonic  motors  are  introduced.  Driving  principles  and 
motor  characteristics  are  explained  in  comparison  with  the  conventional 
electromagnetic  motors. 

CLASSIFICATION  OF  ULTRASONIC  MOTORS 
Historical  Background 

Electromagnetic  motors  were  invented  more  than  a  hundred  years  ago.  While 
these  motors  still  dominate  the  industry,  a  drastic  improvement  cannot  be 
expected  except  through  new  discoveries  in  magnetic  or  superconducting 
materials.  Regarding  conventional  electromagnetic  motors,  tiny  motors  smaller 
than  1cm  long  are  rather  difficult  to  produce  with  sufficient  energy  efficiency. 
Therefore,  a  new  class  of  motors  using  high  power  ultrasonic  energy  -- 
ultrasonic  motor,  is  gaining  wide  spread  attention.  Ultrasonic  motors  made  with 
piezoceramics  whose  efficiency  is  insensitive  to  size  are  superior  in  the  mini¬ 
motor  area.  Figure  1  shows  the  basic  construction  of  an  ultrasonic  motor. 


COMPACT  PIEZOELECTRIC  ULTRASONIC  MOTORS 


[377]/75 


which  consists  of  a  high-frequency  power  supply,  a  vibrator  and  a  slider. 
Further,  the  vibrator  is  composed  of  a  piezoelectric  driving  component  and  an 
elastic  vibratory  part,  and  the  slider  is  composed  of  an  elastic  moving  part  and  a 
friction  coat. 


V! bntor 


P  lesce  loc  t  r  I  e  Untie  Vibrator 
D  r  l  vo  r  P  ic  cq 


Fig.  1  Fundamental  construction  of  ultrasonic  motors. 


Though  there  had  been  some  earlier  attempts,  the  practical  ultrasonic 
motor  was  proposed  firstly  by  H.  V.  Barth  of  IBM  in  1973. »  A  rotor  was 
pressed  against  two  homs  placed  at  different  locations.  By  exciting  one  of  the 
horns,  the  rotor  was  driven  in  one  direction,  and  by  exciting  the  other  hom,  the 
rotation  direction  was  reversed.  Various  mechanisms  based  on  virtually  the 
same  principle  were  proposed  by  V.  V.  Lavrinenko2>  and  P.  E.  Vasiliev3>  in  the 
former  USSR.  Because  of  difficulty  in  maintaining  a  constant  vibration 
amplitude  with  temperature  rise,  wear  and  tear,  the  motors  were  not  of  much 
practical  use  at  that  time. 

In  1980's,  with  increasing  chip  pattern  density,  the  semiconductor 
industry  began  to  request  much  more  precise  and  sophisticated  positioners 
which  do  not  generate  magnetic  field  noise.  This  urgent  request  has  accelerated 
the  developments  in  ultrasonic  motors.  Another  advantage  of  ultrasonic  motors 
over  the  conventional  electromagnetic  motors  with  expensive  copper  coils,  is 
the  improved  availability  of  piezoelectric  ceramics  at  reasonable  cost.  Japanese 
manufacturers  are  producing  piezoelectric  buzzers  around  30  -  40  cent  price 
range  at  the  moment. 
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Let  us  summarize  the  merits  and  demerits  of  the  ultrasonic  motors: 


Merits 

1 .  Low  speed  and  high  torque  --  Direct  drive 

2.  Quick  response,  wide  velocity  range,  hard  brake  and  no  backlash 

--  Excellent  controllability 
--  Fine  position  resolution 

3.  High  power  /  weight  ratio  and  high  efficiency 

4.  Quiet  drive 

5.  Compact  size  and  light  weight 

6.  Simple  structure  and  easy  production  process 

7.  Negligible  effect  from  external  magnetic  or  radioactive  fields, 
and  also  no  generation  of  these  fields 

Demerits 

8.  Necessity  for  a  high  frequency  power  supply 

9.  Less  durability  due  to  frictional  drive 

10.  Drooping  torque  vs.  speed  characteristics 


Classification  and  Principles  of  Ultrasonic  Motors 
From  a  customer's  point  of  view,  there  are  rotary  and  linear  type  motors.  If  we 
categorize  them  from  the  vibrator  shape,  there  are  rod  type,  7t-shaped,  ring 
(square)  and  cylinder  types.  Two  categories  are  being  investigated  for 
ultrasonic  motors  from  a  vibration  characteristic  viewpoint:  a  standing-wave 
type  and  a  propagating-wave  type.  The  standing  wave  is  expressed  by 

us(x,t)  =  A  cos  kx  •  cos  tot,  ( ^ 

while  the  propagating  wave  is  expressed  as 

up(x,t)  =  A  cos  (kx  -  cot).  W 

Using  a  trigonometric  relation,  Eq.  (2)  can  be  transformed  as 

Up(x,t)  =  A  cos  kx  •  cos  cot  +  A  cos  (kx  -  Jt/2)-  cos  (cot  -  k/2).  (3) 
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This  leads  to  an  important  result,  i.  e.  a  propagating  wave  can  be  generated  by 
superimposing  two  standing  waves  whose  phases  differ  by  90  degree  to  each 
other  both  in  time  and  in  space.  This  principle  is  necessary  to  generate  a 
propagating  wave  on  a  limited  volume/size  substance,  because  only  standing 
waves  can  be  excited  stably  in  a  finite  size. 

The  standing-wave  type  is  sometimes  referred  to  as  a  vibratory -coupler 
type  or  a  "woodpecker"  type,  where  a  vibratory  piece  is  connected  to  a 
piezoelectric  driver  and  the  tip  portion  generates  flat-elliptical  movement. 
Figure  2  shows  a  simple  model  proposed  by  T.  Sashida.4)  A  vibratory  piece  is 
connected  to  a  piezoelectric  driver  and  the  tip  portion  generates  flat-elliptical 
movement.  Attached  to  a  rotor  or  a  slider,  the  vibratory  piece  provides 
intermittent  rotational  torque  or  thrust.  The  standing-wave  type  has,  in  general, 
high  efficiency,  but  lack  of  control  in  both  clockwise  and  counterclockwise 
directions  is  a  problem. 

By  comparison,  the  propagating-wave  type  (a  surface-wave  or  "surfing" 
type)  combines  two  standing  waves  with  a  90  degree  phase  difference  both  in 
time  and  in  space.  The  principle  is  shown  in  Fig.  3.  A  surface  particle  of  the 
elastic  body  draws  an  elliptical  locus  due  to  the  coupling  of  longitudinal  and 
transverse  waves.  This  type  requires,  in  general,  two  vibration  sources  to 
generate  one  propagating  wave,  leading  to  low  efficiency  (not  more  than  50  %), 
'  but  it  is  controllable  in  both  the  rotational  directions. 


Fig.  2  Vibratory  coupler  type  motor.  Fig.  3  Principle  of  the  propagating 

wave  type  motor. 
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CONVENTIONAL  MOTOR  DESIGNS 
Sashida  Motor 

Figure  4  shows  the  famous  Sashida  motor.5'  By  means  of  the  traveling  elastic 
wave  induced  by  a  thin  piezoelectric  ring,  a  ring-type  slider  in  contact  with  the 
"rippled"  surface  of  the  elastic  body  bonded  onto  the  piezoelectric  is  driven  in 
both  directions  by  exchanging  the  sine  and  cosine  voltage  inputs.  Another 
advantage  is  its  thin  design,  which  makes  it  suitable  for  installation  in  cameras 
as  an  automatic  focusing  device.  Eighty  percent  of  the  exchange  lenses  in 
Canon's  "EOS"  camera  series  have  already  been  replaced  by  the  ultrasonic 
motor  mechanism.  Most  of  the  studies  on  ultrasonic  motors  done  in  the  US  and 
Japan  have  been  modifications  of  Sashida's  type. 

The  PZT  piezoelectric  ring  is  divided  into  16  positively  and  negatively 
poled  regions  and  two  asymmetric  electrode  gap  regions  so  as  to  generate  a  9th 
mode  propagating  wave  at  44  kHz.  A  proto-type  was  composed  of  a  brass  ring 
of  60  mm  in  outer  diameter,  45  mm  in  inner  diameter  and  2.5  mm  in  thickness, 
bonded  onto  a  PZT  ceramic  ring  of  0.5  mm  in  thickness  with  divided  electrodes 
on  the  back-side.  The  rotor  was  made  of  polymer  coated  with  hard  rubber  or 
polyurethane.  Figure  5  shows  Sashida’s  motor  characteristics. 

Canon  utilized  the  "surfing"  motor  for  a  camera  automatic  focusing 
mechanism,  installing  the  ring  motor  compactly  in  the  lens  frame.  It  is 
noteworthy  that  the  stator  elastic  ring  has  many  teeth,  which  can  magnify  the 
transverse  elliptical  displacement  and  improve  the  speed.  The  lens  position  can 
be  shifted  back  and  forth  through  a  screw  mechanism.  The  advantages  of  this 
motor  over  the  conventional  electromagnetic  motor  are: 

1 .  Silent  drive  due  to  the  ultrasonic  frequency  drive  and  no  gear  mechanism 
(i.  e.  more  suitable  to  video  cameras  with  microphones). 

2.  Thin  motor  design  and  no  speed  reduction  mechanism  such  as  gears, 
leading  to  space  saving. 

3.  Energy  saving. 

A  general  problem  encountered  in  these  traveling  wave  type  motors  is 
the  support  of  the  stator.  In  the  case  of  a  standing  wave  motor,  the  nodal  points 
or  lines  are  generally  supported;  this  causes  minimum  effects  on  the  resonance 
vibration.  To  the  contrary,  a  traveling  wave  does  not  have  such  steady  nodal 
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points  or  lines.  Thus,  special  considerations  are  necessary.  In  Fig.  4,  the 
stator  is  basically  fixed  very  gently  along  the  axial  direction  through  felt  so  as 
not  to  suppress  the  bending  vibration.  It  is  important  to  note  that  the  stop  pins 
which  latch  onto  the  stator  teeth  only  provide  high  rigidity  against  the  rotation. 


Fig.  5  Motor  characteristics  of  Sashida’s  motor. 


Inchworm  Devices 

Although  the  motion  principle  is  different,  inchworm  devices  move  fast  in 
apparent  similarity  to  ultrasonic  motors.  The  inchworm  is  driven  by  a 
rectangular  wave  below  the  resonance  frequency,  and  moves  intermittently  and 
discretely. 
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Fig.  6  Philips’  inchworm. 


Sophisticated  linear  walking  machines  have  been  developed  by  two 
German  companies.  Philips  fabricated  a  linear  drive  inchworm  using  two  d33 
(longitudinal  mode)  and  two  d3i  (transverse  mode)  multilayer  actuators  (Fig. 
6).6>  Very  precise  positioning  of  less  than  1  nm  was  reported.  The  problems 
with  this  type  of  device  are:  (1)  audible  noise,  and  (2)  heat  generation,  when 
driven  at  high  frequency.  Physik  Instrumente  manufactured  a  two-leg 
inchworm. 7)  A  pair  of  inchworm  units  consisting  of  two  multilayer  actuators, 
are  coupled  with  90°  phase  difference  in  time  so  as  to  produce  a  smooth  motion 
instead  of  a  discrete  step  motion. 


COMPACT  MOTOR  DESIGNS 
Traveling  Wave  Types 

Using  basically  the  same  principle  as  Sashida's,  Seiko  Instruments  miniaturized 
the  ultrasonic  motor  to  as  tiny  as  10  mm  in  diameter.8)  Figure  7  shows  the 
construcion  of  this  small  motor  with  10  mm  diameter  and  4.5  mm  thickness.  A 
driving  voltage  of  3  V  and  a  current  60  mA  provides  6000  rev/min  (no-load) 
with  torque  of  0. 1  mN  m.  Seiko  installed  this  tiny  motor  into  a  wrist  watch  as 
a  silent  alarm.  Rotating  an  imbalanced  mass  provides  enough  hand  shake  to  a 
human  without  generating  audible  noise.  AlliedSignal  developed  ultrasonic 
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motors  similar  to  Shinsei’s,  which  would  be  utilized  as  mechanical  switches  for 
launching  missiles.9) 

A  significant  problem  in  miniaturizing  this  sort  of  traveling  wave  motor 
can  be  found  in  the  ceramic  manufacturing  process;  without  providing  a 
sufficient  buffer  gap  between  the  adjacent  electrodes,  the  electrical  poling 
process  (upward  and  downward)  easily  initiates  the  crack  on  the  electrode  gap 
due  to  the  residual  stress  concentration.  This  may  restrict  the  further 
miniaturization  of  the  traveling  wave  type  motros.  To  the  contrary,  standing 
wave  type  motors,  the  structure  of  which  is  less  complicated,  are  more  suitable 
for  miniaturization  as  we  discuss  in  the  following.  They  require  only  one 
uniformly  poled  piezo-element,  less  electric  lead  wires  and  one  power  supply. 


ROTOR 


SPRING  FOR 
PRESSURE 


STATOR 

VIBRATOR 


LEAD  SUPPORT  PLATE  PIEZOELECTRIC 

WIRE  FOR  STATOR  CERAMIC 


Fig.  7  Construction  of  Seiko’s  motor. 


Standing  Wave  Types 

Unear  motors 

Uchino  et  al.  invented  a  7t-shaped  linear  motor.10)  This  linear  motor  is  equipped 
with  a  multilayer  piezoelectric  actuator  and  fork-shaped  metallic  legs  as  shown 
in  Fig.  8.  Since  there  is  a  slight  difference  in  the  mechanical  resonance 
frequency  between  the  two  legs,  the  phase  difference  between  the  bending 
vibrations  of  both  legs  can  be  controlled  by  changing  the  drive  frequency.  The 
walking  slider  moves  in  a  way  similar  to  a  horse  using  its  fore  and  hind  legs 
when  trotting.  A  test  motor  20  x  20  x  5  mm3  in  dimension  exhibited  a 
maximum  speed  of  30  cm/s  and  a  maximum  thrust  of  0.9  kgf  with  a  maximum 
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efficiency  of  20%,  when  driven  at  98kHz  at  6V  (actual  power  —  0.7  W).  This 
motor  has  been  employed  in  a  precision  X-Y  stage. 

(•)  <k) 
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Fig.  8  ;c-shaped  linear  ultrasonic  motor,  (a)  construction  and  (b)  walking 
principle.  Note  the  90  degree  phase  difference  like  human  walk. 


Fig.  9  Miniature  ultrasonic  linear  motor. 


We  further  miniaturized  this  "Shepherd"  shape  by  1/10  into  a 
"Dachshund"  shape  by  reducing  the  leg  length  (Fig.  9).n)  According  to  this 
miniaturization,  we  utilized  a  unimorph  type  drive  mechanism  in  conjunction 
with  a  coupling  mode  between  1st  longitudinal  and  4th  bending  modes.  Under 
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100  Vp_p  applied  (0.6  W),  this  linear  motor  exhibited  the  maximum  speed  of 
1 60  mm/sec  and  the  thrust  of  0.25  -  1 .00  N. 

One  ceramic  multilayer  component  actuator  was  proposed  by  Mitsui 
Chemical. I2>  Figure  10  shows  the  electrode  pattern.  Only  by  the  external 
connection,  a  combined  vibration  of  the  longitudinal  Li  and  bending  B2  modes 
could  be  excited. 


J  115  ym 
-I  x  40  layers 


4.8mm 

Fig.  1 0  Multilayer  ceramic  simple  linear  motor  (Mitsui  Chemical). 

Rotanj  motors 

Hitachi  Maxel  significantly  improved  the  torque  and  efficiency  by  using  a 
torsional  coupler,  and  by  the  increasing  pressing  force  with  a  bolt.13)  The 
torsional  coupler  looks  like  an  old  fashioned  TV  channel  knob,  consisting  of 
two  legs  which  transform  longitudinal  vibration  generated  by  the  Langevin 
vibrator  to  a  bending  mode  of  the  knob  disk,  and  a  vibratory  extruder.  Notice 
that  this  extruder  is  aligned  with  a  certain  cant  angle  to  the  legs,  which 
transforms  the  bending  to  a  torsion  vibration.  This  transverse  moment  coupled 
with  the  bending  up-down  motion  leads  to  an  elliptical  rotation  on  the  tip 
portion.  A  motor  with  30mm  x  60mm  in  size  and  20  -  30°  in  cant  angle 
between  a  leg  and  a  vibratory  piece  provided  the  torque  as  high  as  1.3  N-m  and 
the  efficiency  of  80%.  However,  this  type  provides  only  unidirectional 
rotation. 

The  Penn  State  University  has  developed  a  compact  ultrasonic  rotory 
motor  as  tiny  as  3  mm  in  diameter.  As  shown  in  Fig.  1 1 ,  the  stator  consists 
basically  of  a  piezoelectric  ring  and  two  concave/convex  metal  endcaps  with 
‘’windmill"  shaped  slots  bonded  together,  so  as  to  generate  a  coupled  vibration 
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of  up-down  and  torsional  type.*4)  Since  the  component  number  and  the 
fabrication  process  were  minimized,  the  fabrication  price  would  be  decreased 
remarkably,  and  it  would  be  adaptive  to  the  disposable  usage.  When  driven  at 
160  kHz,  the  maximum  revolution  2000rpm  and  the  maximum  torque 
0.8mNm  were  obtained  for  a  5  mm<{)  motor.  Notice  that  even  the  drive  of  the 
motor  is  intermittent,  the  output  rotation  becomes  very  smooth  because  of  the 
inertia  of  the  rotor.  Figure  12  shows  motor  characteristics  plotted  as  a  function 
of  motor  size  for  modified  "windmill"  motors.15) 


Rotor 


Piezoelectric  Ring 


Fig.  1 1  "Windmill"  motor  with  a  disk-shaped  torsional  coupler. 


100  'a 

M 


0.1 

0.01 


Fig.  12  Radial  mode  resonance  frequency,  no-load  speed  and  starting  torque 
vs.  diameter  of  the  stator.  Speed  and  torque  were  measured  at  1 5.7  V. 
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INTEGRATED  MOTOR  DESIGNS 

We  will  introduce  an  ultrasonic  motor  fabricated  on  a  Si  substrate,  jointly 
developed  by  MIT  and  Penn  State. ,6>  After  coating  a  PZT  thin  film  on  a  Si 
membrane  (2.2mm  x  2.2mm),  an  8-pole  stator  (1.2mm  and  2.0mm  inner  and 
outer  diameters)  was  patterned  on  the  surface  electrode.  The  8  segmented 
electrode  pads  were  driven  in  a  four  phase  sequence  repeated  twice.  A  small 
contact  glass  lens  was  used  as  a  rotor.  This  simple  structure  provided  already 
103  time  higher  torque  than  the  conventional  electrostactic  MEMS  motors. 

Another  intriguing  surface  acoustic  wave  motor  has  been  proposed  by 
Kurosawa  and  Higuchi.17^  Rayleigh  waves  were  excited  in  two  crossed 
directions  on  a  127.8°-rotation  Y-LiNbC>3  plate  with  two  pairs  of  interdigital 
electrode  patterns.  A  slider  was  composed  of  three  balls  as  legs.  The  driving 
vibration  amplitude  and  the  wave  velocity  of  the  Rayleigh  waves  were  adjusted 
to  6. 1  nm  and  22  cm/sec  for  both  x  and  y  directions.  It  is  important  to  note  that 
even  though  the  up-down  vibrational  amplitude  is  much  smaller  (<  1/10)  than 
the  surface  roughness  of  the  LiNbC>3,  the  slider  was  transferred  smoothly.  The 
mechanism  has  not  been  clarified  yet,  it  might  be  due  to  the  locally  enhanced 
friction  force  through  a  ball-point  contact. 

SUMMARY 

Ultrasonic  motors  are  characterized  by  "low  speed  and  high  torque,"  which  are 
contrasted  with  "high  speed  and  low  torque"  of  the  conventional 
electromagnetic  motors.  Thus,  the  ultrasonic  motors  do  not  require  gear 
mechanisms,  leading  to  very  quiet  operation  and  space  saving.  Negligible 
effects  from  external  magnetic  or  radioactive  fields,  and  no  generation  of  these 
fields  are  suitable  for  the  application  to  electron  beam  lithography  etc.  relevant 
to  the  semiconductor  technology.  Moreover,  high  power  /  weight  ratio,  high 
efficiency,  compact  size  and  light  weight  are  very  promising  for  the  future 
micro  actuators  adopted  to  catheter  or  tele-surgery. 

For  the  further  applications  of  the  ultrasonic  motors,  systematic 
investigations  on  the  following  issues  will  be  required: 


86/[388] 


KENJI UCHINO  and  BURHANETTIN  KOC 


( 1 )  development  of  low  loss  &  high  vibration  velocity  piezo-ceramics, 

(2)  piezo-actuator  component  designs  with  high  resistance  to  fracture  and 
good  heat  dissipation, 

(3)  ultrasonic  motor  designs; 

a.  motor  types  (standing-wave  type,  traveling-wave  type,  hybrid  type, 
integrated  type), 

b.  simple  displacement  magnification  mechanisms  of  vibratory  piece 

(horn,  hinge-lever), 

c.  frictional  contact  part, 

(4)  inexpensive  and  efficient  high  frequency/high  power  supplies. 
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This  paper  deals  with  a  new  piezoelectric  ultrasonic  motor  operated  by  single-phase  AC.  The 
ultrasonic  motor  of  windmill  type  using  a  piezoelectric  ceramic  and  a  metal  elastic  body  was 
fabricated.  To  prepare  the  piezoelectric  ceramics  vibrator,  0.05PMN-0.95PZT  composition 
was  mixed  with  0.7 wt%  MnC^-  Dielectric  loss,  the  electromechanical  coupling  factor  of 
radial  mode  and  mechanical  quality  factor  were  0.2%,  62.8%  and  820,  respectively.  The  slot¬ 
ted  metal  elastic  bodies  made  from  brass  sheet  were  machined  at  various  thickness  of  0.15, 
0.20,  0.25, 0.30mm,  and  slots  of  the  elastic  body  were  2, 3, 4, 6,  respectively.  No-load  revolu¬ 
tion  speed  increases  with  increasing  applied  voltage,  and  it  increases  with  decreasing  the 
thickness  of  the  elastic  body.  The  revolution  speed  of  the  ultrasonic  motor  is  510rpm  under  an 
applied  voltage  of  100  Vmax 

Keywords:  Ultrasonic  motor;  piezoelectric  ceramics;  windmill  type;  single-phase  AC;  vibra¬ 
tion  velocity;  revolution  speed 


INTRODUCTION 

It  is  generally  known  that  the  first  practical  ultrasonic  motor  was 
invented  by  T.  Sashida  in  1982.  and  then  further  study  and 
development  of  the  ultrasonic  motor  were  performed1 'l  Most 
ultrasonic  motors12"'1  are  operated  by  two-phase  AC  differing  by  90° . 
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but  in  this  paper  only  single-phase  AC  is  needed  to  operate  the 
ultrasonic  motor  of  windmill  type.  We  have  studied  the  ultrasonic 
motor  considering  the  use  of  a  combination  of  radial  and  torsional 
vibration  modes14'31.  The  aims  of  the  study  are  to  develop  a  new  type 
of  ultrasonic  motor  of  windmill  type  operated  by  single-phase  AC, 
and  to  study  the  revolution  speed  characteristics  of  the  ultrasonic 
motor.  But  the  torque  characteristics  and  three-dimensional  vibration 
mode  of  real-time  animation  will  be  presented  next  time. 


EXPERIMENTAL 

Composition  of  the  piezoelectric  ceramics  is  Pb[(Mgi/3Nb2/3)oo5Zr0475 
Tio475]C>3+0.7w/%  MnC>2  and  the  elastic  body  is  70Cu-30Zn  brass. 


The  powder  was  calcined  at 
850 °C  for  2hrs,  followed  by 
sintering  at  1200°C  for  2hrs. 
The  samples  were  polished  to 
0.99mm  thickness,  and  silver 
paste  was  attached  to  both 
sides  at  700  °C.  Poling 
conditions  were  DC  1 .5 
kV/mm.  10  minute  and  115°C 
in  the  silicone  oil  bath. 
Piezoelectric  and  dielectric 
properties  were  measured 
using  an  impedance  gain- 
phase  analyzer  after  48  hrs. 
For  making  the  ultrasonic 
motor,  the  sample  dimensions 
were  1 1 .25  mm  in  diameter 
and  0.99mm  in  thickness. 


FIGURE  1.  Structure  of  the  ultrasonic  motor 
of  windmill  type. 

The  structure  of  the  ultrasonic  motor  is  shown  in  Fig.l.  where  the 
stator  made  by  piezoelectric  ceramics  and  slotted  elastic  body  vibrates 
the  radial  resonant  mode.  The  radial  vibration  mode  driven  by 
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single-phase  sinusoidal  AC  is  transformed  into  the  torsional  mode.  To 
make  a  stator,  a  piezoelectric  ceramic  vibrator  was  bonded  with  two 
pieces  of  metal  elastic  body  using  conductive  epoxy,  as  shown  in 
Fig.l.  The  elastic  bodies  with  0.15,  0.20,  0.25,  0.30mm  in  thickness 
and  with  2,3 ,4,6  in  slot  of  windmill  type  were  fabricated  using  a  saw 
machine,  respectively.  The  contact  surface  of  the  rotor  was  also  made 
of  brass.  The  characteristics  of  the  ultrasonic  motor  were  measured 
using  the  system  shown  in  Fig.2.  To  operate  the  ultrasonic  motor,  we 
used  a  high  speed  .power  amplifier(NF4015,  Japan)  and  function 
waveform  generator(HP33120A,  USA).  The  resonance  frequency  and 
vibration  velocity  of  the  stator  were  measured  with  vibrometer(Polytec 
OFV3000,  Germany)  inclusive  sensor  head(Polytec  OFV302,  Germany) 
and  digital  storage  oscilloscope(HC5802,  Korea).  The  revolution  speed 
of  the  ultrasonic  motor  was  measured  with  analog  tachometer(Ametek 
C-809,  USA). 


FIGURE  2.  Block  diagram  for  measuring  system  of  vibration 
velocity  and  revolution  speed. 


RESULTS  AND  DISCUSSION 

Table  1  shows  the  optimal  specifications  of  the  ultrasonic  motor  of 
windmill  type  which  has  6  slots  and  0.15mm  in  thickness  of  elastic 
body.  The  revolution  speed  characteristics  of  the  ultrasonic  motor  are 
basically  dependent  on  the  driving  voltage  and  frequency.  Here,  the 
resonance  frequency  of  the  stator  measured  by  vibrometer  was 
observed  at  73.6  kHz. 
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TABLE  1.  Specifications  of  fabricated  ultrasonic  motor 
of  windmill  type. 


Unit 

Constant 

Driving  frequency 

kHz 

73.6 

Driving  voltage 

vmax 

100 

Rated  speed 

rpm 

510 

Stator  weight 

gm 

1.33 

Rotor  weight 

gm 

0.18 

Resistance  of  stator 

Q 

2.53 

Fig.3  shows  the  relation  between  vibration  velocity  of  the  elastic  body 
and  the  applied  voltage  at  various  thickness  of  the  elastic  body  with  6 
slots.  From  this  result,  it  was  found  that  the  vibration  velocity  of  the 
elastic  body  increases  with  increasing  applied  voltage,  and  it  increases 
with  decreasing  the  thickness  of  the  elastic  body.  The  reason  for  this 
behavior  is  that  the  flexibility  of  the  elastic  body  increases  as 
thickness  decreases.  The  stator  of  the  ultrasonic  motor  with  0.15mm  in 
thickness  of  the  elastic  body  exhibits  a  vibration  velocity  of  2.0  m/s 
under  an  applied  voltage  of  lOOVmax. 

Fig.4  shows  the  effect  of  applied  voltage  on  no-load  revolution  speed 


Applied  voltage  (VmJ 

FIGURE  3.  Effect  of  applied  voltage  on  vibration  velocity  of 
elastic  body  at  various  thickness  of  elastic  body  with 
6  slots 
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of  the  ultrasonic  motor  at  various  thickness  of  elastic  body  with  6 
slots.  In  this  experiment,  the  elastic  body  of  brass  and  the  rotor  of 
brass  were  used.  Single-phase  AC  voltage  was  applied  to  the 
ultrasonic  motor  using  an  high  speed  power  amplifier  whose  voltage 
can  be  varied  in  the  range  of  30Vmax  to  100Vmax. 


FIGURE  4.  Effect  of  applied  voltage  on  revolution  speed  of  the 
ultrasonic  motor  at  various  thickness  of  elastic  body 
with  6  slots. 


FIGURE  5.  Effect  of  applied  voltage  on  revolution  speed  of  the 
ultrasonic  motor  at  various  slots  of  elastic 
body  for  0.15mm  in  thickness. 
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The  ultrasonic  motor  was  not  operated  below  an  applied  voltage  of 

30Vmax.  No-load  revolution  speed  increases  with  increasing  applied 
voltage,  and  it  increases  with  decreasing  the  thickness  of  elastic  body. 
The  revolution  speed  of  the  ultrasonic  motor  is  510rpm  under  an 

applied  voltage  of  lOOVma*.  Above  90Vmax,  the  revolution  speed  of  the 
ultrasonic  motor  was  almost  constant. 

Fig.5  indicates  the  effect  of  applied  voltage  on  no-load  revolution 

speed  of  the  ultrasonic  motor  at  various  slots  of  elastic  body  for 
0.15mm  in  thickness.  When  the  applied  voltage  is  100  Vmax,  the 

maximum  revolution  speed  of  the  ultrasonic  motor  is  510rpm.  The 
revolution  speed  increases  with  increasing  applied  voltage,  and  it 
increases  with  increasing  the  slots  of  the  elastic  body,  which  is  due  to 
the  generation  of  torsional  mode  at  the  elastic  body. 


CONCLUSIONS 

The  ultrasonic  motor  of  windmill  type  operated  by  single-phase  AC 
could  simply  be  fabricated  using  a  piezoelectric  ceramic  vibrator  and 
two  metal  elastic  body.  When  the  applied  voltage  to  the  ultrasonic 
motor  is  100Vmax  at  a  resonance  frequency  of  73.6kHz,  the  maximum 
revolution  speed  of  the  ultrasonic  motor  is  510rpm.  In  order  to 
improve  the  reproducibility  in  terms  of  housing  technique,  torque 
characteristics,  efficiency  and  3-dimensional  vibration  mode  of 
real-time  of  the  elastic  body  using  scanning  vibrometer  and  FEA 
method  are  now  under  investigation. 
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PIEZOELECTRIC  MICROMOTOR 


USING  A  METAL-CERAMIC  COMPOSITE  STRUCTURE 

Burhanettin  KOC,  Student  Member,  IEEE,  Philippe  BOUCHILLOUX 
and  Kenji  UCHINO,  Member,  IEEE 

Abstract-  This  paper  presents  a  new  piezoelectric  micromotor  design,  where  a 
uniformly-electroded  piezoelectric  ring  bonded  to  a  metal  ring  is  used  as  the  stator. 
Four  inward  arms  at  the  inner  circumference  of  the  metal  ring  transfer  radial 
displacements  into  tangential  displacements.  The  rotor  ends  in  a  truncated  cone 
shape  and  touches  the  tips  of  the  arms.  A  rotation  takes  place  by  exciting  coupled 
modes  of  the  stator  element,  such  as  a  radial  mode  and  a  second  bending  mode  of 
the  arms.  The  behavior  of  the  free  stator  was  analyzed  using  the  ATILA  finite 
element  software.  Torque  vs.  speed  relationship  was  measured  from  the  transient 
speed  change  with  a  motor  load.  A  starting  torque  of  17  pNm  was  obtained  at  20 
Vrms.  The  main  features  of  this  motor  are  low  cost  and  easy  assembly  owing  to  a 
simple  structure  and  small  number  of  components. 
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I.  INTRODUCTION 


Ultrasonic  motors  are  generally  categorized  as  travelling-wave  or  standing-wave 
types  according  to  the  type  of  vibration  wave  used  on  the  stator.  Currently,  most 
miniaturized  motors  utilize  a  travelling  elastic  wave  [1],  [2]  and  present  two  major 
drawbacks.  First,  the  piezoelectric  element  of  a  travelling  wave  motor,  typically  an 
annular  ring,  must  be  divided  in  sectors  poled  in  alternately  opposite  directions.  This 
makes  the  poling  process  complicated  and  damages  many  samples.  Second,  to  produce  a 
travelling  wave  on  a  vibrating  piezoelectric  element,  at  least  two  ac  power  sources  with  a 
90-degree  phase  difference  both  in  time  and  space  are  required.  This  results  in  a  more 
complex  drive  and  wiring. 

Another  important  factor  in  the  miniaturization  of  piezoelectric  ultrasonic  motors 
comes  from  the  maximum  level  of  vibrations  that  can  be  obtained  from  a  piezoelectric 
element.  The  maximum  vibration  velocity  of  the  piezoceramic  limits  the  vibration 
amplitude  [3].  As  the  size  of  a  piezoelectric  ultrasonic  motor  decreases,  the  resonance 
frequency  of  its  stator  increases.  The  vibration  level,  however,  cannot  be  increased 
beyond  a  certain  vibration  velocity  due  to  the  limitation  of  the  piezoelectric  material. 
Therefore,  a  vibration  magnification  and  transformation  mechanism  (vibration 
concentrator)  is  necessary  on  a  stator  element  to  amplify  vibration  amplitude  on  the 
piezoceramic  material  similar  to  projection  teeth  on  the  stator  elements  of  traveling  wave 
motors. 

An  alternative  design  to  miniature  piezoelectric  motors  consists  of  PZT  thin  films 
[4],  [5],  even  though  their  properties  are  inferior  to  those  of  bulk  piezoelectric  ceramics. 
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However,  the  performance  of  thin-film  motors  is  on  the  way  of  improvement  in  parallel 
with  quality  of  thin  films  [6],  [7]. 

In  order  to  miniaturize  piezoelectric  motors  we  used  a  bulk  piezoelectric  material, 
and  took  the  following  approaches:  i)  the  structure  of  the  active  piezoelectric  element, 
including  its  poling  configuration  should  be  simple  so  as  to  be  manufactured  in  small 
size;  ii)  the  number  of  components  used  in  the  motor  should  be  as  few  as  possible  to 
decrease  the  production  cost.  We  proposed  such  an  ultrasonic  motor  design  using  a 
metal-ceramic  composite  structure  [8].  The  stator  in  that  motor  consisted  of  a 
piezoelectric  ring,  poled  uniformly  in  its  thickness  direction,  bonded  to  two  metal 
endcaps.  Due  to  the  multi-mode  excitation  (radial  mode  of  the  piezoelectric  ring  and 
bending  mode  of  the  metal  endcaps),  an  elliptical  motion  was  generated  at  the  center  of 
the  top  endcap  where  the  rotor  was  placed. 

The  operating  principle  and  structure  of  the  motor  presented  in  this  work  is 
similar  to  our  previous  design.  However,  in  this  new  design,  only  one  metal  endcap  was 
used.  It  was  bonded  to  the  piezoelectric  ring  and  its  center  part  was  removed.  Finally, 
the  size  of  the  stator  vibrator  was  decreased  to  a  diameter  of  3.0  mm. 

The  structure  and  operating  principle  of  the  motor  are  provided  first,  then  the 
finite  element  analysis  of  the  stator  vibrator  is  presented.  The  motor  characteristics 
obtained  with  the  transient  response  method  are  described  last. 

n.  STRUCTURE  AND  OPERATING  PRINCIPLE  OF  THE  MOTOR 

The  motor  is  composed  of  four  components:  stator,  rotor,  ball-bearing  and 
housing  unit  (Fig.  1).  The  key  element  to  this  motor  is  a  multifunctional  stator,  which  has 
a  metal-piezoelectric  composite  structure.  The  piezoelectric  part  has  a  simple  structure  of 
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a  ring  electroded  on  its  top  and  bottom  surfaces  (<J>  3.0mm)  poled  uniformly  in  the 
thickness  direction  (Table  I).  The  metal  ring,  whose  dimensions  are  shown  in  Fig.  2,  is 
machined  using  Electric  Discharge  Machining  (EDM).  It  has  four  inward  arms  placed 
90°  apart  on  its  inner  circumference.  The  metal  and  piezoelectric  rings  are  bonded 
together,  but  the  arms  remain  free;  they  thus  behave  like  cantilever  beams.  The  length 
and  cross-sectional  area  of  each  arm  were  selected  such  that  the  resonance  frequency  of 
the  second  bending  mode  of  the  arms  is  close  to  the  resonance  frequency  of  the  radial 
mode  of  the  stator.  The  rotor  is  placed  at  the  center  of  the  stator  and  rotates  when  an 
electric  field  is  applied  at  a  frequency  between  the  radial  and  bending  resonance  modes. 
The  truncated  cone  shape  at  the  end  of  the  rotor  guarantees  a  permanent  contact  with  the 
tips  of  the  arms. 

The  operating  principle  of  this  motor  is  as  follows:  in  the  contraction  cycle  of  the 
stator,  the  four  arms  at  the  center  of  the  metal  endcap  clamp  the  rotor  and  push  it  in  the 
tangential  direction.  Since  the  radial  mode  frequency  of  the  stator  is  close  to  the  second 
bending  mode  frequency  of  the  arms,  the  respective  deformations  are  added  and  the  tips 
of  the  arms  bend  down.  In  the  expansion  cycle,  the  arms  release  the  rotor  from  a 
different  path  such  that  their  tips  describe  an  elliptical  trajectory  on  the  surface  of  the 

rotor. 

A.  Design  of  the  arms  at  the  inner  circumference  of  metal  ring 

The  length  of  the  inward  arms  was  chosen  such  that  the  resonance  frequency  of 
the  second  bending  mode  of  each  individual  arm  neighbors  the  resonance  frequency  of 
the  first  radial  mode  of  the  stator.  Using  this  strategy,  it  becomes  possible  to  produce 
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multi-mode  excitations  with  only  a  single  electrical  source  [12],  [13].  If  we  assume  that 
each  arm  behaves  like  a  simply  cantilevered  elastic  beam,  then  their  bending  frequencies 
can  be  estimated  according  to  the  following  expression  [9]: 


(1) 


where,  L,  tm,  p  and  E  represent  the  length,  thickness,  density  and  Young  Modulus  of  the 
beam,  respectively.  The  frequency  constants  Xn's  are  determined  from  the  following 
equation: 

CosX„CoshXn  =  -1  (2) 

The  first  three  values  of  the  X„2's  are  3.516,  22.03  and  61.7.  Using  these  values 
and,  the  dimensions,  and  mechanical  properties  of  the  arms,  the  first  three  bending  mode 
frequencies  of  the  arms  are  estimated  at  90.2, 494.1  and  1380  kHz.  The  second  bending 
mode  frequency  of  the  arms  is  close  to  the  stator's  radial  mode  resonance  frequency 

which  is  461  kHz. 


B.  Force  Transformation  Mechanism 

The  force  in  the  radial  direction  generated  by  the  stator  is  transferred  to  the 
tangential  direction  at  the  periphery  of  the  rotor  via  a  mechanism,  which  is  analogous  to  a 
slider-crank  mechanism;  assuming  that  the  inward  arms  and  the  rotor  are  rigid  (Fig.  2). 
When  the  stator  is  electrically  excited  to  generate  ultrasonic  vibrations  at  its  radial  mode 
resonance  frequency,  the  piezoelectric  and  metal  rings  expand  and  contract  in  the  radial 
direction.  The  stator  force  due  to  the  piezoelectric  effect  in  the  radial  direction  is 
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transferred  to  the  rotor  via  the  inward  arms.  Because  of  the  angular  orientation  of  the 


arms  (angle  cp  with  respect  to  the  radial  direction),  there  are  normal  and  tangential  force 
components  at  the  periphery  of  the  rotor.  The  normal  (N2)  and  tangential  (Ti)  force 
components  at  the  rotor  surface  are  calculated  as  a  function  of  the  angles  0  (the  angle 
between  the  origin  of  an  arm  and  the  contact  point  with  the  rotor),  (p  and  the  radial  force 


Fa. 


T2=  FAcos{(p)sm(9+(p) 

(3a) 

N2=  FAcos(<p)cos(8+<p) 

(3b) 

Finally,  the  normal  force  on  the  rotor  surface  is: 

JVj=  Facos(  p)cos(  &+  <p)  cos(  a ) 

(4) 

Since  the  length,  width  and  thickness  of  the  arms  as  well  as  the  radius  of  the  rotor,  which 
touches  tangentially  to  the  arms,  are  fixed  to  match  the  bending  frequency  of  the  arms  to 
the  radial  mode  resonance  frequency  of  the  stator,  the  angles  0  and  cp  were  decided 
automatically.  Small  cp  and  large  0  values  are  desirable  to  increase  tangential  force  T2, 
because  the  normal  force  between  stator  and  rotor  can  be  controlled  by  static  pressing 
force  Fn. 

The  mavimnm  torque  that  the  motor  can  generate  is  limited  by  the  static  pressing 
force,  and  is  provided  by  iVsina,  where  a  is  the  taper  angle  at  the  tips  of  the  arms 
ranging  from  10  to  30  degree  ( 17  degree  for  this  experiment).  The  stator  and  the  rotor 
contact  surfaces  are  two  mating  cones  that  are  kept  in  contact  by  the  normal  force.  They 
thus  constitute  a  conical  clutch  mechanism.  The  relation  between  the  normal  force  FN  and 
torque  capacity  T is  known  by  the  following  equation  [11]. 


T 


Ejl  rJL±Ii 
sin  a  2 


(5) 
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where,  ju,  ra,  and  r,  are  friction  coefficient,  outer  and  inner  radii  of  the  rotor  contact 
surface,  respectively.  The  above  torque,  which  is  also  called  blocking  torque,  is  the 
maximum  torque  that  this  motor  can  generate. 

Since  the  normal  force  and  the  geometrical  parameters  used  in  the  above 
expressions  are  all  known,  the  friction  constant  can  be  estimated  after  measuring  the 
stopping  time  of  the  motor,  that  is,  after  measuring  the  blocking  torque  of  the  motor. 

m.  FINITE  ELEMENT  ANALYSIS 

In  order  to  verify  the  conceptual  operation  principle  described  in  the  previous 
section,  the  ATILA  finite  element  code  was  used  extensively  to  analyze  the  behavior  of 
the  stator  and  the  vibration  transformation  mechanism  at  the  tips  of  the  arms.  Because  of 
the  asymmetrical  behavior  of  the  structure  and  the  angular  orientation  of  the  inward  arms, 
we  modeled  the  entire  stator.  Figure  3  shows  the  3D  mesh  of  the  free  stator.  The  total 
number  of  nodes  and  second  order  brick  elements  used  in  the  model  were  2036  and  256, 
respectively,  which  caused  computation  times  to  be  long  (20  min  per  frequency, 
approximately).  The  dimensions  and  material  properties  used  for  modeling  are  given  in 
Table  I  and  Table  II.  The  epoxy  layer  between  the  metal  and  piezoelectric  ring  was  not 
included  in  the  model. 

A.  Admittance  Spectrum 

The  computed  and  experimental  admittance  spectra  of  the  stator  are  compared  in 
Figures  4a  and  4b,  respectively.  This  figure  shows  that  the  three-dimensional  finite 
element  model  estimates  all  modes  generated  by  the  stator.  It  also  makes  it  possible  to 
identify  some  of  these  resonance  modes  (Fig.  4a).  For  instance,  the  highest  peak  at 
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477  kHz  is  the  radial  mode  of  the  stator,  whose  mode  shape  is  shown  in  Fig.  5a.  Since 
only  one  side  of  the  piezoelectric  ring  is  mechanically  loaded  (the  metal  ring  was  bonded 
to  one  side  of  the  piezoelectric  ring),  the  displacement  is  not  purely  in  plane.  The  stator 
takes  a  convex/concave  shape  as  it  contracts  and  expands  in  the  radial  direction.  Also,  a 
flexure  mode  is  generated  on  the  ring  because  the  arms  are  located  in  four  discrete 
locations  (Fig.  5b).  This  mode,  which  is  very  close  to  the  radial  mode  of  the  stator,  is 
excited  at  489  kHz  and  is  the  second  peak  in  Fig.  4a.  The  third  peak  on  the  admittance 
spectrum  corresponds  to  the  second  bending  mode  of  the  arms,  which  is  shown  in  Fig.  5  c. 
A  (4,l)-like  flexure  mode  is  also  excited  in  the  ring  at  the  bending  mode  frequency  of  the 
arms  which  bend  in  the  width  direction. 

Although  the  fundamental  radial  resonance  mode  was  estimated  within  less  than  5 
%  error,  magnitudes  of  the  motional  admittance  values  are  not  so  good.  Since  the  epoxy 
layer  between  the  piezoelectric  ring  and  the  metal  ring,  as  well  as  the  damping  caused  by 
the  cable  connections,  were  not  included  in  the  finite  element  model,  the  losses  in  the 
FEA  are  less  than  in  the  experimental  structure.  As  a  result,  the  calculated  motional 
admittance  values  are  larger  than  the  measured  ones.  Nevertheless,  the  magnitude  of 
calculated  off-resonance  admittance  matches  very  closely  to  the  measured  admittance. 

B.  Displacement  Transformation  Mechanism 

In  order  to  better  understand  the  mechanism  of  transformation  of  the 
displacements  that  takes  place  at  the  stator/rotor  interface  (from  the  radial  to  the 
tangential  direction),  displacements  versus  frequency  were  also  calculated  using 
harmonic  analyses  of  the  ATILA  code.  The  calculated  displacement  spectrum  at  nodes 
101  and  854  are  shown  in  Figures  6a  and  6b.  The  node  numbers  lOland  854  are  nodes 
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on  the  inner  circumference  of  the  metal  ring,  where  the  inward  arms  are  connected,  and 
on  the  tip  of  one  arm,  respectively  (Fig  3).  Figures  6a  and  6b  illustrate  the  displacement 
magnitudes  of  the  arms  as  a  function  of  frequency.  For  instance,  the  displacements  in  the 
z-axis  are  amplified  for  the  second  and  third  peaks  of  Figures  6a  and  6b.  This  is  due  to 
the  flexural  mode  of  the  stator  and  the  bending  modes  of  the  metal  arms.  Finally, 
elliptical  trajectories  at  nodes  101  and  854  were  calculated  for  three  different  frequencies 
(477, 490  and  500  kHz)  and  plotted  in  Fig.  7.  This  figure  also  demonstrates  that  the 
magnitudes  of  the  elliptical  trajectories  are  amplified  at  the  tips  of  the  arms. 

IV.  MOTOR  CHARACTERISTICS 

When  a  motor  size  is  decreased,  characterizing  the  motor  becomes  more  difficult. 
A  non-contact  characterization  method,  which  was  initially  proposed  by  Nakamura  [10], 
was  used  to  characterize  the  motor.  The  principle  of  this  method  consists  of  mounting  a 
load  (usually  a  disk  whose  moment  of  inertia  is  known)  onto  the  motor,  running  the 
motor,  and,  finally,  analyzing  the  transient  speed  obtained  as  a  function  of  time.  More 
explicitly,  the  angular  acceleration  of  the  motor  can  be  calculated  by  taking  derivative  of 
the  measured  speed.  The  transient  torque  is  then  calculated  by  multiplying  the  angular 
acceleration  by  the  moment  of  inertia  of  the  load,  which  is  known.  Using  this  method, 
the  starting  transient  response  of  the  motor  gives  the  speed-torque  relation.  Similarly,  the 
friction  coefficient  between  the  rotor  and  the  stator  is  estimated  from  the  transient 
response  for  stopping  the  motor. 

Instantaneous  output  mechanical  power  (Pout)  can  also  be  calculated  by 
multiplying  torque  (T)  with  angular  speed  (Q).  Input  electrical  power  (Pi„  )  is 
calculated  from  measured  input  voltage,  current  and  phase  angle  between  them.  From 
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input  and  output  powers  the  efficiency  can  be  calculated  by  dividing  output  mechanical 
power  to  input  electrical  power. 

A.  Estimation  of  Load  Characteristics 

The  load  consists  of  a  metal  disk  (lOg)  with  a  high  moment  of  inertia  (0.2 
kg*mm2)  compared  to  that  of  the  rotor.  It  is  mounted  onto  the  stator  (0.25g)  and  guided 
by  two  ball  bearings.  The  motor,  whose  frequency  spectrum  is  shown  in  Fig.  4b  (thick 
line),  was  driven  with  an  AC  voltage  of  20  V  at  468  kHz. 

The  position  of  the  loaded  motor  is  detected  in  a  form  of  square  pulses  through  an 
optical  sensor  (OMRON  EE-SPZ401 Y)  that  utilizes  a  photocell  pair.  The  frequency  of 
the  square  wave,  which  is  proportional  to  the  position,  was  converted  into  voltage  using  a 
frequency-to- voltage  converter.  Since  the  output  voltage  of  the  converter  is  proportional 
to  the  input  frequency,  the  transient  speed  of  the  motor  was  obtained  with  a  simple  gain 
factor.  Second  derivative  of  the  recorded  data  gives  the  angular  acceleration.  Finally,  the 
product  of  the  angular  acceleration  and  the  moment  of  inertia  of  the  rotating  disk  (I)  give 
the  motor  transient  torque. 

A  typical  transient  response  of  the  motor  under  loaded  condition  is  shown  in  Fig. 
8.  An  exponential  curve  fits  the  transient  speed  of  the  rotor: 

Q=Qo(l-exp(t/xr))  (6) 

where,  Qo  is  the  steady-state  speed  and  xr  is  the  time  constant,  which  can  be  found  by  trial 
and  error.  The  initial  value  of  xr  can  be  chosen  as  a  quarter  of  the  rise  time.  Then,  it  can 
be  adjusted  to  obtain  the  best  fit.  For  this  particular  measurement,  Qq  and  xr  were  found 
as  62  rad/sec  and  0.74  sec,  respectively. 
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A  blocking  torque  of  25  |xNm  was  obtained  from  the  stopping  time  of  the  motor. 
Substituting  25  pNm  in  the  torque  capacity  equation  (5),  gives  an  estimate  of  the  friction 
constant  as  0.15.  This  number  is  almost  half  of  the  expected  friction  constant  between 
brass  (stator)  and  steel  (rotor).  Finally,  the  load  characteristics  obtained  from  the 
transient  speed  are  shown  in  Fig.  9.  A  starting  torque  of  17  (iNm  is  one  order  of 
magnitude  higher  than  that  of  a  thin  film  motor  with  a  similar  size  [7]. 

The  difference  between  blocking  torque  (25  pNm)  and  starting  torque  (17  pNm) 
states  that  the  arms,  which  were  chosen  thin  enough  to  match  their  bending  mode 
frequency  to  the  stator  radial  mode  frequency,  are  not  strong  enough  to  push  the  rotor,  but 
they  deform.  On  the  other  hand,  an  efficiency  of  more  than  10  %  is  sufficient  for  real 
application  with  an  output  power  of  around  150  p  Watts. 

V.  CONCLUSION 

This  paper  presented  a  small  piezoelectric  ultrasonic  motor  whose  active  element 
is  a  bulk  piezoelectric  ring  3.0  mm  in  diameter  and  0.5  mm  in  thickness.  The  stator 
produces  a  displacement  transformation  such  that  multi-mode  resonance  frequencies  can 
be  excited  with  only  one  ac  source.  The  stator  was  analyzed  with  a  three-dimensional 
model  of  the  ATILA  finite  element  analysis  code,  and  its  dynamic  behavior  was 
predicted.  In  particular,  the  FEA  illustrated  the  bending  and  asymmetrical  flexure  modes 
that  are  generated  close  to  the  fundamental  radial  resonance  frequency  of  the  stator. 

The  rotation  in  this  motor  occurs  by  exciting  the  stator  with  a  combination  of  two 
modes:  the  fundamental  radial  mode  of  the  stator  and  the  second  bending  mode  of  the 
arms. 
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The  major  applications  of  this  motor  will  be  micromechanical  devices  for  the 


medical  and  microrobotic  areas. 
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TABLE  I 

Dimensions  of  the  stator 


Outer  Diameter  of  Piezoelectric  ring 

Rop 

3.0  mm 

Inner  Diameter  of  Piezoelectric  ring 

Rn» 

1.5  mm 

Thickness  of  Piezoelectric  ring 

tp 

0.5  mm 

Outer  Diameter  of  Metal  (Brass)  ring 

Rom 

3.0  mm 

Inner  Diameter  of  Metal  ring 

Rim 

2.4  mm 

Thickness  of  Metal  ring 

t  m 

0.15  mm 

Length  of  the  arm 

L 

1.02  mm 

Radius  of  circular  part  at  the  tip  of  the  arms 

r 

0.5  mm 

Width  of  the  arm 

Wm 

0.4  mm 

Angle  between  arm  and  radial  force  direction 

<P 

25° 

Angle  between  tips  of  the  arm  and  radial  force  direction 

0 

28.7° 

Tapering  angle  at  the  tips  of  the  arms 

a 

17° 

TABLE  II. 


Physical  and  piezoelectric  properties  of  ceramic  material 
_ (APC  International  Ltd.)  _ 


Parameter 

Symbol 

Units 

APC  841 

Relative  Permittivity 

eT33/eo 

1 

1380 

STii/80 

1 

1490 

Dielectric  Loss 

Tan  8 

% 

.35 

Coupling  Factors 

kp 

I 

.56 

k3i 

I 

.32 

k33 

I 

.65 

kis 

I 

.65 

kL 

I 

.48 

Piezoelectric  Strain  Constants 

-d3i 

10’^C/N 

109 

d33 

323 

Stiffness 

SE,i 

10'lzmz/N 

11.7 

Constants 

sE33 

17.3 

Mechanical  Quality  Factor 

Qm 

I 

1400 

Density 

-ft _ 

G/cc 

7.9 

Figure  Captions 


Fig.  1(a).  Encapsulated  motor  next  to  a  0.5  mm  metallic  pencil  (Steel  housing  was 
coated  with  an  epoxy  insulator),  (b)  Structure  of  the  motor  showing  the  number  of 
components. 

Fig.  2.  Top  and  cross-sectional  view  of  the  stator  vibrator.  The  metal  ring  is  bonded  to 
piezoelectric  ring  from  its  outer  circumference.  Interface  between  inwardly  directed  arms 
and  piezoelectric  ring  needs  to  be  empty.  Force  is  transferred  from  radial  to  tangential 
direction  through  the  arms. 

Fig.  3.  ATILA  3D  mesh  of  the  stator  vibrator  (256  brick  elements  were  used  in  the  model 
with  2036  nodes) 

Fig.  4.  (a)  calculated  admittance  spectrums  of  the  free  stator ,  b)  measured  admittance 
spectrums  of  the  free  stator  and  the  motor  (thick  line). 

Fig.  5  Top  view  of  calculated  mode  shapes  of  the  free  stator,  (a)  Mode  shapes  at 
fundamental  radial  resonance  frequency  at  477kHz.  Non-symmetric  structure  of  the  metal 
ring  distorts  radial  mode  shape  and  causes  the  stator  to  take  a  convex/concave  shape,  (b) 
A  flexure  mode  close  to  radial  mode  resonance  at  490kHz.  (c)  Mode  shape  corresponds 
to  bending  mode  resonance  frequency  of  the  arms  at  504kHz.  (d)  Another  flexural  mode 
at  569  kHz  of  the  stator  in  lateral  direction  corresponding  to  the  bending  mode  of  the 
arms  in  the  width  direction. 
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Fig.  6.  Magnitude  of  calculated  displacement  spectrums.  (a)  at  node  101,  and  (b)  at  node 
854  (for  location  of  nodes  101  and  854  see  Fig.  3.) 

Fig  7  Transformation  and  amplification  of  elliptical  trajectories  from  circumference  of 
the  stator  (node  #  101)  to  the  tips  of  the  arms  (node  #  854  )  at  frequencies:  477, 490  and 
500  kHz.  (for  location  of  node  101  and  854  see  Fig.  3.). 

Fig.  8.  Experiment  setup  to  characterize  small  size  motor  using  transient  response 
method 


Fig.  9.  Transient  response  of  the  motor  at  20V.  Transient  speed  was  curve  fitted  to  an 
exponential  function. 

Fig.  10.  Load  characteristics  of  the  motor. 
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Fig.  la  and  b.  Koc  et  al. 
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Using  a  two-dimensional  Ginzburg-Landau  model  and  the  finite-element  computational  method, 
we  have  calculated  stable  domain  configurations  resulting  from  a  second-order  ferroelectric  phase 
transition  for  a  finite-sized  system.  The  boundary  conditions  applied  here  correspond  to  fully  charge 
compensated  situations,  either  by  surface  electrodes  or  by  the  injection  of  charges  (or  defects)  near 
the  sample  surface.  The  domain  wall  thickness  of  a  finite  system  without  surface  electrodes  was 
found  to  become  thinner  as  it  approaches  sample  surfaces.  This  is  distinctively  different  from  that 
of  an  infinite  system  for  which  a  planar  wall  assumption  can  be  used.  The  orientation  of  the 
macroscopic  polarization  of  a  finite  system  without  surface  electrodes  was  found  to  be  determined 
by  its  aspect  ratio.  A  size  effect  was  observed  when  all  the  dimensions  were  reduced  simultaneously. 

The  relaxation  process  in  the  formation  of  domains  and  the  switching  process  have  also  been 
simulated  for  charge  neutral  boundary  conditions  using  a  time  dependent  Ginzburg-Landau  model. 

The  simulation  results  verified  that  the  surfaces  are  the  favored  nucleation  sites  for  domain 
switching.  ©  1999  American  Institute  of  Physics.  [S0021-8979(99)02522-0] 


I.  INTRODUCTION 

Many  physical  properties  of  technologically  important 
materials  are  determined  by  mesoscopic  structures,  such  as 
domains  and  domain  walls.  For  example,  it  has  been  deter¬ 
mined  experimentally  that  70%  of  the  dielectric  and  piezo¬ 
electric  contribution  in  ferroelectric  materials  can  be  attrib¬ 
uted  to  domain  related  activities.1’2  An  understanding  of  the 
formation  processes  of  such  mesoscopic  structures  is  there¬ 
fore  crucial  for  designing  and  engineering  better  active  ma¬ 
terials.  A  great  deal  of  effort  has  already  been  invested  in  this 
area,  including  first  principle  calculations,  analytical  ap¬ 
proaches,  and  computer  simulations,3-7  and  some  notable 
successes  have  been  achieved.  However,  analytical  solutions 
are  generally  limited  to  quasi-one  dimensional  (ID)  solu¬ 
tions,  and  most  simulations  have  used  periodic  boundary 
conditions  for  reasons  of  computational  efficiency.  While 
these  procedures  may  be  appropriate  when  dealing  with  very 
large  systems,  the  growing  trend  towards  making  smaller  and 
smaller  particles  and  thin  films  provides  new  challenges.  In 
particular,  size  effects  have  attracted  more  and  more  atten¬ 
tion  in  the  past  few  years  due  to  the  fast  growing  field  of  thin 
films  and  microelectromechanical  devices.  Static  analyses 
using  an  isotropic  surface  layer  for  particles  and  a  one¬ 
dimensional  model  for  thin  films  have  been  carried  out  by 

o  i  | 

several  authors  for  ferroelectrics. 

An  established  way  of  treating  finite  sized  materials  in 
simulation  studies  is  to  use  layered  vacuum— material— 
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vacuum  composites  and  periodic  boundary  conditions.  The 
periodic  boundary  conditions  allow  the  use  of  computational 
efficient  fast  Fourier  transform  techniques.  While  this  ap¬ 
proach  may  be  appropriate  to  approximate  free  boundary 
conditions  if  the  volume  of  the  vacuum  region  is  sufficiently 
large,  it  cannot  describe  practical  situations  in  which  external 
boundary  conditions  are  prescribed.  For  this  reason,  it  is  nec¬ 
essary  to  explore  models  that  allow  the  direct  application  of 
external  boundary  conditions  and  numerical  methods  that 
make  the  computation  size  manageable.  In  this  article  we 
have  used  scale  normalization  and  the  finite-element  method 
in  an  effort  to  simulate  domain  formation  in  a  finite  ferro¬ 
electric  system  with  high  computational  efficiency.  We  focus 
on  the  effects  of  surfaces  in  a  finite  system  under  prescribed 
boundary  conditions.  In  order  for  the  problem  to  be  well 
posed,  full  electric  charge  compensation  is  assumed,  which 
can  be  achieved  either  through  surface  electrodes  or  by  in¬ 
jection  of  charges  (or  the  distribution  of  charged  defects)  to 
neutralize  the  polarization  gradient  near  the  sample  surfaces. 

The  time  evolution  of  the  domain  pattern  under  specified 
boundary  conditions  was  also  investigated  using  a  time- 
dependent  Ginzburg-Landau  (TDGL)  theory.12  A  variable 
time-step  integration  technique  was  used  that  allows  the  en¬ 
ergy  minima  to  be  obtained  efficiently.  In  this  technique,  the 
time  step  was  increased  by  several  orders  of  magnitude  as 
the  system  approached  equilibrium.  Without  prescribed  po¬ 
larization  at  the  surfaces,  a  single  domain  state  was  found  to 
be  the  global  energy  minimum,  as  expected,  since  domain 
walls  contain  positive  energy  and  will  be  driven  out  of  the 
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system  without  surface  pinning.  Because  there  is  no  depolar¬ 
ization  energy  for  the  neutralized  system,  multidomain  states 
can  exist  only  with  prescribed  polarization  at  the  boundaries. 
For  the  case  of  charge  injection  near  the  surface  region,  po¬ 
larization  became  geometry  dependent  and  the  aspect  ratio 
was  found  to  break  the  orientational  degeneracy  of  the  sys¬ 
tem  due  to  the  gradient  energy.  Our  model  does  not  involve 
the  introduction  of  a  surface  layer;  the  surface  effect  is  in¬ 
trinsic  and  anisotropic. 

In  this  preliminary  exploration,  we  did  not  include  the 
strain  coupling  in  order  to  isolate  the  effects  of  electrical 
boundary  conditions.  The  model  system  has  a  square  to  rect¬ 
angular  phase  transition  [which  may  also  be  used  to  describe 
the  tetragonal  to  orthorhombic  phase  transition  in  a  three- 
dimensional  (3D)  system  with  the  tetragonal  dimension 
fixed],  and  a  square  to  rhombic  transition. 

An  outline  of  this  article  is  as  follows.  In  Sec.  II  we 
introduce  the  model,  in  Sec.  HI  we  discuss  the  effects  of 
boundary  conditions  on  the  formation  of  domains  in  equilib¬ 
rium  states,  and  in  Sec.  IV  we  discuss  the  time  evolution  of 
domain  processes. 


II.  GINZBURG-LANDAU  MODEL  FOR 
SQUARE-RECTANGULAR  AND  SQUARE-RHOMBIC 
FERROELECTRIC  PHASE  TRANSITIONS 


The  model  system  we  investigated  is  a  two-dimensional 
(2D)  analog  of  the  perovskite  structure  with  the  positive  and 
negative  ions  sitting  at  the  center  and  comer  of  a  square  unit 
cell,  respectively.  The  symmetry  of  the  high  temperature 
paraelectric  phase  is  square.  We  assume  that  the  system  un¬ 
dergoes  a  proper  second-order  ferroelectric  phase  transition 
so  that  the  system  can  be  described  by  a  fourth-order 
Ginzburg -Landau  theory.  Using  the  polarization  as  the  order 
parameter,  the  free  energy  density  for  a  second-order  phase 
transition  can  be  written  as 


/=  ax{P2x  +  P\)  +  an(P]  +  P22)2+  anP2xP\ 


+  d{ 
+  dy 


dxx) 


fpl) 

dx2J 


fPA 

dx2) 


2 


(1) 


where  only  the  first  coefficient  ax  =  a0(T—  Tc)(a0>0)  is 
temperature  dependent;  all  the  other  coefficients  are  assumed 
to  be  independent  of  the  temperature.  Unlike  previously  pro¬ 
posed  models,  the  gradient  energy  in  Eq.  (1)  is  anisotropic. 

Minimizing  the  free  energy  density  described  by  Eq.  (1) 
gives  three  types  of  stable  homogeneous  solutions  as  listed 
below. 


A.  Solution  1 

For  temperature  T>Tc,(ax>  0)  a  paraelectric  solution 
exists  for  which 


PX  =  P2= 0.  (2) 

The  symmetry  of  this  paraelectric  state  is  square. 


B.  Solution  2 


For  temperature  T<Tc(ax< 0),  if  <ar11>0  and  ax2>0,  a 
rectangular  ferroelectric  phase  is  stable,  and  has  four  vari¬ 
ants:  (±PR, 0)  and  (0,±PR),  where  the  spontaneous  polar¬ 
ization  Pr  is  given  by 


(3) 


The  symmetry  of  this  low-temperature  state  is  rectangular. 


C.  Solution  3 

For  temperature  T<Tc(a i<0),  if  <2'n>0  and  -A axx 
<ax2< 0,  a  rhombic  ferroelectric  phase  is  stable.  It  also  has 
four  variants:  (±PRh,±PRh),  with  the  spontaneous  polar¬ 
ization  component,  PRh ,  given  by 

PRh=  ^ A  an  +  ax2'  ^ 

In  order  to  compute  an  inhomogeneous  static  solution 
below  Tc ,  one  needs  to  solve  the  Euler  equations,2,3’13  that  in 
this  case  are  two  coupled  second-order  nonlinear  partial  dif¬ 
ferential  equations  of  the  following  form: 

2^1^U1  +  ^2^2,21  +  2J3P  iv22 

=  2axPx  +  Aaxx{P'\  +  P\)Px  +  2anPxP\,  (5a) 

2^1^2,22+  d2P  1,12+2^3^2,11 

=2axP2+ A  axx(P2x  +  P\)P2+2aX2P  2P2X.  (5b) 

Under  a  quasi- ID  approximation,  these  equations  can  be  re¬ 
duced  to  two  coupled  ordinary  differential  equations  which 
could  be  addressed  by  the  method  given  in  Refs.  4  and  14.  In 
this  article,  however,  we  will  directly  solve  this  set  of  partial 
differential  equations  using  the  finite-element  method  (see 
the  Appendix)  with  several  prescribed  boundary  conditions. 

For  demonstration  purposes  we  only  study  the  case  of 
#12>0,  for  which  the  rectangular  phase  is  stable.  The  case  of 
^12<0  has  a  very  similar  physical  appearance  if  one  rotates 
the  coordinate  system  45°.  For  generality,  Eq.  (5)  is  con¬ 
verted  into  dimensionless  form  by  making  the  following  sub¬ 
stitutions: 


where  a(>0)  is  a  dimensionless  scaling  factor  included  here 
to  give  the  flexibility  to  adjust  the  strength  of  the  gradient 
energy  with  respect  to  that  of  the  Landau  energy.  In  addition, 
we  define  three  new  dimensionless  coefficients  specifying 
the  anisotropy  of  the  gradient  energy  and  the  coupling 
strength  between  different  polarization  components, 
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c— 


and 


(8b) 


d= 


&12 

2an' 


(8c) 


The  final  dimensionless  form  of  Eq.  (5)  is 

au  b  cu  vv=  -  u  +  u3  +  (  \  +  d)ui?‘,  (9a) 

i?  +  (l+d)ru2.  (9b) 

For  a  given  system,  the  gradient  coefficients  a,  b>  and  c 
may  be  determined  by  inelastic  neutron  scattering 
experiments,14,15  while  the  nonlinearity  coefficient,  d ,  can  be 
obtained  from  measurements  of  the  nonlinear  dielectric 
constant.16  In  the  current  2D  model  system,  we  will  investi¬ 
gate  a  particular  set  of  normalized  constants  to  provide  the 
general  trend  of  the  solutions.  For  any  specific  material  one 
could  apply  the  procedure  described  in  this  article  to  obtain 
the  inhomogeneous  solutions  after  determining  these  param¬ 
eters  from  experiments. 

The  coupled  partial  differential  equations,  Eqs.  (9a)  and 
(9b),  were  solved  using  the  finite-element  method  described 
in  the  Appendix.  Nine  node  quadrilateral  elements  were  used 
in  our  calculations  and  the  density  of  the  elements  in  any 
given  direction  reflects  the  real  dimension  in  that  direction. 
The  dimensions  of  the  system  were  normalized  to  1  so  that 
the  computational  procedure  could  be  applied  to  any  size 
through  a  scaling  factor.  The  scaling  factor  was  also  used  to 
convert  the  simulation  results  back  to  the  original  dimen¬ 
sions  after  the  problem  was  solved  in  the  normalized  system. 


III.  EFFECTS  OF  BOUNDARY  CONDITIONS  ON  THE 
FORMATION  OF  DOMAIN  STATES 

A.  Surface  electrodes  versus  charge  injection 

The  natural  boundary  conditions  arising  in  the  finite- 
element  solution  technique  correspond  to  a  system  with  sur¬ 
face  electrodes  that  is  short  circuited  during  the  phase  tran¬ 
sition.  The  depolarization  field  is  screened  by  the  free  charge 
so  that  the  magnitude  of  the  polarization  can  be  maintained 
as  a  constant  throughout  the  sample.  Under  this  circum¬ 
stance,  a  monodomain  state  is  the  global  energy  minimum.17 
Our  model  simulation  quickly  converges  to  such  a  single 
domain  state.  The  initial  guess  used  when  solving  Eq.  (9) 
was  generated  by  the  time-dependent  model  described  in 
Sec.  IV. 

The  boundary  condition  P*n=0  corresponds  to  the  situ¬ 
ation  of  charge  injection  (or  charged  defect  concentration) 
near  the  surface  region.  This  case  reflects  many  realistic  situ¬ 
ations  in  which  a  surface  layer  develops  with  reduced 
polarization.8-11  We  found  a  strong  aspect  ratio  effect  when 
this  boundary  condition  was  applied.  The  polarization  of  the 
stable  configuration  was  always  parallel  to  the  long  dimen¬ 
sion.  In  other  words,  the  dimensional  difference  breaks  the 
symmetry  degeneracy  of  the  two  pairs  of  solutions  under  this 
boundary  condition.  Only  at  the  critical  aspect  ratio  of  1  do 
all  four  possible  solutions  (with  orientation  (10))  have  equal 
chances  of  forming. 


FIG.  1.  Reduced  effective  polarization  near  surfaces  causes  a  size  effect. 
There  is  a  critical  size  at  which  the  effective  polarization  vanishes.  The 
parameters  used  in  the  calculations  are  a~  b  —  c~  0.05  and  d=3. 


The  implication  of  this  result  is  that  the  natural  state  for 
a  free  standing  thin  film  without  surface  electrodes  should 
have  polarization  lying  in  the  plane,  although  the  net  polar¬ 
ization  would  still  be  zero  due  to  the  in-plane  random  orien¬ 
tation.  We  must  emphasize  that  this  conclusion  is  only  true 
for  a  free  standing  film,  otherwise,  the  polarization  orienta¬ 
tion  will  be  determined  by  the  constraints  originating  from 
lattice  constant  mismatch  between  the  film  and  the  substrate, 
which  is  beyond  the  scope  of  this  model  system. 

B.  Size  effects  in  a  single  domain  system  with  P.n=0 
boundary  conditions 

Size  effects  become  increasingly  important  as  material 
systems  become  smaller  and  smaller.  Many  experimental  re¬ 
sults  showed  that  material  properties  degrade  as  the  size 
decreases.18-23  To  determine  this  size  effect,  we  used  a 
square  system  of  size  LXL  and  continuously  reduced  the 
value  of  L.  The  net  effective  polarization  is  shown  in  Fig.  1. 
A  critical  size  is  clearly  shown  in  the  computational  results 
although  the  actual  value  of  the  critical  size  depends  on  ma¬ 
terial  parameters,  particularly  on  the  gradient  coefficients. 
Our  results  are  similar  to  those  reported  earlier  by  introduc¬ 
ing  an  isotropic  surface  layer  with  reduced  polarization.11 
The  difference  is  that  we  did  not  make  the  assumption  of  a 
surface  layer,  but,  rather,  incorporated  the  gradient  anisot¬ 
ropy  and  boundary  conditions.  Our  model  predicts  strong 
size  effects  only  when  all  the  dimensions  are  being  reduced 
at  the  same  time.  If  the  reduction  is  only  in  one  of  the  di¬ 
mensions,  such  as  in  thin  films,  the  polarization  can  still  be 
formed  in  the  other  larger  dimension  unless  a  reduced  polar¬ 
ization  layer  is  introduced. 

Generally  speaking,  this  critical  size  is  directly  related  to 
the  strength  of  the  nonlocal  interaction  of  the  polarization, 
i.e.,  the  amplitude  of  the  polarization  gradient  coefficients. 

C.  Twin  structures  and  their  size  dependence 

Under  the  quasi- ID  approximation,  a  single  twin  solu¬ 
tion  can  be  derived  using  the  method  described  in  Refs.  4 
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FIG.  2.  Equilibrium  configuration  of  a  180°  twin  with  top  and  bottom  sur¬ 
faces  electroded  and  the  left  and  right  boundaries  being  prescribed  with  the 
boundary  condition  of  ±  1  for  the  vertical  component  and  zero  for  the  hori¬ 
zontal  component,  respectively. 

and  14  for  infinite  systems.  Here,  we  study  a  finite  system 
with  surface  electrodes  at  the  top  and  the  bottom  surfaces 
and  prescribe  the  polarization  on  the  two  side  surfaces.  The 
180°  twin  and  90°  twin  solutions  are  shown  in  Figs.  2  and  3, 
respectively.  In  order  to  compare  them  with  the  quasi- ID 
solutions  of  Ref.  4,  we  look  at  the  solution  along  the  middle 
of  the  vertical  dimension  in  the  2D  structure  (maximum  do¬ 
main  wall  width).  The  contour  lines  in  Figs.  2  and  3  repre¬ 
sent  the  position  where  the  polarization  amplitude  has 
reached  79%  [  =  tanh(l)]  of  its  maximum  value.  We  take  the 
region  in  between  the  two  contour  lines  as  the  “domain 
wall.”  One  can  see  that  the  domain  wall  thickness  for  the 
90°  twin  decreases  as  it  approaches  the  top  and  bottom  sur¬ 
faces.  The  corresponding  parameters  used  for  the  simulation 
are  given  in  the  Figs.  2  and  3  captions. 

The  results  imply  that  the  quasi- ID  assumption  (plane 
wall  assumption)  is  inadequate  if  the  systems  is  very  small. 
As  shown  in  Fig.  3,  the  symmetry  inside  a  90°  wall  is  lower 


FIG.  3.  2D  system  with  the  dimensions  parallel  to  [11]  and  [11].  The 
equilibrium  configuration  of  a  90°  twin  is  obtained  with  the  top  and  bottom 
surfaces  electroded  and  the  left  and  right  boundaries  prescribed  with  the 
boundary  condition  of  the  polarization  pointing  to  [10]  and  [01]  directions, 
respectively  (45°  with  respect  to  the  surface  orientations).  The  contour  lines 
illustrate  the  change  of  domain  wall  width  from  the  center  to  the  top  and 
bottom  surfaces. 


FIG.  4.  Dependence  of  the  domain  wall  thickness  on  the  amplitude  of  the 
gradient  coefficients.  We  chose  a  =  b=c  for  convenience  to  compare  these 
results  with  that  of  ID  models.  The  dashed  line  is  the  square-root  relation¬ 
ship  obtained  from  a  ID  model. 

than  that  of  the  single  domain  regions,  and  the  direction  of 
the  polarization  at  the  center  of  the  domain  wall  is  along  the 
diagonal  direction  of  the  unit  cell  (i.e.,  45°  from  the  single 
domain  regions).  In  Fig.  4  we  plot  the  domain  wall  thickness 
measured  at  the  center  of  the  simulated  systems  versus  the 
gradient  coefficient  (for  convenience,  we  have  set  a  a  =  b 
=  c)  while  keeping  all  the  other  parameters  fixed.  The 
dashed  line  is  the  square  root  relationship  predicted  by  the 
ID  model,  and  the  dots  are  our  computational  results.  One 
can  see  that  the  simulation  results  agree  well  with  the  ID 
analysis  even  when  the  wall  thickness  reached  50%  of  the 
system  size.  The  deviation  for  the  large  gradient  coefficient 
is  caused  by  the  boundary  constraints  to  the  domain  wall. 

IV.  RELAXATION  IN  DOMAIN  FORMATION  AND  FIELD 
INDUCED  SWITCHING  PROCESS 

A  time-dependent  Ginzburg-Landau  model  can  provide 
information  about  the  time  evolution  of  domain  formation 
from  a  quenched  state  below  Tc .  The  TDGL  used  here  is  a 
relaxation  equation  for  a  nonconserved  order  parameter,  cor¬ 
responding  to  model  A  described  by  Hohenberg  and 
Halperin.12  Simulations  using  periodic  boundary  conditions 
have  revealed  that  the  domain  coarsening  slows  down  con¬ 
siderably  upon  approaching  equilibrium.6,7  In  order  to  simu¬ 
late  a  finite  system  in  direct  space  within  a  realizable  time 
scale  without  having  to  invoke  the  periodic  boundary  condi¬ 
tion,  we  used  a  variable  time-step  scheme  that  is  described  in 
the  Appendix,  The  computational  efficiency  of  this  method 
allows  us  to  perform  real  space  finite  size  simulations  and  to 
explore  the  effects  of  boundary  conditions. 

The  master  equation  used  here  is  of  the  same  form  as 
model  A  described  in  Ref.  12: 

dPi  SF 

^=-r  —  +R(r9t)  +  TEg9  (i  =  l,2),  (10) 

where  T  is  the  kinetic  coefficient,  F  is  the  volume  integration 
of  Eq.  (1),  R(r,t)  represents  the  Gaussian  random  noise,  and 
Ej  is  the  external  electric  field. 
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When  carrying  out  the  simulation,  R(r,t)  [<1CT4],  was 
only  used  for  initiating  the  nucleation  process  at  the  begin¬ 
ning  from  a  paraelectric  state,  and  was  neglected  at  later 
stages  because  it  does  not  affect  the  ultimate  equilibrium 
states.  Using  a  normalized  time  rand  a  frequency  fl  defined 
by 

r=  -2  tfjFr,  (11) 


CO 

-2axY' 


(12) 


and  substituting  the  dimensionless  variables  and  coefficients 
defined  by  Eqs.  (6)-(8),  we  can  derive  the  normalized  gov¬ 
erning  equations  for  the  temporal  evolution  of  the  normal¬ 
ized  polarization  under  an  external  AC  electric  field. 

-e^cos(ftr),  (13a) 

v,r=ai',tf+buJfv+ci',vv+  v-  i?-(l+d)w2 

-evcos(Yl  r),  (13b) 

where  e  =  E/(-2alPR)  is  the  amplitude  of  the  dimension¬ 
less  form  of  the  external  electric  field. 

Figure  5  shows  the  domain  formation  process  with  P.n 
=  0  boundary  conditions  from  a  quenched  state  below  T c . 
The  simulation  was  initiated  from  a  small  amplitude  random 
disturbance.  Because  the  amplitude  is  too  small  to  see  at  the 
beginning,  we  have  multiplied  the  magnitude  of  the  local 
polarization  vectors  by  an  amplification  factor.  For  example, 
the  amplitude  of  the  polarization  vectors  drawn  in  Fig.  5(d) 
has  been  amplified  80  000  times.  It  is  interesting  to  see  that 
the  amplitude  of  local  polarization  first  shrank  to  almost  zero 
before  re-nucleating  into  small  clusters  which  grew  into 
many  small  domains.  As  the  coarsening  progressed,  the 
dominant  domains  grew  at  the  expense  of  the  smaller  ones. 
The  final  equilibrium  state  was  determined  by  the  nucleus 
whose  local  polarization  amplitude  reached  the  equilibrium 
value  first.  In  order  to  compare  the  efficiency  of  the  variable 
time  step  finite-element  method  used  in  our  simulation,  we 
define  a  unified  dimensionless  time  — 10  3  as  the  unit. 
Using  this  unit  as  a  measure,  the  polarization  evolution  from 
Figs.  5(a)-5(f)  will  take  42  000  time  steps.  While  using  the 
variable  time  step  method,  it  only  took  140  variable  time 
steps  to  reach  the  equilibrium  state  in  Fig.  5(f),  and  this 
represents  a  tremendous  saving  of  computational  time.  The 
time  steps  of  both  measures  are  given  in  the  Fig.  5  caption. 

Figure  6  is  another  simulation  sequence  with  surface 
electrode  boundary  conditions.  In  this  case,  the  equilibrium 
was  reached  in  14  000  time  steps  in  the  unit  of  A  t^,  but  our 
simulation  only  required  110  variable  time  steps.  It  is  impor¬ 
tant  to  point  out  that  our  simulation  showed  no  size  effect 
with  the  electroded  surfaces  in  terms  of  the  polarization  for 
this  case.  This  is  the  main  difference  between  our  model  and 
previous  models,  which  assumed  a  surface  layer  with  re¬ 
duced  polarization. 

Letting  fl->0  in  Eq.  (13)  we  can  simulate  the  quasistatic 
switching  process.  It  is  much  easier  to  form  nuclei  of  the 
opposite  domain  near  the  surface  region  because  of  the  in- 
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FIG.  5.  Time  evolution  of  the  domain  formation  process  from  a  random  to 
a  single  domain  state  for  the  P.n=0  boundary  condition,  (a)  r=  0  and  the 
amplitude  of  the  local  polarization  has  been  amplified  3333  times;  (b)  30 
variable  time  steps=35Af0 ,  and  the  amplitude  of  the  local  polarization  has 
been  amplified  20  000  times;  (c)  35  variable  time  steps=97A/0  and  the 
amplitude  of  the  local  polarization  has  been  amplified  20  000  times;  (d)  40 
variable  time  steps  =  350Ar0,  and  the  amplitude  of  the  local  polarization  has 
been  amplified  80  000  times;  (e)  45  variable  time  steps=  1600A/0,  and  the 
amplitude  of  the  local  polarization  has  been  amplified  6667  times;  (f)  140 
variable  time  steps =42  000 A  r0  and  there  was  no  amplification  at  the  local 
polarization  vectors.  The  parameters  used  in  the  simulation  are  a  =  0.045, 
6  =  0.055,  c  =  0.05,  d—2.5,  and  e=0. 


fluence  of  the  boundary  conditions.  In  our  model  the  switch¬ 
ing  does  not  happen  homogeneously  in  the  whole  system 
but,  rather,  nucleates  at  the  surfaces  and  then  progresses  via 
the  movement  of  domain  walls.  The  time  evolution  of  the 
domain  switching  process  for  a  90°  twin  structure  is  illus¬ 
trated  in  Fig.  7.  The  net  polarization  of  the  twin  was  pointing 
to  the  right  before  application  of  an  electric  field,  which  is 
opposite  to  the  initial  polarization  direction.  Two  domain 
walls  formed  near  the  left  and  the  right  surface  regions  at  the 
beginning,  then  they  moved  toward  each  other  and  finally 
annihilated  one  another  in  the  middle  of  the  structure,  form¬ 
ing  a  reversed  twin  with  the  effective  polarization  pointing  to 
the  left  (the  field  direction).  It  is  interesting  to  see  that  a  swirl 


5744  J.  Appi.  Phys.,  Vol.  86,  No.  10,  15  November  1999 


Cao,  Tavener,  and  Xie 


i  \ \w \^^-~/ / / / /  7  ~ 

►  \  W  /  » 

•  s / /  /  /  t  ' 

♦  \ \\ /  /  /  S  /  /  ' 

♦  \\ /  T 

*  /  • 

»  \\\\\^^/////  < 

*  S  /  X  /  /  • 

*  \ /  • 

♦  WWW'"//////  t 

*  \\ 'aW^^'S/X  /  /  /  t 

•  \ • 

*  \ \v\ /  ♦ 

*  WWW'"//////  t 

*  WW\^'"/////  /  • 
i  \  wwww////  /  " 


(a) 


0>) 


(a) 


(b) 


■*  *  *  • 

\  *  >  • 

<  .  *  V  ^  V  4  '  ‘ 

*  *  *  »  4  *  *  '  ' 

1  ....  . 

*******  1  J 

*  *  *  *  »  /  i  *  X 

*  *  ^  V  *  i  i  V 

. . 

*  '  >»  V  W  *  '  - 

A 

T 

TVT 

J* 

TTT 

T 

V 

^7 

^zryTTT 

* 

T 

\  K  Vx 

V  - 

✓ 

t 

l  X  V 

-*■ 

/  t  t 

t 

N 

••  a 

HH 

/  /  / 

< 

X  X  x  X 

y  y 

V  - 

i 

* 

*  X  \ 

/  /  t 

t 

V 

-  jr 

v  a  y  y  / 

' 

X  X  X  X 

—  ^  y 

V  - 

* 

* 

*  X  \ 

/  t  t 

f 

X 

--  A 

■ 

\"///  y 

-- 

X  X  X  x 

~y  y 

V  - 

4 

x 

\  x  X 

-*r 

x  t  t 

f 

* 

-  y 

—  y  y  y  / 

x  X  X  x 

y 

V  - 

4 

\ 

XXX 

/ ft 

t 

\ 

-  y 

/  /  /  A 

-* 

' 

x  x  X  x 

■+-  y  y 

\  - 

4 

\ 

X  X  X 

"ft 

t 

* 

-  y 

Vx  -  ✓  /  /  / 

< 

v  X  X  V 

•*-  y  y 

V  * 

/ 

V 

X  X  X 

-♦ 

*  /  / 

t 

* 

-  y 

■ 

W-  ^  /  /  * 

* 

'  X  X  X 

~  y  y 

\  ' 

* 

x 

X  Xx 

j* 

y  1 1 

t 

-  y 

■ 

-  -r  /  i  * 

* 

* 

*  \  X  X 

-  *  y 

V  ' 

4 

\ 

XXX 

* 

/ft 

t 

* 

'  y 

\Xw  /  f  * 

~  y  y 

V  ' 

4 

x 

XXX 

.*■ 

*  f  t 

t 

* 

-  y 

■ 

\  X  ^  '  i  1  v 

* 

/  t  *  ' 

*  y  y 

\ ' 

4 

* 

X  X  x 

jr 

/ 1  t 

t 

*  y 

\  X  '  *  *  *  \ 

■*> 

/  t  *  ' 

*  y  y 

\  V 

- 

* 

XXX 

J" 

/  f  t 

t 

- 

*  y 

x  -  •  x  x 

/  /  *  - 

<  y  y 

\  ' 

x 

xxx 

/ft 

4 

- 

>  / 

XX*  'Xx 

-- 

A  /  '  • 

/  /  / 

\  * 

XXX 

•* 

/ft 

/ 

/  / 

\  X  \  '  x  x 

/  /  *  * 

/  /  / 

UL 

xxx 

/// 

> 

* 

±jL 

v  x  ♦  *-  '  X  X  ^ 

/  A  -  V 

*  // 

(C) 


(d) 


(c) 


(d) 


\ 

X 

X 

X 

\ 

X 

X 

X 

X  X 

X  X  X  X  4  t 

\ 

X 

X 

X 

X 

X 

X 

X 

X  X 

X  X  X  X  X  X 

X 

X 

X 

X 

X 

X 

X 

X 

X  X 

X  X  X  X  X  X 

X 

X 

X 

X 

X 

X 

X 

X 

X  X 

X  X  X  X  X  X 

X 

X 

X 

X 

X 

X 

X 

X 

X  X 

X  X  X  X  X  X 

X 

X 

V 

X 

X 

X 

X 

X 

X  X 

X  X  X  X  X  X 

V 

X 

X 

X 

X 

X 

X 

X 

X  X 

X  X  X  X  X  X 

X 

X 

X 

X 

X 

X 

X 

X 

X  X 

X  X  X  X  X  X 

V 

X 

X 

X 

X 

X 

X 

X 

X  X 

X  X  X  X  X  X 

V 

V 

X 

X 

X 

X 

X 

X 

X  X 

X  X  X  X  X  X 

X 

X 

X 

X 

X 

X 

X 

X 

X  X 

X  X  X  X  X  X 

X 

X 

X 

X 

X 

X 

X 

X 

X  X 

X  X  X  X  X  X 

« 

X 

X 

X 

X 

X 

X 

X 

X  X  X  X  X  X  X  t 

X 

X 

X 

X 

X 

X 

X 

X  X 

X  X  X  X  X  X 

X 

X 

X 

X 

X 

X  X  X  X  X  x ;  x  u 

1 

X 

X 

X 

X 

\AJl±±±±±±A . k) 

(•) 


(0 


^VV^T/  -  ^  V  V'"/// 

VW"*  — 

VWX— — 

\V\\rx  — ^*./'///// 

(e) 


[VAVAVVv-^aVTa'TTI 

VWWW'-'/////// 
^  X  X  X  X  X  X  X 
WWWV'"/////// 
wxxwv^-'vyAA'A'A'y 

XWWXV - /////// 

WWVW^-y////// 

W\V\V^--^>'yVA/A'A/ 

VWWWW////// 

WWVW XXX  X X X X 

VWWW'"/////// 

(0 


FIG.  6.  Time  evolution  of  the  domain  formation  process  from  a  random  to 
a  single  domain  state  for  a  system  with  shortened  surface  electrodes  (a) 
t=  0  and  the  amplitude  of  the  local  polarization  has  been  amplified  4000 
times;  (b)  35  variable  time  steps=78Ar0  and  the  amplitude  of  the  local 
polarization  has  been  amplified  24  000  times;  (c)  40  variable  time  steps 
=  230Ar0  and  the  amplitude  of  the  local  polarization  has  been  amplified 
40  000  times;  (d)  45  variable  time  steps=780Af0.  and  the  amplitude  of  the 
local  polarization  has  been  amplified  40  000  times;  (e)  50  variable  time 
steps =3300Af0  and  the  amplitude  of  the  local  polarization  has  been  ampli¬ 
fied  600  times;  (f)  110  variable  time  steps=  14000Ar0,  and  there  was  no 
amplification  of  the  local  polarization  vectors.  The  parameters  used  in  the 
simulation  a =0.06,  6  =  0.04,  0.05,  d= 2,  and  e=0. 

pattern  developed  in  Fig.  7(d).  Such  a  phenomenon  can  oc¬ 
cur  because  our  current  model  does  not  include  elastic  strain 
coupling.  The  walls  are  expected  to  be  planar  while  moving 
under  the  force  of  the  electric  field  if  the  elastic  constraints 
are  included. 

V.  SUMMARY  AND  CONCLUSIONS 

A  2D  ferroelectric  model  system  has  been  investigated 
using  a  finite-element  technique  combined  with  a  variable 
time  step  numerical  scheme.  The  unique  feature  of  our 
model  is  the  introduction  of  the  anisotropic  gradient  energy 
which  is  compatible  with  the  structural  symmetry.  We  also 
directly  solved  the  coupled  partial  nonlinear  differential 
equations  rather  than  trying  to  reduce  the  problem  to  one 


FIG.  7.  Switching  process  of  a  90°  twin.  The  starting  net  polarization  was 
pointing  to  the  right.  The  time  steps  are  (a)  0;  (b)  80  variable  time  steps 
=  170A/0;  (c)  140  variable  time  steps  =6200Ar0;  (d)  180  variable  time 
steps=17  000A/o;  (e)  280  variable  time  steps=30  000A/0;  (f)  330  vari¬ 
able  time  steps  =  37  000A/C.  The  parameters  used  for  the  simulations  are 
a  =  0.014,  6=0.01,  c  =  0.012,  d~ 2.0,  and  field  amplitude  e~  10~2.  The 
switching  process  started  at  both  surfaces  on  the  left  and  right  sides  and 
propagates  to  the  interior.  Two  domains  walls  form  near  the  surfaces  and 
they  move  toward  each  other  until  annihilation. 

dimension.  Due  to  the  high  efficiency  of  the  numerical 
method,  we  were  able  to  simulate  both  the  equilibrium  po¬ 
larization  configuration  and  the  temporal  evolution  of  the 
local  polarization  formation  in  a  finite  system  with  different 
boundary  conditions.  The  simulation  results  agree  well  with 
all  the  known  facts  for  a  second-order  ferroelectric  phase 
transition,  such  as  size  effects,  switching  via  surface 
nucleation-domain  wall  movement,  and  aspect  ratio  effects. 
In  particular,  we  have  demonstrated  the  temporal  evolution 
of  domain  formation  and  the  switching  process  in  a  90°  twin 
system. 
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APPENDIX 

We  used  the  finite-element  code  ENTWIFE24  to  solve 
the  following  coupled  system  of  partial  differential  equations 
[see  Eq.  (13)]  in  a  rectangular  domain  D : 

u, T=au,tf+bpyfv+cuiVV+u-u3-(l  +d)u?? 

e  ^cos(fl  7*),  (Ala) 

v, T=av,(f->rbu,fv+ci',vv+  v-  S~(\  +d)m2 

-e^cos^r).  (Alb) 

The  code  ENTWIFE  incorporates  sophisticated  parameter 
continuation  capabilities  and  the  ability  to  locate  bifurcation 
points  at  which  qualitative  changes  in  the  solution  set  occurs, 
making  it  particularly  suitable  to  solve  multiparameter  prob¬ 
lems.  The  numerical  bifurcation  techniques  implemented  in 
the  code  have  proven  to  be  useful  to  study  a  wide  range  of 
problems  in  fluid  and  other  continuous  systems. 

In  order  to  write  Eqs.  (Ala)  and  (Alb)  in  a  form  to 
which  we  can  apply  the  standard  finite-element  approach 
(see,  e.g.,  Refs.  25  and  26),  it  is  convenient  to  introduce  the 
notation, 

p  f.  *2  v . 

so  that 

Pu~u,{  311(1  Pl.i-f'.v'  etc- 

First,  consider  the  steady  state  problem  without  an  exter¬ 
nal  field.  Let  f(p)  be  a  vector  valued  function, 

fi(p)  =  apijj+/3pj,ij+  TPk.ki  (*  =  1.2),  (A2) 

where  A:  =  3  —  /  with  a  summation  over  the  repeated  index  j 
=  1,2  (but  not  over  k).  Expanding  Eq.  (A2)  using  the  origi¬ 
nal  variables  we  have 

f(  \=(  *(“•#+«•  we*] 

W  \  V ,  f(+  V,  vv) + 0{u ,  fv+  V,  vv) +  yu , £v) 

_  l  (a+fi)u,  &+(/?+  y)v,£r>+  au,vv\ 

\  on/,(f¥{0+  y)u,£v+{a+0)  ^'yyl 

Defining  a  second  vector-valued  function  g(p)  by 

\  p—  — (1  +  d)u2rj9 

the  static  problem  without  external  field  becomes 

f(p)+g(p)=°.  (A3) 

provided  constants  a,  /?,  and  y  are  related  to  constants  a,  b, 
and  c  as 

a=c,  p=a—c ,  y=b-c. 

Let  the  boundary  of  domain  D  to  be  divided  into  two 
parts,  JD  =  <?D  j  +  <?D2 ,  where  u  and  ^are  specified  on  3D\ , 
i.e.,  p=q  on  SD\ ,  and  are  not  specified  on  dD2 . 

We  seek  a  solution  peV?  of  Eq.  (A3)  where  Vq= 
{p:peH1(L>),p=qon^iD1}.  Here  H'(£>)  is  the  space  of 
vector  functions  defined  on  D  whose  function  values  and  first 
derivatives  lie  in  L2(£>).  In  the  usual  way,  we  construct  the 


weak  form  of  Eq.  (A3)  by  taking  the  dot  product  of  (f+g) 
with  a  (vector)  test  function  r,  and  requiring  that  the  integral 
of  the  dot  product  over  the  domain  D  vanishes  for  all  test 
functions  reV0«  where 

V0e{P’Ps//1(L,i):P=0  on  }• 

Summing  over  the  repeated  indexes  i  and  j,  notice  that 

"TSTj^ + ^  bun) + (Pk,krj) 

=  apijjri+/3pj'ijri+  ypkMr,+  aptJrtJ+ PpUirW) 

+  yPk,krj,i>  (A4) 

where  k  =  3~j. 

By  applying  the  divergence  theorem  in  the  plane, 

J  Jj.a'Pi,jjr i  +  fiPj,ijr ;+  YPk,kjrj)d£d£ 

=  (apijrinj+j3pjjrinj+  yp^n^ds 
J  #D 

J  ^PiJ^iJ^"  PP yPk,kT j,j)d£d£,  (A5) 

where  n  is  the  outward  surface  normal.  By  including  the 
terms  in  g(p),  the  weak  form  of  the  steady  state  equations  is 

|  ^[P\-p\-^+d)pxp\\rx+\_p1-p\-{\+d)p2p\]r1 

-[<zPijrij+J3Pj,irij+  7Pk,kT j,j]}dfdf=0.  (A6) 

Since  the  test  functions  r,  ,  z  =  1,2  vanish  on  dDx  by  con¬ 
struction,  the  “natural”  boundary  conditions  that  apply  on 
<?D2  arise  from  the  requirement  that  the  boundary  integral 
term  in  Eq.  (A5)  vanishes  for  all  possible  test  functions.  The 
natural  boundary  conditions  on  <9D2  are,  therefore, 

<xPijnj+flpJ,inl+  ?ttiP2.2=0,  (A7a) 

aP2Jnj  +  fiPj.2nj  +  1.1  =  (A7b) 

Rectangular  elements  with  biquadratic  interpolation  were 
used  to  construct  finite-dimensional  subspaces  VqCVq  and 
VqCVq.  The  finite-element  method  seeks  to  determine  a 
(finite-dimensional)  function  ph  e  Vhq  that  satisfies  Eq.  (A6) 
for  all  rA  e  Vq  . 

A  Gaussian  quadrature  was  used  to  approximate  the 
double  integral  in  Eq.  (A6)  and  the  resulting  nonlinear  sys¬ 
tem  of  equations  for  the  nodal  degrees  of  freedom  was 
solved  via  Newton’s  method.  A  direct  method  was  used  to 
solve  the  linear  system  at  each  Newton  iteration. 

The  time-dependent  problem  in  Eqs.  (Al)  can  be  written 
as  an  initial  value  problem  of  the  following  form: 

^  =  h[p(r)],  t  e  (0,7)  with  p(0)  =  p0.  (A8) 
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In  our  simulation,  the  method  of  Byrne  and  Hindmarsh27 
was  used,  which  was  designed  for  “stiff 5  initial  value  prob¬ 
lems  of  the  form 

~F[Y(f),/],  t e (0,T)  with  Y(0)  =  Y0,  (A9) 

where  F :RNXRl 2-+RN  is  continuous  and  satisfies  a  Lipschitz 
condition.  Solutions  to  stiff  differential  equations  develop 
over  a  wide  range  of  time  scales.  Unlike  fixed  time-step  in¬ 
tegration  techniques,  the  method  of  Byrne  and  Hindmarsh  is 
a  variable  time-step  method  based  on  backward  difference 
formulas  of  the  order  1-5.  In  such  variable  time-step  meth¬ 
ods,  the  interpolation  knots  of  interpolating  polynomials 
must  be  distinct  but  they  are  otherwise  arbitrary,  allowing 
considerable  flexibility  in  the  size  of  the  time  step.  It  has 
changed  six  orders  of  magnitude  during  our  calculations  of 
equilibrium  states. 

An  initial  estimate  for  the  solution  at  the  nth  time  step, 
yn0  is  determined  from  prior  solutions  y ,  (i  =  1,. and 
is  subsequently  corrected  using  the  residual  function 
yno~F(y/Io,rn).  The  order  of  interpolation  q  and  the  stepsize 
h  are  varied  in  response  to  estimates  of  local  discretization 
error  and  user  provided  tolerances.  The  order  of  the  method 
q  is  restricted  to  change  in  steps  of  ±  1  only. 
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The  transmission  coefficient  of  a  layered  structure  made  of  glass  and  water  was  calculated  using 
transfer  matrix  method  and  also  measured  as  a  function  of  frequency.  It  was  found  that  acoustic 
band  gaps  can  be  created  using  only  3-4  cells  of  a  two-phase  layered  structure.  By  introducing  two 
or  more  periods  into  the  layered  structure,  very  sharp  passbands  and  very  broad  stopbands  can  be 
engineered  for  acoustic  waves.  Such  acoustic  band-gap  materials  could  be  used  for  making 
high-quality  acoustic  filters,  acoustic  mirrors  and  vibration  insulation  devices  in  selective  frequency 
range.  ©  1999  American  Institute  of  Physics.  [S0003-695 1(99)0 1549-1] 


Photonic  band-gap  engineering  in  mesoscopic  systems 
has  been  successfully  realized  by  using  crystals  with  periodic 
structures.1,2  The  fundamental  principle  is  to  use  the  periodic 
medium  to  regulate  the  propagation  of  electromagnetic 
waves.  In  the  frequency  window  of  the  band  gaps,  electro¬ 
magnetic  waves  cannot  go  through  the  crystal.  In  principle, 
the  band-gap  phenomena  can  be  produced  when  the  wave¬ 
length  of  the  propagating  waves  approaches  the  period  of  the 
medium.  As  an  analogue  to  photonic  band-gap  engineering, 
such  band-gap  phenomena  also  occur  in  the  case  of  acoustic 
waves  propagating  in  a  periodic  medium  with  the  period 
comparable  to  the  wavelength.  Passbands  and  stopbands  are 
produced  in  the  frequency  spectrum.  In  the  frequency  range 
of  an  acoustic  (or  phonotic)  stopband,  sound  waves  or  ultra¬ 
sonic  vibrations  cannot  propagate.  Such  acoustic  band-gap 
materials  can  have  many  practical  applications,  such  as 
elastic/acoustic  wave  filters,  ultrasonic  silent  blocks,  acoustic 
mirrors,  and  ultrasonic  array  transducers.  In  the  past,  some 
interesting  microscopic  physical  phenomena  have  been  du¬ 
plicated  using  ultrasonic  waves  with  scale-up  designs,  such 
as  the  Anderson  localization  of  sound  and  vibrations;3”5  it  is 
therefore  natural  to  extend  the  photonic  band-gap  concept  to 
acoustic  waves.  This  letter  reports  an  interesting  study  of 
using  multiple  periodicity  and  acoustic  impedance  difference 
to  engineer  large  band  gaps  and/or  narrow  passbands  using 
layered  structure  of  finite  size. 

The  materials  under  study  are  comprised  of  repeated 
cells  made  of  two  different  materials  with  large  acoustic  im¬ 
pedance  difference.  The  length  scale  of  the  period  is  the 
same  as  the  wavelength  of  the  acoustic  waves.  There  has 
been  growing  interest  in  the  band-gap  materials  in  recent 
years.  Several  theoretical  methods  were  developed  in  the 
past  few  years  to  study  the  band-gap  formation,  for  example, 
the  transfer  matrix  method,6  plane  wave  method,7,8  and  ef¬ 
fective  medium  method.9  The  first  experimental  study  on  ul¬ 
trasonic  band  gaps  reported  in  the  literature  was  the  measure¬ 
ment  of  sound  attenuation  in  a  sculpture.10  Ultrasonic  band 
gaps  were  also  experimentally  observed  recently  in  one¬ 
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dimensional  arrays  and  two-dimensional  periodic  composites 
for  longitudinal  waves.11,12 

In  this  letter,  we  report  a  combined  theoretical  and  ex¬ 
perimental  study  on  a  finite  one-dimensional  layered  medium 
containing  two  different  periods.  The  main  objective  is  to 
use  the  idea  of  band-gap  engineering  to  fabricate  materials 
that  can  have  narrow  passbands  and/or  very  broad  stopbands. 

As  shown  in  Fig.  1,  the  structure  is  made  of  two  sub¬ 
structures  with  period  dx  and  d2,  respectively,  and  each  of 
them  consists  of  four  cells.  Each  substructure  can  produce  its 
own  band  structures.  If  the  passbands  of  the  two  substruc¬ 
tures  are  not  overlapped  in  certain  frequency  range,  acoustic 
waves  cannot  pass  through  the  entire  structure  creating  a 
broader  stopband;  if  they  are  partially  overlapped,  only  the 
overlapped  frequencies  can  pass  through,  thus,  a  filter  of  nar¬ 
row  passband  could  be  obtained. 

In  Fig.  1,  element  type  1  and  type  3  are  glass  plates 
while  type  2  and  type  4  are  both  water.  A  longitudinal  acous¬ 
tic  plane  wave  was  sent  from  the  left  at  x=0.  At  time  L  the 
wave  function  in  the  nth  layer  can  be  written  as 

(j  =  1  -4,  n  =  l-N),  (1) 


H-di-H  h—d2'*{ 


FIG.  1.  Schematic  of  the  layered  structure  studied,  which  consists  of  two 
periodic  substructures  with  period  d\  and  d2>  respectively.  Materials  1  and 
3  are  two  different  glasses  and  materials  2  and  4  are  both  water, 
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FIG.  2.  Calculated  and  measured  transmission  coefficient  as  a  function  of 
frequency  for  a  layered  medium  containing  glass  (1)  and  water  (2).  The 
thickness  of  glass  is  1.23  mm  and  the  water  gap  is  1.27  mm. 

where  the  first  and  second  terms  on  the  right-hand  side  of 
Eq.  (1)  represent  the  forward  and  reflected  waves,  respec¬ 
tively,  kj=277f/cj  is  the  wave  vector,  and /and  Cj  are  the 
frequency  and  phase  velocity,  respectively.  The  indices  n 
and  j  are  the  cell  number  and  material  type,  respectively. 

The  transfer  matrix  technique  of  Ref.  6  was  used  to  cal¬ 
culate  the  acoustic  band  structures  for  the  design  given  in 
Fig.  1.  The  properties  of  the  materials  used  in  the  calcula¬ 
tions  are  as  follows: 

Material  1:  glass  with  density  2459  kg/m3  and  phase  velocity 
5660  m/s. 

Material  3:  glass  with  density  2767  kg/m3  and  phase  velocity 
5784  m/s. 

Materials  2  and  4:  distilled  water  with  density  1000  kg/m3 
and  phase  velocity  1480  m/s. 

Dissipation  was  included  in  the  numerical  calculations 
by  adding  a  small  imaginary  component  to  the  phase  velocity 
in  the  glasses,  i.e.,  c1  =  (5660+ 10/)  m/s  and  c3  =  (5784 
+ 100  m/s. 

The  experiments  were  conducted  in  a  water  tank  using  a 
setup  similar  to  the  one  described  in  Refs.  12  and  13.  Two 
broadband  ultrasonic  transducers  were  used  for  the  measure¬ 
ments,  one  as  transmitter  and  the  other  as  receiver.  In  order 
to  cover  the  interested  frequency  range  of  0.5-2.6  MHz,  two 
pairs  of  broadband  ultrasonic  transducers  were  used.  The 
center  frequencies  of  the  transducers  are  2.5  and  1.5  MHz, 
respectively,  and  the  nominal  active  diameter  of  the  trans¬ 
ducers  is  12  mm.  The  transmitting  transducer  was  driven  by 
a  DPR35  pulser/receiver  and  the  transmitted  signal  was  re¬ 
ceived  by  a  digital  oscilloscope  (Tektronix  TDS  460A  with 
fast  Fourier  transform  analysis  capabilities),  then  down¬ 
loaded  to  a  personal  computer.  Ten-signal  average  scheme 
was  used  to  improve  the  signal-to-noise  ratio  during  the  mea¬ 
surements. 

Shown  in  Fig.  2  are  the  calculated  and  measured  fre¬ 
quency  spectra  of  the  transmission  coefficient  for  a  periodic 
layered  structure  with  dimensions  given  in  the  figure  caption. 
There  are  four  passbands  between  the  frequencies  of  0.5  and 
2.6  MHz.  Gibbs  type  oscillations  can  be  clearly  seen  in  the 
calculated  passbands  and  they  are  also  confirmed  by  experi¬ 
ments  within  the  experimental  error.  The  Gibbs  type  oscilla¬ 
tions  here  is  similar  to  that  reported  in  Ref.  12. 
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Frequency  (MHz) 

FIG.  3.  Calculated  and  measured  transmission  coefficient  for  the  structure 
shown  in  Fig.  1  as  a  function  of  frequency.  The  thickness  of  materials  1,  2, 
3,  and  4  are  1.23,  1.14,  3.00,  and  0.51  mm,  respectively. 

Figure  3  shows  the  calculated  and  measured  frequency 
spectra  of  the  transmission  coefficient  for  the  structure 
shown  in  Fig.  1,  which  has  two  different  periods.  It  is  in¬ 
triguing  to  see  that  a  very  narrow  passband  is  indeed  ob¬ 
tained  near  1.9  MHz.  The  agreement  between  the  theoretical 
and  experimental  results  is  very  good.  Figure  4  shows  the 
frequency  spectra  of  the  transmission  coefficient  of  another 
design,  for  which  the  two  substructures  do  not  have  over¬ 
lapped  passbands.  As  a  result,  no  significant  signal  can  pass 
through  in  the  whole  measured  frequency  range.  Again,  the 
calculated  and  experimental  measured  results  agreed  very 
well. 

In  conclusion,  we  have  successfully  engineered  acoustic 
band-gap  materials  with  either  a  very  broad  stopband  or  a 
very  sharp  passband.  By  changing  the  dimensions  and  type 
of  materials  used  in  the  structure,  the  acoustic  band-gap  en¬ 
gineering  concept  of  using  multiple  periodicity  can  be  gen¬ 
eralized  to  higher  or  lower  frequency  ranges.  The  produced 
acoustic  band-gap  materials  have  many  potential  applica¬ 
tions,  such  as  acoustic  filters  allowing  only  selective  frequen¬ 
cies  to  pass  through,  ultrasonic  silent  block  which  can  pro¬ 
vide  a  vibration-free  environment  for  high  precision  devices, 
and  acoustic  mirrors. 


Frequency  (MHz) 

FIG.  4.  Calculated  and  measured  transmission  coefficient  as  a  function  of 
frequency  for  another  design  of  the  structure  shown  in  Fig.  1.  The  thickness 
of  materials  1,  2,  3,  and  4  are  1.23,  1.27,  3.00,  and  0.51  mm,  respectively. 
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Abstract  --  Electric  field  induced  acoustic 
emission  (AE)  in  ferroelectric  PZT  ceramics  has 
been  studied.  It  was  shown  that  the  field  induced 
AE  originated  from  where  the  displacement  was 
induced  and  ended  when  the  applied  field  was 
being  reduced  during  poling.  A  larger  electric 
field  than  the  previously  applied  maximum  field 
was  required  for  new  AE  generation  (Kaiser 
effect).  From  the  AE  event  rate  and  the  AE 
signal  amplitude  distribution,  it  was  assumed 
that  there  were  two  origins  of  the  bipolar  field 
induced  AE  in  the  ferroelectric  ceramics: 
deformation  related  to  domain  reorientation  and 
piezoelectric  deformation  without  domain 
reorientation. 

§1.  Introduction 

In  ferroelectric  materials,  the  AE  method  has 
been  used  to  determine  phase  change, to 
detect  domain  reorientation^)  and  to  monitor 
crack  propagation.^)  It  has  been  believed  that 
ferroelectric  domain  motion  is  a  dominant  AE 
source  in  the  ferroelectrics.  However,  our  recent 
research  has  shown  a  significant  measurement 
problem  in  the  most  of  previous  field  induced 

AE  studies.)®)  A  mechanical  vibration  of  sample 
is  easily  excited  by  an  electrical  coupling 
between  the  power  supply  feedback  and  the 
piezoelectricity  of  sample,  resulting  in  a  vibro- 
acoustic  emission.  The  vibro-acoustic  emission 
could  vary  with  electric  field,  and  was  similar  to 
the  reported  field  induced  AE  results.  The 
sample  vibration  could  be  eliminated  by 
increasing  the  time  constant  T  of  the  voltage 
(field)  application  system.  It  was  shown  that 
this  modified  voltage  (field)  application  method 
makes  the  field  induced  AE  measurements  at 
high  sensitivities  possible  without  vibro-acoustic 
emissions.  Considerably  different  AE  results 
from  the  reported  field  induced  AE  were  obtained 
using  the  modified  method.  The  AE  appeared  to 
be  observed  from  where  the  displacement  was 
induced  during  poling.  New  AE  will  not  be 
generated  after  sufficient  poling. 

This  paper  deals  with  the  electric  field  induced 
AE  generations  in  a  ferroelectric  lead  zirconate 


titanate  Pb(ZrJ,Ti1.I)03  (PZT)  ceramics.  The 
field  induced  AE  activities  were  analyzed  by 
comparing  them  with  obtained  induced 
displacement  A  fractal  dimension  analysis  of 
the  AE  signal  amplitude  distribution  was  also 
applied. 

§2.  Experimental  Setup 

Figure  1  shows  the  field  induced  AE  and 
displacement  measurement  system.  The  sample 
was  placed  on  the  AE  sensor  (NF  Corporation, 
AE-904E)  with  its  resonant  frequency  of  450 
kHz.  The  induced  displacement  was  also 
observed  using  a  L.V.D.T  (Millitoron,  Nr. 
1301).  The  AE  signal  was  amplified  by  40  cB 
through  a  low  noise  pre-amplifier  (NF 
Corporation,  AE-9913)  and  again  up  to  60  cB 
with  a  main  amplifier  (NF  Corporation,  AE- 
922).  The  amplified  AE  signals  were  counted 
after  passing  through  a  high-pass  filter  (f  >100 
kHz)  and  a  discriminator. 


Discriminator 


Figure  1  Field  induced  AE  and  displacement 
measurement  system. 

Disc  specimens  of  commercialized  ferro¬ 
electric  lead  zirconate  titanate  PZT  ceramic 
(Morgan  Matroc,  PZT-5A)  with  12.7  mm  in 
diameter  and  0.4  mm  in  thickness  w ere 
employed  for  this  study.  Gold  electrodes  were 
formed  on  both  surfaces  by  sputtering.  The  AE 
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sensor  was  attached  to  the  sample  through 
brass-foil  with  0.025  mm  in  thickness.  A 
silicone  grease  was  used  as  a  couplant  to  make  a 
better  acoustical  contact.  An  external  resistor  R 
(R  =  20  MQ)  and  a  capacitor  C  (C=  0.01  |aF) 
were  installed  to  the  power  supply  in  order  to 
eliminate  the  vibro-acoustic  emission. 

§3.  Results  and  Discussions 

3.1  The  Kaiser  Effect  in  Field  Induced  AE 

Figure  2  shows  the  AE  event,  induced 
displacement  and  applied  electric  field  as  a 
function  of  time  for  a  cyclic  application  of 
various  electric  fields.  The  threshold  level  of  the 
AE  signal  were  set  to  be  400  mV  at  100  dB.  An 
unpoled  sample  was  employed  for  this 
measurement.  At  the  first  cycle  (poling),  an 
electric  field  of  20  kV/cm  was  applied  to  the 
sample.  The  amplitude  of  the  field  was  increased 
to  be  25  kV/cm  at  the  second  cycle,  and  again 
up  to  30  kV/cm  at  the  third  cycle. 


Figure  2  The  AE  event,  induced  displacement 
and  applied  electric  field  as  a  function  of  time 
for  a  cyclic  application  of  various  electric  fields. 


The  sample  was  poled  at  the  first  cycle.  The 
induced  displacement  and  total  residual 
displacement  observed  as  a  zero  point  shift 
increased  with  amplitude  of  applied  field, 
because  the  higher  degree  of  poling  is  achieved 
by  the  applications  of  the  larger  electric  fields. 
At  first  cycle,  the  AE  event  occurred  from  where 
the  displacement  was  induced  and  increased 
with  field.  The  AE  was  not  generated  during  the 
field  was  being  reduced.  In  the  second  cycle,  no 
AE  events  were  observed  until  the  electric  field 
reached  E  =  20  kV/cm.  New  AE  generation  was 
observed  above  E  =  20  kV/cm,  and  increased  with 
the  electric  field.  In  the  third  cycle,  this  AE 
onset  field  again  shifted  to  be  E  =  25  kV/cm. 


The  AE  was  generated  only  when  the  applied 
field  exceeded  the  previous  maximum  value. 
This  pre- applied  electric  field  (induced 
displacement)  dependence  of  the  AE  generation 
implies  a  Kaiser  effect  in  the  electric  field  induced 
deformation  of  the  PZT  ceramics. 

The  domain  motion  has  been  believed  to  be 
origin  of  the  field  induced  AE.  However,  it  is 
shown  that  all  domain  motions  can  not  be  the 
origin  of  the  field  induced  AE,  because  the  field 
induced  AE  depends  on  the  pre-applied  electric 
field  (Kaiser  effect).  Therefore,  particular  domain 
reorientation  processes  rather  than  domain 
motion  should  be  considered  to  be  the  origins 
of  the  AE.  Another  possible  AE  source  is  a 
mechanical  stress,  since  the  mechanical  stress  is 
supposed  to  be  induced  in  the  sample  at  high 
strain  (electric  field)  state. 


3.2  The  Electric  Field  Induced  AE 

Figure  3  shows  the  AE  event  count  rate  and 
the  induced  displacement  as  a  function  of  applied 
field  at  0.0015  Hz  with  field  of  E  =  ±35  kV/cm. 
The  AE  measurement  conditions  were  set  to  be 
at  100  dB  with  a  threshold  level  of  400  mV.  A 
butterfly  shape  induced  displacement  due  to  the 
domain  reorientation  was  obtained.  The  critical 
electric  field,  where  the  AE  started  to  be 
generated,  corresponded  to  the  point  of  inflection 
of  the  displacement  with  respect  to  the  field.  (In 
Fig  3,  d( displacement ydE  is  shown  for  the  point 
of  inflection.) 

The  maximum  of  AE  event  count  rate  was 
not  observed  at  the  maximum  applied  field.  The 
maximum  AE  event  rate  was  found  to  be  around 
E  =  27  kV/cm.  Considering  that  internal  stress 
increases  with  the  applied  field  in  the 
ferroelec tries,  this  decrease  in  the  AE  event  rate 
might  imply  that  the  internal  stress  could  be  the 
origin  of  the  AE,  but  not  only  one  source.  The 
induced  displacement  in  the  ferroelectric  ceramics 
consists  of  the  deformation  related  to  domain 
reorientation  and  piezoelectric  deformation 
without  the  domain  reorientation.  Therefore,  it 
was  supposed  that  the  field  induced  AE  under 
bipolar  electric  fields  in  the  PZT  ceramics  was 
generated  first  through  deformation  related  to 
domain  reorientation.  After  the  domain 
reorientation  was  completed,  the  piezoelectric 
deformation  unrelated  to  domain  reorientation, 
which  accompanied  with  the  induced  internal 
stress,  was  expected  to  be  the  origin  of  the  AE. 
The  decrease  in  the  AE  event  rate  might  indicate 
the  completion  of  the  domain  reorientation 
related  deformation. 
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Figure  3  AE  event  count  and  the  induced 
displacement  as  a  function  of  applied  field. 


3.3  The  Fractal  Dimension  of  the  Electric  Field 
Induced  AE 

When  the  AE  event  rate  /(x)is  the  minus 
m-th  power  of  X ,  where  X  is  the  AE  signal 
amplitude, 


existence  of  ferroelectric  domain  and  the  effect  of 
the  deformation  related  to  domain  reorientation, 
instead  of  the  crack  in  the  material,  since 
repeatable  and  stable  field  induced  AE  events 

were  observed. 


Figure  4  AE  event  count  per  cycle  as  a  function 
of  the  AE  signal  threshold  level  for  various 
bipolar  electric  field  amplitudes. 


f(x)  =  oTm  (c:  constant).  (1) 

The  number  of  m  is  defined  as  a  fractal 
dimension.11)  The  fractal  dimension  m  is  used 
to  estimate  the  damage  in  the  materials. 

The  integrated  AE  event  F(x ), 

F(x)  =  \°°  f(x)  dx  =  -l—cx~m+1  (m>l)  (2) 
Jx  m  —  1 

was  observed  through  changing  the  AE  signal 
threshold  level  in  the  measurements. 

Figure  4  shows  the  AE  event  count  per  cycle 
as  a  function  of  the  AE  signal  threshold  level. 
The  observed  AE  event  count  logarithmically 
decreased  with  the  AE  signal  threshold  level 
when  a  field  of  E  =  ±  25  kV/cm  was  applied. 
Thus,  a  fractal  dimension  of  m=l 


the  logarithmic  decrease.  In  general,  the  fractal 
dimension  is  found  to  be  more  than  2  for  the  AE 
generated  from  the  plastic  deformation.  The 
fractal  dimension  of  m  =  2  is  a  critical  number 
to  evaluate  the  condition  of  materials.  If  there  is 
a  damage  (i.e.,  crack)  in  the  material,  the  fractal 
dimension  becomes  lower  than  2.  From  this 
concept  of  the  fractal  dimension,  the  obtained 
dimension  indicated  that  the  sample  had  been 
damaged.  However,  it  is  supposed  that  the 
lowered  fractal  dimension  was  due  to  the 


When  the  amplitude  of  electric  field  was 
increased,  the  higher  AE  signal  threshold  level 
parts  did  not  follow  to  the  logarithmic  change. 
A  critical  amplitude  of  the  electric  field,  above 
which  non-logarithmic  AE  distributions  were 
observed,  was  found  to  be  around  E  =  27  kV/cm. 
The  lower  signal  threshold  level  region  still 
could  be  fitted  to  a  logarithmic  curve  and  the 
extended  line  also  ended  around  the  AE  signal 
threshold  level  obtained  for  ±  25  kV/cm.  The 
higher  AE  signal  threshold  level  parts  (>10 
mV)  showed  a  linear  decrease  with  the  threshold 
level.  When  an  electric  field  of  E  =  ±35  kV/cm 
was  applied,  a  fractal  dimension  of  m  =  2.8  was 
obtained.  Figure  5  shows  the  induced 
displacement  under  unipolar  and  bipolar  fields. 
Although  the  butterfly  shape  induced 
displacement  for  the  bipolar  field  is  due  to  the 
domain  reorientation,  the  bipolar  field  induced 
displacement  coincided  with  the  unipolar  field 
induced  displacement  above  the  electric  field  E  = 
27  kV/cm.  Thus,  it  is  assumed  that  the  domain 
reorientation  was  completed  at  higher  electric 
fields  than  the  critical  electric  field  E  =  27 
kV/cm.  The  origin  of  the  AE  with  the  fractal 
dimension  of  2.8  was  assumed  not  to  be  the 
deformation  related  to  domain  reorientation,  but 
the  piezoelectric  deformation  without  domain 
reorientation. 
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Applied  electric  field  (kV/cm) 


Figure  5  Field  induced  displacement  under 
unipolar  and  bipolar  fields. 


domain  reorientation  and  piezoelectric 
deformation  without  domain  reorientation.  It  is 
assumed  that  the  fractal  dimension  of  the 
ferroelectric  PZT  ceramics  consisted  of  the 
domain  reorientation  related  m  =  1  and  the 
domain  reorientation  unrelated  m  =  2.8.  A 
critical  bipolar  electric  field  for  the  domain 
reorientation  unrelated  AE  generation  was  found 
to  be  E  =  37  kV/cm. 
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It  was  shown  that  the  fractal  dimension  of  the 
field  induced  AE  for  the  deformation  related  to 
domain  reorientation  (m  =  1)  was  lower  than  that 
for  the  piezoelectric  deformation  without  domain 
reorientation  (m  =  2.8).  The  fractal  dimension 
greater  than  two  corresponded  to  the  plastic 
deformation.  However,  it  was  difficult  to 
categorize  the  electrically  induced  deformation  of 
ferroelectric  ceramics  as  the  plastic  deformation, 
because  the  induced  displacement  recovers  to  the 
initial  state  after  removing  the  field.  Hence,  it 
should  be  noted  that  the  AE  generation  process 
in  the  ferroelectric  ceramics  could  be  different 
from  that  of  reported  AE  in  the  materials  such  as 
metals  and  structural  ceramics. 

§4.  Conclusion 

The  electric  field  induced  acoustic  emissions 
(AE)  in  the  feiToelectric  PZT  ceramics  woe 
studied.  It  is  shown  that  the  Kaiser  effect  in 
terms  of  applied  electric  field  (displacement) 
took  place  in  the  electric  field  induced 
deformation  of  ferroelectric  PZT  ceramics.  The 
AE  event  rate  as  a  function  of  applied  field  and 
AE  signal  threshold  level  distribution  were 
employed  to  identify  the  origin  of  field  induced 
AE.  The  AE  event  rate  in  ferroelectric  PZT 
showed  two  origins:  deformation  related  to 
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ABSTRACT 

Losses  in  piezoelectrics  are  considered  in  general  to  have  three  different  mechanisms;  dielectric, 
mechanical  and  piezoelectric  losses.  This  paper  deals  with  the  phenomenology  of  losses,  first,  then 
how  to  measure  these  losses  separately  in  experiments. 

INTRODUCTION 

Loss  or  hysteresis  in  piezoelectrics  exhibits  both  merits  and  demerits.  For  positioning  actuator 
applications,  hysteresis  in  the  field-induced  strain  provides  a  serious  problem,  and  for  resonance 
actuation  such  as  ultrasonic  motors,  loss  generates  significant  heat  in  the  piezoelectric  materials. 
Further,  in  consideration  of  the  resonant  strain  amplified  in  proportion  to  a  mechanical  quality 
factor,  low  (extrinsic)  mechanical  loss  materials  are  preferred  for  ultrasonic  motors  [1-3].  On  the 
contrary,  for  force  sensors  and  acoustic  transducers,  high  mechanical  loss,  which  corresponds  to  a 
low  mechanical  quality  factor  Qm,  is  essential  to  widen  a  frequency  range  for  receiving  signals. 

However,  not  much  research  effort  has  been  put  into  systematic  studies  of  the  loss  mechanisms 
in  piezoelectrics,  particularly  in  high  voltage  and  high  power  range  [4-6] .  Since  not  many 
comprehensive  descriptions  can  be  found  in  previous  reports  [7],  this  paper  will  clarify  the  loss 
mechanisms  in  piezoelectrics  phenomenologically,  first,  then  how  to  measure  these  losses 
separately  in  experiments. 

GENERAL  CONSIDERATION  OF  LOSS  AND  HYSTERESIS 
Extrinsic  Losses 

We  will  start  from  the  Gibbs  free  energy  G  in  a  piezoelectric  material  expressed  by 

dG  =  -  x  dX  -  D  dE  -  S  dT,  (1) 

and 

G  =  -  (1/2)  sE  X2  -  d  X  E  -  (1/2)  exeo  E2.  (2) 

Here,  x  is  strain,  X,  stress,  D,  electric  displacement,  E,  electric  field,  S,  enthalpy  and  T  is  temperature. 
Equation  (2)  is  the  energy  expression  in  terms  of  intensive  physical  parameters  (externally 
controllable  parameters)  X  and  E.  Temperature  dependence  is  carried  into  the  elastic  compliance 

SE,  the  dielectric  constant  ex  and  the  piezoelectric  constant  d.  We  will  obtain  the  following  two 
piezoelectric  equations: 

x  =  -  0G/9X)  =  sE  X  +  d  E ,  (3) 

D  =  -  0G/9E)  =  d  X  +  exeo  E .  (4) 

Note  that  thermodynamical  equations  and  the  consequent  piezoelectric  equations  (Eqs.  (l)-(4)) 
cannot  yield  a  loss,  without  taking  into  account  irreversible  thermodynamic  equations  or  dissipation 
functions,  in  general.  However,  the  latter  considerations  are  mathematically  equivalent  to  the 
introduction  of  complex  physical  constants  into  the  phenomenological  equations,  if  the  loss  can  be 
treated  as  a  perturbation. 


Therefore,  we  will  introduce  complex  parameters  fr^*,  s^*,  and  d*  in  order  to  consider  the 
hysteresis  losses  in  electric,  elastic  and  piezoelectric  coupling  energy : 


e^*  =  eX(l  -  j  tan  8'), 

(5) 

SE*  =  sE(i  -  j  tan  <t>), 

(6) 

d*  =  d  (1  -  j  tan  0’). 

(7) 

Note  that  the  negative  connection  in  Eqs.  (5)-(7)  come  from  the  time  “delay.”  0’  is  the  phase 
delay  of  the  strain  under  an  applied  electric  field,  or  the  phase  delay  of  the  electric  displacement 
under  an  applied  stress.  Both  delay  phases  should  be  exactly  the  same  if  we  introduce  the  same 

complex  piezoelectric  constant  d*  into  Eqs.(3)  and  (4).  8’  is  the  phase  delay  of  the  electric 
displacement  to  an  applied  electric  field  under  a  constant  stress  (i.e.,  zero  stress)  condition,  and  (|)’  is 
the  phase  delay  of  the  strain  to  an  applied  stress  under  a  constant  electric  field  (i.e.,  short-circuit) 
condition.  We  will  consider  these  phase  delays  as  "extrinsic"  losses. 

Figures  1(a)  -  1(d)  correspond  to  the  model  hysteresis  curves  for  practical  experiments:  D  vs.  E 
curve  under  a  stress-free  condition,  x  vs.  X  under  a  short-circuit  condition,  x  vs.  E  under  a  stress- 
free  condition  and  D  vs.  X  under  an  open-circuit  condition  for  measuring  charge  (or  under  a  short- 
circuit  condition  for  measuring  current),  respectively. 

The  stored  energies  and  hysteresis  losses  for  pure  electrical  and  mechanical  energies  can  be 


calculated  as: 

Ue  =  (1/2)  £Xeo  Eo2,  (8) 

we  =  n  e^eo  Eo^  tan  8’ ;  (9) 

Um  =  (1/2)  sExo2,  (10) 

wm  =  tcsE  Xq2  tan  <)>’ .  (11) 


The  electromechanical  loss,  when  measuring  the  induced  strain  under  an  electric  field,  is  more 
complicated,  because  the  field  vs.  strain  domain  cannot  be  used  for  directly  calculating  the  energy. 
Let  us  calculate  the  stored  energy  Uem  during  a  quarter  electric  field  cycle  (i.e.,  0  to  Eq)  ,  first : 


Fig.  1  (a)  D  vs.  E  (stress  free),  (b)  x  vs.  X  (short-circuit),  (c)  x  vs.  E  (stress  free)  and 
(d)  D  vs.  X  (open-circuit)  curves  with  a  slight  hysteresis  in  each  relation. 


(12) 


Uem  =  -  J  x  dX  =  (1/2)  (xo2^)  =  (1/2)  (dEo)2/sE 
=  (1/2)  (d2/sE)Eo2 

Replacing  d  and  sE  by  d*  =  d  (1  -  j  tan  9’)  and  sE*  =  sE(l  -  j  tan  <J>’),  we  obtain 

Uem=(l/2)(d2/sE)Eo2,  (13> 

and 

wem  =  Jt  (d2/  sE)  Eo2  (2  tan0’  -  tan<t>’).  ( 14) 

Note  that  the  strain  vs.  electric  field  measurement  seems  to  provide  the  piezoelectric  loss  tan0 
directly,  however,  the  observed  loss  should  include  an  additional  elastic  loss  because  the  strain 
should  be  delayed  to  the  initial  stress,  which  is  needed  to  calculate  energy. 

Similarly,  when  we  measure  the  induced  charge  under  stress,  the  stored  energy  Ume  and  the 

hysteresis  loss  wme  during  a  quarter  and  a  full  stress  cycle,  respectively,  are  obtained  as 

Ume  =  (1/2)  (d2/  eoex)  Xo2,  (15) 

and 

wme  =  Jt  (d2/  E0£X)  Xo2  (2  tan0’  -  tan8’).  (16) 

Hence,  from  the  measurements  of  D  vs.  E  and  x  vs.  X,  we  obtain  tan5’  and  tan<j> ,  respectively, 
and  either  the  piezoelectric  (D  vs.  X)  or  converse  piezoelectric  measurement  (x  vs.  E)  provides  tan 
0’  through  a  numerical  subtraction. 

Intrinsic  Losses 

So  far,  we  discussed  the  “extrinsic”  electric,  mechanical  and  piezoelectric  losses.  In  order  to 
consider  real  physical  meanings  of  the  losses,  we  will  introduce  the  “intrinsic”  losses. ,  When  we 
start  from  the  energy  expression  in  terms  of  extensive  physical  parameters  (material  s  intrinsic 


parameters)  x  and  D,  that  is, 

dA  =  Xdx  +  EdD-S  dT,  (17) 

we  can  obtain  the  piezoelectric  equations  as  follows: 

X  =  (9A/9x)  =  cDx-hD,  (18) 

E  =  0A/9D)  =  -  h  x  +  kxko  D  .  (19) 

We  introduce  the  intrinsic  electric,  mechanical  and  piezoelectric  losses  as 

Kx*  =  kx  (1  +j  tan  8),  (20) 

cD*  =  cD(l  +j  tan<f>),  (21) 

h*  =  h  (1  +  j  tan  0).  (22) 


It  is  notable  that  the  permittivity  under  a  constant  strain  (i.  e.,  zero  strain  or  completely  clamped) 
condition,  ex*  and  the  elastic  compiance  under  a  constant  electric  displacement  (i.  e.,  open-circuit) 
condition,  s^*  can  be  provided  as  an  inverse  value  of  kx*  and  c^*,  respectively.  Thus,  using  the 
exactly  the  same  losses  in  Eqs.  (20)  and  (21), 

ex*  =  ex(l  -  j  tan  8), 

SD*  =  sD(i  -  j  tan  <j>). 

We  will  consider  these  phase  delays  again  as  "intrinsic"  losses. 


(23) 

(24) 


Here,  we  consider  the  physical  property  difference  between  the  boundary  conditions;  E  constant 
and  D  constant,  or  X  constant  and  x  constant.  When  an  elecric  field  is  applied  on  a  piezoelectric 
sample  as  illustrated  in  the  top  of  Fig.  2,  this  state  will  be  equivalent  to  the  superposition  of  the 
following  two  steps:  first,  the  sample  is  completely  clamped  and  the  field  Eo  is  applied  (pure 

electrical  energy  (1/2)  ex£0  Eo2  is  input);  second,  keeping  the  field  at  Eo,  the  mechanical  constraint 

is  released  (additional  mechanical  energy  (1/2)  (d2/s^)  Eo2  is  necessary).  The  total  energy  should 

correspond  to  the  total  input  electrical  energy  (1/2)  £^£Q  Eo2;  thus,  we  obtain  the  relation, 


eQgX  =  £QeX  +  (<j2/sE)5  (25) 

Similarly,  from  the  bottom  of  Fig.  2, 

SE  =  SD  +  (d^/eQgX)  (26) 

Hence,  we  obtain  the  following  equations: 

ex/eX  =  (i_k2),  (27) 

sD/sE  =  (l -k2),  (28) 

where 

k2  =  d2/(sE  eoeX).  (29) 


Total  Input  Stored  Stored 

Electrical  Energy  Electrical  Energy  Mechanical  Energy 

=  (1/2)  eoeX  Eo2  =  (1/2)  eoex  Eg2  =  (1/2)  (d2/sE)Eo2 


clamped  Electrically 

short-circuited 


Total  Input 
Mechanical  Energy 
=  (1/2)sEXq2 


Electrically 

short-circuited 


Stored 

Mechanical  Energy 
=  (1/2)  sD  Xg2 


♦ 


I 


Electrically 

open-circuited 


Stored 

Electrical  Energy 
=  (1/2)  (d2/e0ex)Xo2 


Polarization 
P  =  dXo 
Inverse  field 
E  =  dXo/EOE* 


Fig.  2  Conceptual  figure  for  explaining  the  relation  between  eX  and  £x,  s^  and  sP. 


Similarly, 

kx/kx  =  (1  -k2), 
cE/cD=  (1  -  k2), 

where 

k2  =  h2/(cD  kxkq). 


(30) 

(31) 

(32) 


This  k  is  called  the  electromechanical  coupling  factor,  and  is  the  same  as  the  k  in  Eq.  (29),  because 
the  equation  d2/(sE  eoex)  =  h2/(c^  kxko)  can  be  verified  mathematically.  We  define  the  k  as  a 
real  number  in  this  manuscript 

In  order  to  obtain  the  relationships  between  the  extrinsic  and  intrinsic  losses,  the  following  three 
equations  are  essential: 


£()£X  =  [kxko  (1  -  h2/(cD  kxko))]"1, 

sE  =  [CD  (i .  h2/(cD  Kxko))]-1, 

d  =  [h2/(cD  kxko)]  [h  (1  -  h2/(cD  kxkq))]  -1. 


(33) 

(34) 

(35) 


Replacing  the  parameters  in  Eqs.  (33)  -  (35)  by  the  complex  parameters  in  Eqs.  (5)  -  (7),  (20)  - 
(22),  we  obtain  the  relationships  between  the  extrinsic  and  intrinsic  losses: 


tan  5'=  (1/(1  -  k2))[tan  8  +  k2(tan  <(>  -  2  tan  0)], 
tan  (j)'=  (1/(1  -  k2))[tan  <()  +  k2(tan  8  -  2  tan  0)], 
tan  0'=  (1/(1  -  k2))[tan  8  +  tan  <|>  +  (1  +k2)  tan  0], 


(36) 

(37) 

(38) 


where  k  is  the  electromechanical  coupling  factor  defined  by  either  Eq.  (29)  or  Eq.  (32),  and  here  as 
a  real  number.  It  is  important  that  the  extrinsic  dielectric  and  elastic  losses  are  mutually  correlated 
with  the  intrinsic  dielectric,  elastic  and  piezoelectric  losses  through  the  electromechanical  coupling 

k2,  and  that  the  denominator  (1  -  k2)  comes  basically  from  the  ratios,  ex  /ex  =  (1  -  k2)  and  s^/sE  = 

(1  -  k2),  and  this  real  part  reflects  to  the  dissipation  factor  when  the  imaginary  part  is  devided  by  the 
real  part 

EXPERIMENTAL  RESULTS 

Figure  3  shows  “extrinsic”  dissipation  factors  determined  from  (a)  D  vs.  E  (stress  fiee),  (b)  x 
vs.  X  (short-circuit),  (c)  x  vs.  E  (stress  fiee)  and  (d)  D  vs.  X  (open-circuit)  curves  for  a  PZT  based 
soft  multilayer  actuator.  Figure  4  shows  the  result  for  the  piezoelectric  loss  tan0\  We  used  the 

correlation  factor  between  electric  field  and  compressive  stress  given  averagely  by  X  =  (eoex/sE)l/2 
E. 

From  Figs.  3  and  4,  we  can  calculate  the  “intrinsic”  losses  as  shown  in  Fig.  5.  Note  that  the 
piezoelectric  losses  tan0’  and  tan0  are  not  so  small  as  previously  believed,  but  comparable  to  the 
dielectric  and  elastic  losses,  and  increase  gradually  with  the  field  or  stress.  The  experimental  details 
will  be  reported  in  the  successive  papers. 

When  similar  measurements  to  Figs.  1(a)  and  l.(b),  but  under  constrained  conditions;  that  is,  D 
vs.  E  under  a  completely  clamped  state,  and  x  vs.  X  under  an  open-circuit  state,  respectively,  we  can 
expect  smaller  hystereses;  that  is,  intrinsic  losses,  tan8  and  tan<|).  These  measurements  seem  to  be 
alternative  methods  to  determine  the  three  losses  separately,  however,  they  are  rather  difficult  in 
practice. 


(a) 


(b) 


(c) 


(d) 


Fig.3  Dissipation  factors  determined  from  (a)  D  vs.  E  (stress  free),  (b)  x  vs.  X  (short-circuit),  (c)  x 
vs.  E  (stress  free)  and  (d)  D  vs.  X  (open-circuit)  curves  for  a  PZT  based  actuator. 
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Fig.4  Extrinsic  piezoelectric  dissipation  factor  tan0’  as  a  function  of  electric  field  or  compressive 

stress,  measured  for  a  PZT  based  actuator. 


0.15 


0.15 


to 

s 

-M 


Electric  Field  (kV/mm) 


Compressive  Stress  (MPa) 


Electric  Field  (kV/mm) 


Fig.5  Intrinsic  loss  factors,  tan8,  tan<j),  and  tan0  as  a  function  of  electric  field  or  compressive  stress, 

measured  for  a  PZT  based  actuator. 


CONCLUSIONS 

Theoretical  and  experimental  techniques  for  measuring  the  extrinsic  and  intrinsic  electric, 
mechanical  and  piezoelectric  coupling  losses  separately  have  been  discussed.  We  found  that  the 
piezoelectric  loss  is  not  so  small  as  believed  by  the  previous  researchers. 

This  research  was  sponsored  by  the  Office  of  Naval  Research  through  the  MURI  program, 
Grant  No.  N00014-96-1-1 173. 
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LOSS  MECHANISMS  IN  PIEZOELECTRICS 

-  How  to  Measure  Different  Losses  Separately  - 
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ABSTRACT 

Losses  in  piezoelectrics  are  considered  in  general  to  have  three  different 
mechanisms;  dielectric,  mechanical  and  piezoelectric  losses.  This  paper  deals  with  the 
phenomenology  of  losses,  first,  then  how  to  measure  these  losses  separately  in  experiments. 
We  found  that  heat  generation  at  off-resonance  is  caused  mainly  by  dielectric  loss  tan  8’  (i. 
e.  P-E  hysteresis  loss),  not  by  mechanical  loss,  and  that  a  significant  decrease  in  mechanical 
Qm  with  an  increase  of  vibration  level  was  observed  in  resonant  piezoelectric  ceramic 
devices,  which  is  due  to  an  increase  in  the  intrinsic  dielectric  loss,  not  in  the  intrinsic 
mechanical  loss.  Finally,  we  propose  the  usage  of  the  antiresonance  mode  rather  than  the 
conventional  resonance  mode,  particularly  for  high  power  applications,  since  the  mechanical 
quality  factor  Q0  at  an  antiresonance  frequency  is  larger  than  QA  at  a  resonance  frequency. 

*,  **)  Kenji  Uchino  is  with  the  International  Center  for  Actuators  and  Transducers, 
Materials  Research  Laboratory,  The  Pennsylvania  State  University,  University  Park,  Pa 
16802,  and  Seiji  Hirose  is  with  Faculty  of  Engineering,  Yamagata  University,  Yonezawa 
992,  Japan. 
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INTRODUCTION 


Loss  or  hysteresis  in  piezoelectrics  exhibits  both  merits  and  demerits.  For 
positioning  actuator  applications,  hysteresis  in  the  field-induced  strain  provides  a  serious 
problem,  and  for  resonance  actuation  such  as  ultrasonic  motors,  loss  generates  significant 
heat  in  the  piezoelectric  materials.  Further,  in  consideration  of  the  resonant  strain  amplified 
in  proportion  to  a  mechanical  quality  factor,  low  (extrinsic)  mechanical  loss  materials  are 
preferred  for  ultrasonic  motors.  On  the  contrary,  for  force  sensors  and  acoustic 
transducers,  a  low  mechanical  quality  factor  Qm  (which  corresponds  to  high  mechanical 
loss)  is  essential  to  widen  a  frequency  range  for  receiving  signals. 

K.  H.  Haerdtl  wrote  a  review  article  on  electrical  and  mechanical  losses  in 
ferroelectric  ceramics.0  Losses  are  considered  to  consist  of  four  portions:  (1)  domain  wall 
motion,  (2)  fundamental  lattice  portion,  which  should  also  occur  in  domain-free 
monocrystals,  (3)  microstructure  portion,  which  occurs  typically  in  polycrystalline  samples, 
and  (4)  conductivity  portion  in  highly-ohmic  samples.  However,  in  the  typical  piezoelectric 
ceramic  case,  the  loss  due  to  the  domain  wall  motion  exceeds  the  other  three  contributions 
significantly.  They  reported  interesting  experimental  results  on  the  relationship  between 
electrical  and  mechanical  losses  in  piezoceramics,  Pb0  9La0  |(Zr05TiOJ)1  xMex03,  where  Me 

represents  the  doped  ions  Mn,  Fe  or  A1  and  x  varied  between  0  and  0.09.  However,  they 
measured  the  mechanical  losses  on  poled  ceramic  samples,  while  the  electrical  losses  on 
unpoled  samples,  i.e.,  in  a  different  polarization  state.  Thus,  they  completely  neglected 
piezoelectric  losses. 

As  far  as  the  authors’  knowledge,  not  much  research  effort  has  been  put  into 
systematic  studies  of  the  loss  mechanisms  in  piezoelectrics,  particularly  in  high  voltage  and 
high  power  ranges.  Since  not  many  comprehensive  descriptions  can  be  found  in  previous 
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reports,  this  paper  will  clarify  the  loss  mechanisms  in  piezoelectrics  phenomenologically, 
describe  heat  generation  processes  and  high  power  characteristics,  and  finally,  discuss  the 
resonance  and  antiresonance  vibration  modes  from  a  viewpoint  of  a  quality  factor. 

Although  T.  Ikeda  described  part  of  the  formulas  of  this  paper  in  his  textbook,2’  he 
totally  neglected  the  piezoelectric  losses,  which  have  been  found  not  to  be  neglected  in  our 
investigtions.  We  derive  the  full  descriptions  of  the  all  losses  in  this  paper. 

LOSS  AND  HYSTERESIS  IN  THE  POLARIZATION  CURVE 


Relation  between  Hysteresis  and  Dissipation  Factor 

Let  us  start  first  with  loss  and  hysteresis  in  the  electric  displacement  D  (nearly  equal 
to  polarization  P)  vs.  electric  field  E  curve  without  considering  the  electromechanical 
coupling.  Figure  1(a)  shows  an  example  P-E  hysteresis  curve.  When  the  D  (or  P)  traces  a 
different  line  with  increased  and  decreased  applied  electric  field  E,  it  is  called  hysteresis. 


When  the  hysteresis  is  not  very  large,  the  electric  displacement  D  can  be  expressed 
by  using  a  slight  phase  lag  to  the  applied  electric  field.  Assuming  that  the  electric  field 
oscillates  at  a  frequency  f  (=  (0/2tt)  as 


E*  =  £0 


jo* 
e  , 


(1) 


the  induced  electric  displacement  oscillates  also  at  the  same  frequency  under  the  steady 
state,  but  with  some  time  phase  delay  8’: 


£>*  =  £> 

0 


j(ox  -  8’) 

e 


(2) 


If  we  express  the  relation  between  D*  and  E*  as 


D*  =  £*£0  E*, 


(3) 
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(4) 


where  the  complex  dielectric  constant  e*  is 

e  =e  -j  e  , 
and  where 

e”/  e’  =  tan  S’  .  (5) 

Note  that  the  negative  connection  in  Eq.(4)  comes  from  the  time  “delay”,  and  that  £  £0  = 
(Dq/Eo)  cos  S’  and  e”eo  =  {Dq/Eq)  sin  8’ . 

The  area  we  corresponds  to  the  consumed  loss  energy  during  an  electric  field  cycle 

per  unit  volume  of  the  dielectrics,  and  can  be  related  in  isotropic  dielectrics  with  £  or  tan  8’ 
as  follows: 

2ji/g) 

we  =  -JZ)d£  =  -.f  D  (d£/dt)  dt  =  n  Eq  Dq  sin  8’ 
o 

=  n  e’e  £ 2  =  k  e’e  £ 2  tan  8’  (6) 

0  0  0  0 

When  there  is  no  phase  delay  (8’  =  0),  we  =  0;  i.  e.  the  electrostatic  energy  stored  in 
the  dielectric  will  be  recovered  completely  after  a  full  cycle  (100%  efficiency).  However, 
when  there  is  a  phase  delay,  the  loss  we  will  be  accompanied  per  cycle,  and  the  dielectric 
material  generates  heat.  The  tan  8’  is  called  dielectric  dissipation  factor. 

In  consideration  of  the  stored  electrostatic  energy  during  a  half  cycle  from  -  £o  to 
£0  [=  4  Ue,  which  is  illustrated  as  an  area  in  Fig.  1(a)]  provided  by 

4  Ue  =  (l/2)(2  £0)(2  e’e#  £q)  =  2  e’eo  E0\  (7) 

the  dissipation  factor  tan  8’  can  be  experimentally  obtained  by 
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tan  8’  =  ( l/2:t)  (we  /  Ue). 


(8) 


Note  that  we  is  the  hysteresis  in  a  full  cycle  and  Ue  is  the  stored  energy  in  a  quarter  of 
cycle. 

Temperature,  Electric  Field  and  Frequency  Dependence  of  P-E  Hysteresis 

Figures  2, 3  and  4  show  temperature,  electric  field  and  frequency  dependence  of  the 
dissipation  factor  tan  8’  calculated  from  the  P-E  hysteresis  loss  measured  under  stress  free 
condition  for  a  PZT-based  ceramic.  Refer  to  our  paper3)  for  the  experimental  details.  The 
loss  tan  8’  decreases  gradually  with  increasing  temperature,  but  is  rather  insensitive  to 
frequency.  On  the  contrary,  the  tan  8’  increases  initially  in  proportion  to  the  applied  electric 
field,  exhibiting  a  saturation  above  a  certain  electric  field.  After  reaching  the  saturation,  this 
complex  physical  quantity  treatment  should  not  be  employed.  The  value  for  E  =  0  (solid 
triangle  mark  in  the  figure)  was  obtained  with  an  impedance  analyzer. 

GENERAL  CONSIDERATION  OF  LOSS  AND  HYSTERESIS 

Theoretical  Formulas 

Let  us  expand  the  above  discussion  into  more  general  cases;  i.e.,  piezoelectric 
materials.  We  will  start  from  the  Gibbs  free  energy  G  expressed  by 

dG  =  -xdX-DdE-SdT,  (9) 

or 

G  =  -  (1/2)  sE  X2  -  d  X  £  -  (1/2)  £X£0  E2.  (10) 

Here,  x  is  strain,  X,  stress,  D,  electric  displacement,  E,  electric  field,  S,  enthalpy  and  T  is 
temperature.  Equation  (10)  is  the  energy  expression  in  terms  of  intensive  (i.e.,  externally 
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controllable)  physical  parameters  X  and  E.  Temperature  dependence  is  carried  into  the 

E  x 

elastic  compliance  s  ,  the  dielectric  constant  £  and  the  piezoelectric  constant  d.  We  will 
obtain  the  following  two  piezoelectric  equations: 


*=-(aG/ax)=sEx+d£, 

(11) 

D  =  -(dG/dE)  =  dX  +  ZE0E. 

(12) 

Note  that  thermodynamical  equations  and  the  consequent  piezoelectric  equations 
[Eqs.  (9)-(12)]  cannot  yield  a  delay-time  related  loss,  without  taking  into  account 
irreversible  thermodynamic  equations  or  dissipation  functions,  in  general.  However,  the 
latter  considerations  are  mathematically  equivalent  to  the  introduction  of  complex  physical 
constants  into  the  phenomenological  equations,  if  the  loss  can  be  treated  as  a  perturbation. 

Therefore,  we  will  introduce  complex  parameters  e  ,  s  and  d  in  order  to 
consider  the  hysteresis  losses  in  dielectric,  elastic  and  piezoelectric  coupling  energy : 


eX*  =  eX(l  -  j  tan  5’), 

(13) 

SE*  =  sE(l  -  j  tan  <{>’), 

(14) 

d*  =  d(l  -  j  tan  0’). 

(15) 

0’  is  the  phase  delay  of  the  strain  under  an  applied  electric  field,  or  the  phase  delay  of  the 
electric  displacement  under  an  applied  stress.  Both  delay  phases  should  be  exactly  the  same 
if  we  introduce  the  same  complex  piezoelectric  constant  d*  into  Eqs.(l  1)  and  (12).  8’  is  the 
phase  delay  of  the  electric  displacement  to  an  applied  electric  field  under  a  constant  stress 
(e.g.,  zero  stress)  condition,  and  (I)’  is  the  phase  delay  of  the  strain  to  an  applied  stress  under 
a  constant  electric  field  (e.g.,  short-circuit)  condition.  We  will  consider  these  phase  delays 
as  “extrinsic”  losses. 
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Figures  1(a)  -  1(d)  correspond  to  the  model  hysteresis  curves  for  practical 
experiments:  D  vs.  E  curve  under  a  stress-free  condition,  x  vs.  X  under  a  short-circuit 
condition,  x  vs.  E  under  a  stress-free  condition  and  D  vs.  X  under  an  open-circuit  condition 
for  measuring  charge  (or  under  a  short-circuit  condition  for  measuring  current), 
respectively.  Notice  that  these  measurements  are  easily  conducted  in  practice. 

In  a  similar  fashion  to  the  previous  section,  the  stored  energies  and  hysteresis  losses 
for  pure  dielectric  and  elastic  energies  can  be  calculated  as: 


Ue  =  (1/2)  £X£0  E2, 

(16) 

X  2 

we  =  k  £  eo  Eq  tan  8’ , 

(17) 

Um  =  (l/2)sEX02, 

(18) 

E  2 

wm  =  k  s  XQ  tan  (J)’ . 

(19) 

The  electromechanical  loss,  when  measuring  the  induced  strain  under  an  electric 
field,  is  more  complicated.  Let  us  calculate  the  stored  energy  Uem  during  a  quarter  electric 
field  cycle  (i.e.,  0  to  Eo) ,  first : 

Uem  =  -  J x dX  =  (1/2)  (x*/*)  =  (1/2)  (dEQ)2/sE 

=  (1/2)  (d2/sE)Eo2.  (20) 

E  )k  Ejjc  E 

Replacing  d  and  s  by  d  =  d(  1  -  j  tan  0’)  and  s  =  s  (1  -  j  tan  4>’),  we  obtain 
Uem  =  d/2)(d2/sE)E02,  (21) 

and 

wem  =  n  (d2/  sE)  E*  (2  tan0’  -  tan<j>’ ).  (22) 
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Note  that  the  strain  vs.  electric  field  measurement  seems  to  provide  the  piezoelectric 
loss  tan0’  directly,  however,  the  observed  loss  should  include  an  additional  elastic  loss 
because  the  strain  should  be  delayed  to  the  initial  stress,  which  is  needed  to  calculate  energy. 

Similarly,  when  we  measure  the  induced  charge  under  stress,  the  stored  energy  Ume 
and  the  hysteresis  loss  wme  during  a  quarter  and  a  full  stress  cycle,  respectively,  are 
obtained  as 

Ume  =  (l/2)(d2/e0eX)X02,  (23) 

and 

wme  =  « (d2/  eoeX)  XQ2  (2  tan©’  -  tan8’).  (24) 

Hence,  from  the  measurements  of  D  vs.  E  and  *  vs.  X,  we  obtain  tan8’  and  tan<|)’, 

respectively,  and  either  the  piezoelectric  (D  vs.  X)  or  converse  piezoelectric  measurement  (x 

\ 

vs.  E)  provides  tan  0’  through  a  numerical  subtraction. 

So  far,  we  discussed  the  “extrinsic”  dielectric,  mechanical  and  piezoelectric  losses. 
In  order  to  consider  real  physical  meanings  of  the  losses,  we  will  introduce  the  “intrinsic” 
losses.25  When  we  start  from  the  energy  expression  in  terms  of  extensive  (material’s  own) 


physical  parameters  x  and  D,  that  is, 

dA  =  Xdx  + EdD-SdT,  (25) 

we  can  obtain  the  piezoelectric  equations  as  follows: 

X  =  (dAJdx)  =  c°  x-hD ,  (26) 

E  =  (dA/dD)  =  -hx  +  k*kq  D .  (27) 

We  introduce  the  intrinsic  dielectric,  elastic  and  piezoelectric  losses  as 
KX*  =  K(l+jtan8),  (28) 
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(29) 


Dj|c  D  .  ,v 

c  =  c  (l+jtan<}>), 

h*  =  h  (1  +  j  tan  0).  (30) 

It  is  notable  that  the  permittivity  under  a  constant  strain  (e.g.,  zero  strain  or 
completely  clamped)  condition,  £  *  and  the  elastic  compiance  under  a  constant  electric 
displacement  (e.g.,  open-circuit)  condition,  s°*  can  be  provided  as  an  inverse  value  of  k** 
and  c°*,  respectively.  Thus,  using  the  exactly  the  same  losses  in  Eqs.  (28)  and  (29), 

£**  =  £X(1  -  j  tan  5),  (31) 

s°*  =  s°(  1  -  j  tan  <(>),  (32) 

We  will  consider  these  phase  delays  again  as  "intrinsic"  losses. 

Here,  we  consider  the  physical  property  difference  between  the  boundary 
conditions;  E  constant  and  D  constant,  or  X  constant  and  x  constant.  When  an  elecric  field 
is  applied  on  a  piezoelectric  sample  as  illustrated  in  the  top  of  Fig.  5,  this  state  will  be 
equivalent  to  the  superposition  of  the  following  two  steps:  first,  the  sample  is  completely 

clamped  and  the  field  EQ  is  applied  (pure  electrical  energy  (1/2)  £X£0  E*  is  input);  second, 
keeping  the  field  at  Efl,  the  mechanical  constraint  is  released  (additional  mechanical  energy 

(1/2)  (d2/sE)  £q2  is  necessary).  The  total  energy  should  correspond  to  the  total  input 

X  2 

electrical  energy  (1/2)  £  eq  EQ ;  thus,  we  obtain  the  relation, 

£0£X  =  £Q£X  +  (d2/sE),  (33) 

Similarly,  from  the  bottom  of  Fig.  5, 

s  =  s  +  (d  /£0£  ).  (34) 

Hence,  we  obtain  the  following  equations: 

£/EX=(l-k2),  (35) 
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sD/sE  =  (l-k2),  (36) 

where 

k2  =  d2/(sEe/).  (37) 

Similarly, 

kX/k=(1 -k2),  (38) 

cE/c°=  (1-k2),  (39) 

where 

k!  =  hVlc\).  (40) 


This  k  is  called  the  electromechanical  coupling  factor,  and  is  the  same  as  the  k  in  Eq.  (37), 

2  E  X  2  D  x 

because  the  equation  d  /(s  eQe  )  =  h  /(c  k  kq)  can  be  verified  mathematically.  We  define 
the  k  as  a  real  number  in  this  manuscript. 

In  order  to  obtain  the  relationships  between  the  extrinsic  and  intrinsic  losses,  the 


following  three  equations  are  essential: 

£oeX  =  [kxk0  (1  -  h2/(c°  kxk0))]'',  (41) 

sE  =  [cD(l-h2/(cDKK0))]1,  (42) 

d  =  [h2/(c°  kk0)]  [h  (1  -  h2/(c°  kxk0))]''.  (43) 

Replacing  the  parameters  in  Eqs.  (41)  -  (43)  by  the  complex  parameters  in  Eqs.  (13)  -  (15), 
(28)  -  (30),  we  obtain  the  relationships  between  the  extrinsic  and  intrinsic  losses: 
tan  8'=  (1/(1-  k2))[tan  8  +  k2(tan  <|>  -  2  tan  0)],  (44) 

tan  <j>'=  (1/(1-  k2))[tan  <|>  +  k2(tan  8  -  2  tan  0)],  (45) 

tan  0'=  (1/(1  -  k2))[tan  8  +  tan  <|>  +  ( 1  +k2)  tan  0],  (46) 


where  k  is  the  electromechanical  coupling  factor  defined  by  either  Eq.  (37)  or  Eq.  (40),  and 
here  as  a  real  number.  It  is  important  that  the  extrinsic  dielectric  and  elastic  losses  are 
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mutually  correlated  with  the  intrinsic  dielectric,  elastic  and  piezoelectric  losses  through  the 
electromechanical  coupling  k ,  and  that  the  denominator  (1  -  k  )  comes  basically  from  the 
ratios,  £*/£*  =  (1  -  k2)  and  sD/sE  =  (1  -  k2),  and  this  real  part  reflects  to  the  dissipation  factor 

when  the  imaginary  part  is  devided  by  the  real  part.  Also  note  that  depending  on  the 

2 

vibration  mode,  the  definition  of  electromechanical  coupling  k  can  be  changed  such  as  k  = 
d  /(s  eoe  ). 

Experimental  Example 

Figure  6  shows  “extrinsic”  dissipation  factors  determined  from  (a)  D  vs.  E  (stress 
free),  (b)  x  vs.  X  (short-circuit),  (c)  x  vs.  E  (stress  free)  and  (d)  D  vs.  X  (open-circuit)  curves 
for  a  soft  PZT  based  multilayer  actuator  used  for  Figs.  2,  3  and  4.  The  details  on  the 
experiments  will  be  reported  in  the  successive  papers.  Figure  7  shows  the  result  for  the 
piezoelectric  loss  tan0\  We  used  the  correlation  factor  between  electric  field  and 

•  ,  ,,  .  X ,  1/2  _ 

compressive  stress  given  averagely  by  X  =  (£Q£  /s  )  E. 

From  Figs.  6  and  7,  we  can  calculate  the  “intrinsic”  losses  as  shown  in  Fig.  8. 
Note  that  the  piezoelectric  losses  tanO’  and  tan0  are  not  so  small  as  previously  believed,  but 
comparable  to  the  dielectric  and  elastic  losses,  and  increase  gradually  with  the  field  or  stress. 
Also  it  is  noteworthy  that  the  intrinsic  dielectric  loss  tan8  increases  significantly  with  an 
increase  of  the  intensive  parameter,  i.e.,  the  applied  electric  field,  while  the  intrinsic  elastic 
loss  tan<|>  is  rather  insensitive  to  the  intensive  parameter,  i.e.,  the  applied  compressive  stress. 

When  similar  measurements  to  Figs.  1(a)  and  1(b),  but  under  constrained 
conditions;  that  is,  D  vs.  E  under  a  completely  clamped  state,  and  x  vs.  X  under  an  open- 
circuit  state,  respectively,  we  can  expect  smaller  hystereses;  that  is,  intrinsic  losses,  tanS  and 
tan<{>.  These  measurements  seem  to  be  alternative  methods  to  determine  the  three  losses 
separately,  however,  they  are  rather  difficult  in  practice. 
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LOSS  AND  HEAT  GENERATION 


Heat  generation  in  various  types  of  PZT-based  actuators  has  been  studied  under  a 
relatively  large  electric  field  applied  (1  kV/mm  or  more)  at  an  off-resonance  frequency,  and 
a  simple  analytical  method  was  established  to  evaluate  the  temperature  rise,  which  is  very 
useful  for  the  design  of  piezoelectric  high-power  actuators.  Heat  generation  in  a  resonating 
piezoelectric  sample  is  discussed  in  the  following  Section. 

Zheng  et  al.  reported  the  heat  generation  from  various  sizes  of  multilayer  type 
piezoelectric  ceramic  actuators.3*  Figure  9  shows  the  temperature  change  with  time  in  the 
actuators  when  driven  at  3  kV/mm  and  300  Hz,  and  Fig.  10  plots  the  saturated  temperature 
as  a  function  of  Ve/A,  where  Ve  is  the  effective  volume  (electrode  overlapped  part)  and  A  is 
the  surface  area.  This  linear  relation  is  reasonable  because  the  volume  Ve  generates  the  heat 
and  this  heat  is  dissipated  through  the  area  A.  Thus,  if  we  need  to  suppress  the  temperature 
rise,  a  small  Ve/A  design  is  preferred. 

According  to  the  law  of  energy  conservation,  the  rate  of  heat  storage  in  the 
piezoelectric  resulting  from  heat  generation  and  dissipation  effects  can  be  expressed  as 

qg-qout  =  Vpc(d77dt),  (47) 

assuming  uniform  temperature  distribution  in  the  sample.  V,  p,  c  are  total  volume,  density 
and  specific  heat,  respectively.  The  heat  generation  is  considered  to  be  caused  by  losses. 
Thus,  the  rate  of  heat  generation  (qg)  in  the  piezoelectric  can  be  expressed  as 

qg  =  ufVe,  (48) 

where  u  is  the  loss  of  the  sample  per  driving  cycle  per  unit  volume,  f,  the  driving  frequency, 
and  Ve  is  the  effective  volume  where  the  ceramic  is  activated.  According  to  the  measuring 
condition,  this  u  corresponds  to  the  extrinsic  dielectric  loss  we  of  Eq.  (17),  which  consists 
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of  the  intrinsic  dielectric  loss  tan  8  and  the  electromechanical  and  piezoelectric  combined 
loss  (tan<))  -  2  tan0)  in  the  previous  section: 

X  2 

u  =  we  =  n  e  e0  Eq  tan  8’ 

=  [1/(1  -  k2)][tan  8  +  k2(tan  <)>  -2  tan  0)]  %  e  £„  2?0  .  (49) 

Note  that  we  do  not  need  to  add  wem  explicitly,  because  the  corresponding 
electromechanical  loss  is  already  included  implicitly  in  we- 

When  we  neglect  the  conduction  heat  transfer,  the  rate  of  heat  dissipation  (qout) 
from  the  sample  is  the  sum  of  the  rates  of  heat  flow  by  radiation  (qr)  and  convection  (qc): 

qout  =  qr  +  qc 

=  aeA(r4  -  T*)  +  he  A(r  -  T0),  (50) 

where  c  is  the  Stehan-Boltzmann  constant,  e  is  the  emissivity  of  the  sample,  he  is  the 
average  convective  heat  transfer  coefficient,  and  A  is  the  sample  surface  area. 

Thus,  Eq.(47)  can  be  written  in  the  form: 

ufV-  Ak(7)(r-ro)  =  Vpc(dT/dt),  (51) 

where 

k(T)  =  at(f  +  Tq2)(T  +  TJ  +  he  (52) 

is  defined  as  the  overall  heat  transfer  coefficient.  If  we  assume  that  k (T)  is  relatively 
insensitive  to  temperature  change,  the  solution  to  Eq.(51)  for  the  piezoelectric  sample 
temperature  is  given  as  a  function  of  time  (t): 

7  -  r0  =  [u  f  Ve/  k(7)  A]  [1  -  e-tft],  (53) 
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where  the  time  constant  x  is  expressed  as 


x  =  pcV/k(7)A.  (54) 

As  t  — >  oo ,  the  maximum  temperature  rise  in  the  sample  becomes 

AT = u  f  Ve/k(7)  A.  (55) 

As  t  — >  0 ,  the  initial  rate  of  temperature  rise  is 

(dF/dt)  =  u  f  Ve/p  c  V  =  A  77  x  .  (56) 

Figures  11  and  12  show  the  dependence  of  k (7)  on  applied  electric  field  and 
frequency.  Since  k(T)  is  not  really  constant,  we  can  calculate  the  total  loss  u  of  the 
piezoelectric  more  precisely  through  Eq.(56).  The  calculated  results  are  shown  in  Table  I. 
The  experimental  data  of  P-E  hysteresis  losses  under  a  stress-free  condition  are  also  listed 
for  comparison.  It  is  seen  that  the  P-E  hysteresis  extrinsic  loss  agrees  well  with  the  total 
loss  contributing  to  the  heat  generation  under  an  off-resonance  drive. 

LOSSES  AT  A  PIEZOELECTRIC  RESONANCE 

So  far,  we  have  considered  the  losses  for  a  quasi-static  or  off-resonance  state. 
Problems  in  ultrasonic  motors  which  are  driven  at  the  resonance  frequency  include 
significant  distortion  of  the  admittance  frequency  spectrum  due  to  nonlinear  behavior  of 
elastic  compliance  at  a  high  vibration  amplitude,  and  heat  generation  which  causes  a  serious 
degradation  of  the  motor  characteristics  through  depoling  of  the  piezoceramic.  Therefore, 
the  ultrasonic  motor  requires  a  very  hard  type  piezoelectric  with  a  high  mechanical  quality 
factor  Qn,,  leading  to  the  suppression  of  heat  generation.  It  is  also  notable  that  the  actual 
mechanical  vibration  amplitude  at  the  resonance  frequency  is  directly  proportional  to  this  Qm 
value. 


14 


Losses  at  a  Piezoelectric  Resonance 


Piezoelectric  resonance  without  loss 

Let  us  review  first  the  longitudinal  mechanical  vibration  of  a  piezo-ceramic  plate 
without  loss  through  the  transverse  piezoelectric  effect  (d31)  as  shown  in  Fig.  13.4) 
Assuming  that  the  polarization  is  in  the  z-direction  and  the  x-y  planes  are  the  planes  of  the 
electrodes,  the  extentional  vibration  in  the  x  direction  is  represented  by  the  following 
dynamic  equation: 

(d2u/d  t2)  =  F  =  (3X  Jdx)  +  (dXJdy)  +  (dXJdz),  (57) 

where  u  is  the  displacement  of  the  small  volume  element  in  the  ceramic  plate  in  the  x- 
direction.  When  the  plate  is  very  long  and  thin,  X^  and  X^  may  be  set  equal  to  zero  through 
the  plate,  and  the  relations  between  stress,  electric  field  (only  Ez  exists)  and  the  induced 
strain  are  given  by  : 

X  =x/s  E-(d  /s  E)£z-  (58) 

1  1  11  v  31  11  ^ 

Introducing  Eq.(58)  into  Eq.(57),  and  allowing  for  x=du/dx  and  dEz/dx=0  (due  to 
the  equal  potential  on  each  electrode),  leads  to  a  harmonic  vibration  equation : 

-  ©2 p  s  E  u  =  cfuldx  .  (59) 

r  u 

Here,  ©  is  the  angular  frequency  of  the  drive  field,  and  p  is  the  density.  Substituting  a 
general  solution  u=M)(x)ej<B4«2(x)e-ja)t  into  Eq.(58),  and  with  the  boundary  condition  X  = 
0  at  x  =  0  and  L  (sample  length),  the  following  solutions  can  be  obtained: 

(strain)  du/dx  =  x  =  d ^Ez  [sin©(L-x)/v  +  sin(©x/v)]  /sin(©L/v),  (60) 

L 

(total  displacement)  AL  =  J  x  dx  =  d3  £ZL  (2v/©L)  tan(©L/2v).  (61) 

o 
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Here,  v  is  the  sound  velocity  in  the  piezoceramic  which  is  given  by 


v  =  1  /  V  p  s  E . 


Since  the  total  current  is  given  by : 


i  =  jcow  lDdx=jcow  j  [(eQe  x  -  d312/s  *)EZ  +  (d3/suV]  dx. 


and  using  Eq.  (60),  the  admittance  for  the  mechanically  free  sample  is  calculated  to  be: 

Y  =  (1/Z)  =  (i/V)  =  (i/Ezt) 

=  (jcowL/t)  e„e3LC[l  +  (d^/eoe^Cst  iE)(tan(©L72v)  /  (coL/2v)],  (64) 

where  w  is  the  width,  L  the  length,  t  the  thickness  of  the  sample,  and  V  the  applied  voltage. 

LC 

£3  is  the  permittivity  in  a  longitudinally  clamped  sample,  which  is  given  by 

LC  X  ,j  2,  Es  \ 

e  e  =  e  e  -  (d  /s  )  =  en  e  (1  -  k  ) .  (65) 

03  0  3  31  11  0  3  31  ^  V/ 

The  final  transformation  is  provided  by  the  definition, 


k  =  d  /  V  s  e  e  .  (66) 

31  31  II  o  3  v  ’ 

When  the  drive  frequency  is  much  lower  than  the  resonance,  taking  co  ->  0  in  Eq. 

x 

(64)  leads  to  Y  =  (jcowL/t)  e3  (corresponding  to  the  static  capacitance).  The  piezoelectric 
resonance  is  achieved  where  the  admittance  becomes  infinite  or  the  impedance  is  zero.  The 
resonance  frequency  f  is  calculated  from  Eq.  (64),  and  the  fundamental  fR  is  given  by 


f  =  v/2L  =  1/(2L  Vps  ). 
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On  the  other  hand,  the  antiresonance  state  is  generated  for  zero  admittance  or  infinite 
impedance: 

(coL/2v)  cot  (coL/2v)  =  -  d^I  s^  =  -  k3)2/  ( 1  -  k  j2).  (68) 

The  resonance  and  antiresonance  states  are  described  by  the  following  intuitive 
model.4)  In  a  high  electromechanical  coupling  material  with  k  almost  equal  to  1,  the 
resonance  or  antiresonance  states  appear  for  tan(coL/2v)  =  °°  or  0  [i.  e.,  coL/2v  =  (m- 1/2)71  or 
nut  (m:  integer)],  respectively.  The  strain  amplitude  jc  distribution  for  each  state  [calculated 
using  Eq.  (60)]  is  illustrated  in  Fig.  14.  In  the  resonance  state,  large  strain  amplitudes  and 
large  capacitance  changes  (called  motional  capacitance )  are  induced,  and  the  current  can 
easily  flow  into  the  device.  Note  that  for  a  loss-free  piezoelectric  the  strain  is  calculated  to 
be  infinite  in  Eq.  (60).  On  the  other  hand,  at  the  antiresonance,  the  strain  induced  in  the 
device  compensates  completely,  resulting  in  no  capacitance  change,  and  the  current  cannot 
flow  easily  into  the  sample.  Both  ends  of  the  plate  correspond  to  the  nodal  points,  which  do 
not  generate  any  motion  to  be  used  for  actuators.  Thus,  for  a  high  k  material  the  first 
antiresonance  frequency  fA  should  be  twice  as  large  as  the  first  resonance  frequency  fR. 

In  a  typical  case,  where  k3i  =  0.3,  the  antiresonance  state  varies  from  the  above- 

mentioned  mode  and  becomes  closer  to  the  resonance  mode.  The  low-coupling  material 
exhibits  an  antiresonance  mode  where  capacitance  change  due  to  the  size  change  is 
compensated  completely  by  the  current  required  to  charge  up  the  static  capacitance  (called 
damped  capacitance).  Thus,  the  antiresonance  frequency  fA  will  approach  the  resonance  fR. 

When  (f  -  f  )  is  not  very  large  due  to  a  small  electromechanical  coupling,  we  can 

A  R 

derive  the  following  approximate  expression  for  f  .  Assuming  that  ©A-  cdr  is  much  smaller 
than©  (=  7tv /L),  and 

R 
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(70) 


<Oa  =  (tcv/L)[1  +  {Ain)  K3I2], 
where  we  introduced  a  new  parameter  K31  as 

K3,2  =  k3,2/<1-1S,2>-  <71> 

It  is  notable  that  for  a  piezoelectric  sample  with  a  typical  k3|  value,  the  two  ends  of  the  plate 

are  not  the  nodal  points;  that  is,  we  can  expect  rather  large  displacements,  which  can  be 
applied  for  ultrasonic  motors. 

Piezoelectric  resonance  with  losses 

Now,  we  will  introduce  the  complex  parameters  into  the  admittance  curve  around  the 

X*  X  E* 

resonance  frequency,  in  a  similar  way  to  the  previous  section:  e3  =  e3  (1  -  j  tan  8’),  s 
=  Sj  jE(  1  -  j  tan  <t>’),  and  d3i*  =  d(l  -  j  tan  0’)  into  Eq.  (64): 

Y  =  Y  +  Y 

d  m 

=  0'c»wL/t)eoe3X(l  -k  [2)[1  -  j  (l/(l-k3i2))((tan  S’+  k_2(tan  <(>’  -  2  tan  0’))] 

+  (jQ)wL/t)e0e3Xk3i2  [(1  -  j(2  tan  0’-  tan  <f>’)][(tan(©L/2  v*)  /  (coL/2v*)] 

=  jc»Co  (1  -k3i2)  [1  -  j  (l/(l-k3i2))((tan  8’+  k3]2(tan  (j)’  -  2  tan  0’))] 

+ jo yCQ  k3[2[(l-j(2tan0’-tan  <f>’)][(tan(©L/2v*)/(©L/2v*)] 

=  jcoC  (1  -  j  tan  8)  +  j©C  K  2[(l-j(2tan0’-tan  <t>’)][(tan(©L/2v*)/(©L/2v*)],  (72) 

d  d  31 


where 


(73) 


Co=  (wL/t)  eoe3X, 

C  =  (1  -  k  2)  C  .  (74) 

Note  that  the  loss  for  the  first  term  (damped  conductance)  is  represented  by  the  “intrinsic” 
dielectric  loss  tan  8,  not  by  the  extrinsic  loss  tan  8’.  Taking  into  account 

v*  =  1/Vp  S(iE(l  -  j  tan  <t>’)  =  v  (1  +  (1/2)  j  tan  <|)’),  (75) 

we  further  calculate  l/(tan(coL/2v*)  with  an  expansion-series  approximation  around  (coL/2v) 
=  7C/2.  The  resonance  state  is  defined  in  this  case  for  the  maximum  admittance  point,  rather 
than  the  infinite  Y. 

We  will  use  new  frequency  parematers, 

Q.  =  toL/2  v,  AD  =  £2  -  rc/2  («1) .  (76) 

Since  coL/2  v*=  (n/2+  AD)  [1  -  (1/2)  j  tan  <t>’], 

1  /  (tan(a>L/2v*)  =  -  AD  +  j  (7t/4)tan  <)>’  .  (77) 

Thus,  using  K  2  =  k  2/  (1  -  k  2),  the  motional  admittance  Y  is  approximated  aound  the 

°  31  31  31  m 

first  resonance  frequency  by 

Y  =jcoC  K,,2  [(l-j(2tan0’-tan  (|)’)][(tan((oL/2v*)/(o)L/2v:*:)] 

m  d  31 

= jco0Cd  Km2  [(l-j(2tan9’-tan  <t)’)]/[(-A^+j(7c/4)tan<}),)(7c/2)(  1  -( l/2)j tan<t>’ )] 

=  j(7t2/8)  cooCd  K3|2  [(l+j((3/2)tan<j>’-2tane’)]/[(-  (4/7t)An+jtan(()’).  (78) 

The  maximum  Y  is  obtained  at  AD  =  0: 

m 

Y  m“  =  (k2/8)  co  C  K  2  (tan  <J>’) 

m  u  d  31 


(79) 


In  order  to  obtain  the  mechanical  quality  factor,  let  us  obtain  AH  which  provides  Y  m“/V2. 

m 

Since  AQ  =  (n/4)tan  <t>’  is  obtained, 

Q  =  Q  /2AQ  =  (k/2)/  2  (7t/4)tan  <J>’=  (tan  (80) 

m  0 

This  verifies  the  ready-used  relation,  Q  =  (tan  <J>’) 

m 

Here,  the  displacement  amplification  is  also  considered.  From  Eq.  (61),  also  by 
using  the  complex  parameters: 

w(L)  =  d31*£zL  (2v*/coL)  tan(o)L/  2v*) 

=  2d3](l-j  tanG’)  2szL  [v(l+(l/2)  j  tan<J>’)  /coL]  tan(coL/2  v*) 

=  2d31(l-j  tan6’)  EZL  [v(l+(l/2)  j  tan<}>’)/cDoL]  /  (-  AQ  +  j(rc/4)tan  <)>’).  (81) 

The  maximum  displacement  u  is  obtained  at  AQ  =  0: 

max 

“max =  ^31  (tan  (()’)  (82) 

2 

The  maximum  displacement  at  the  resonance  frequency  is  (k  /8)  Q  times  larger  than  that  at 

m 

a  non-resonance  frequency  (d3]£zL). 

In  a  brief  summary,  we  obtained  a  ready-used  knowledge:  when  we  observe  the 
admittance  or  dislpacement  spectrum  as  a  function  of  drive  frequency,  and  obtain  the 
mechanical  quality  factor  Qm  estimated  from  Qm  =  co/2Aco,  where  2A(0  is  a  full  width  of  the 

3  dB  down  (i.e.,  1/V2)  of  the  maximum  value  at  ©  =  ©q,  we  can  obtain  the  extrinsic 
mechanical  loss  tan  <{>’. 

Equivalent  Circuit 

The  equivalent  circuit  for  the  piezoelectric  acuator  is  represented  by  a  combination 
of  L,  C  and  R.  Figure  15  (a)  shows  an  equivalent  circuit  for  the  resonance  state,  which  has 
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very  low  impedance.  Taking  into  account  Eq.  (72),  we  can  understand  that  Cd  and  Rd 
correspond  to  the  electrostatic  capacitance  (for  a  longitudinally  clamped  sample  in  the 
previous  case,  not  a  free  sample)  and  the  clamped  (or  “intrinsic”)  dielectric  loss  tan  5, 
respectively,  and  the  components  La  and  Ca  in  a  series  resonance  circuit  are  related  to  the 
piezoelectric  motion.  For  example,  in  the  case  of  the  longitudinal  vibration  of  the  above 
rectangular  plate  through  d3[,  these  components  are  represented  approximately  by 


LA  =  (p/8)(Lb/w)(siiE2/d3i2), 

(83) 

CA  =  (8  /  7t2)(Lw  /  b)(d3J2  /  snE). 

(84) 

The  total  resistance  R  (=  R  +  R  )  should  corresponds  to  the  loss  tan  <{>\  which  is 
composed  of  the  intrinsic  mechanical  loss  tan  <J>  and  dielectric/piezoelectric  coupled  loss  (tan 
8  -  2tan  0)  (see  Eq.  (45)).  Thus,  roughly  speaking,  Rd  and  Rm  correspond  to  the  intrinsic 
dielectric  and  mechanical  losses,  respectively.  Note  that  we  introduced  an  additional 
resistance  Rd  to  explain  a  large  contribution  of  the  dielectric  loss  when  a  vibration  velocity 
is  relatively  large.  In  contrast,  the  equivalent  circuit  for  the  antiresonance  state  of  the  same 
actuator  is  shown  in  Fig.  15  (b),  which  has  high  impedance. 

Losses  as  a  Function  of  Vibration  Velocity 

Figure  16  shows  the  mechanical  versus  basic  composition  x  at  two  effective 
vibration  velocities  vq=0.05  m/s  and  0.5  m/s  for  Pb(Zr^Tij  ^)03  doped  with  2.1  at.%  of  Fe.5) 
The  decrease  in  mechanical  with  an  increase  of  vibration  level  is  minimum  around  the 
rhombohedral-tetragonal  morphotropic  phase  boundary  (52/48).  In  other  words,  the 
smallest  Qm  material  at  a  small  vibration  level  becomes  the  best  at  a  large  vibration  level,  and 
the  data  obtained  by  a  conventional  impedance  analyzer  with  a  small  voltage/power  does  not 
provide  data  relevant  to  high  power  materials. 
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Let  us  consider  here  the  degradation  mechanism  of  the  mechanical  quality  factor 
with  increasing  vibration  velocity.  Figure  17  shows  an  important  notion  on  heat  generation 
from  the  piezoelectric  material.®  The  damped  and  motional  resistances,  R  and  R  ,  in  the 
equivalent  electrical  circuit  of  a  PZT  sample  are  separately  plotted  as  a  function  of  vibration 
velocity.  Note  that  Rm,  mainly  related  to  the  intrinsic  mechanical  loss,  is  insensitive  to  the 
vibration  velocity,  while  Rd,  related  to  the  intrinsic  dielectric  loss,  increases  significantly 
around  a  certain  critical  vibration  velocity.  Thus,  the  resonance  loss  at  a  small  vibration 
velocity  is  mainly  determined  by  the  intrinsic  mechanical  loss  which  provides  a  high 
mechanical  quality  factor  and  with  increasing  vibration  velocity,  the  intrinsic  dielectric 
loss  contribution  significantly  increases.  After  RJ  exceeds  R  ,  we  started  to  observe  heat 

d  m 

generation. 

S.  Tashiro  et.  al  observed  the  heat  generation  in  a  rectangular  piezoelectric  plate 
during  a  resonating  drive.7)  Eventhough  the  maximum  electric  filed  is  not  very  large,  heat  is 
generated  due  to  the  large  induced  strain  and  stress  at  the  resonance.  The  maximum  heat 
generation  was  observed  at  the  nodal  point  of  the  resonance  vibration,  where  the  maximum 
strain  and  stress  are  generated.  This  observation  supports  that  the  heat  generation  in  a 
resonating  sample  is  attributed  to  the  “extrinsic”  mechanical  loss  tan  <j>\  This  is  not 
contradictory  to  the  result  in  Section  4,  where  a  high-voltage  drive  was  conducted  at  an  off- 
resonance  frequency.  We  concluded  there  that  the  heat  is  originated  from  the  “extrinsic” 
dielectric  loss  tan  5’.  In  consideration  that  both  the  “extrinsic”  dielectric  and  mechanical 
losses  are  composed  of  the  “intrinsic”  dielectric  and  mechanical  losses,  and  that  the 
intrinsic  dielectric  loss  tan  5  changes  significantly  with  the  external  electric  field  and  stress, 
the  major  contribution  to  the  heat  generation  seems  to  come  from  the  “intrinsic”  dielectric 
loss.  Further  investigations  are  waited  for  the  microscopic  explanations  of  this  phenomon. 
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LOSSES  AT  RESONANCE  AND  ANTIRESONANCE  MODES 


Losses  at  a  Piezoelectric  Antiresonance  State 

We  consider  here  the  losses  at  the  antiresonance  frequency  in  comparison  with  the 
resonance  mode.  The  antiresonance  mode  is  obtained  at  a  frequency  which  provides  the 
minimum  value  of  admittance  Y,  instead  of  zero  of  Y  for  the  loss-free  case.  Taking  an 
approximation  technique  on  Eq.  (72)  around  the  antiresonance  frequency  coa,  similar  to  the 
previous  section,  we  obtain: 

£2a  =  coaL/2  v,AQ  =  Q-Qa  («1) .  (85) 

If  k31  is  not  very  large,  the  following  relationship  is  obtained: 

ft  =  (0aL/  2v  =  (jc/2)  (1  +  (4/tc2)  k  2)  .  (86) 

A  A  31 

In  the  following  approximation,  however,  this  relation  is  not  used,  but  we  will  neglect  the 
higher  order  of  Ml  and  tan  <t>’  in  Eq.  (72): 

Y  =  jcoCJ  1  +  (k  2/(l-k  2))  tan[(Q  +A£2)  (1- j  (1/2)  tan<t>’)]/ 

0  31  31  ^ 

[(Qa+M2)  (1- j  (1/2)  tan<t>’)].  (87) 

Taking  into  account 

tan[(£2A+A£2)  (1- j  (1/2)  tan<t>’)] 

=  [(nA2-K3i2AQ+  j(l/2)i2Atan<|)’K3i2]/[(-K312-  Cl^Ml)  -jQA(l/2)flAtan<i>’],  (88) 

Y  =  j©Cd  [-  (Q\k\k*)AQ.  -  j(l/2)QAtan(|>’  (Q^-K^-K  /)/ 

[-  K3]2Qa  -  (nA2+K3i2)AD  -  j(1/2)Qa  tan<|)’  (Q*-  K^)],  (89) 

where 
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Then,  Y™"  can  be  obtained  at  A £2  =  0: 

Y"”"  =  toCd  (l/2)tan4>’  (Q*-  K*  -  K  /  )/K  (2.  (91) 

V2  Y™"  can  be  obtained  at 

A£2  =  (1/2)  Qa  tan*’  (QA2-  K  ^  K  *)  /(QA  +  K /).  (92) 

Thus,  mechanical  quality  factor  at  the  antiresonance  can  be  obtained  as: 

Qm  =  Oa  /2A£2  =  [(Oa  +  K  (2+  K  /)  /(^A2-  K,2-  K /)]  (tantj)’)''.  (93) 

Since  (Q  +  K  +  K  )/(C2  -  K  -  K  )  is  larger  than  1,  Q  can  be  verified  to  be 

A  31  31  a  31  31  m 

larger  than  Qm  at  the  resonance  (=  (tan<|>’) ').  When  k3|  is  small,  using  Eq.  (86)  and 

4 

neglecting  k3j  or  higer  orders,  Qm  is  approximated  as 

Qm  =  (1  +  (8/tt2)  k3j2)  (tantj)’)*'.  (94) 

Note  that  this  Qm  is  equal  to  the  inverse  value  of  the  intrinsic  mechanical  loss,  (tantf))  '. 

Experimental  Results 

Figure  18  illustrates  mechanical  quality  factors,  Qa,  Qg  and  the  temperature  rise  for 
the  resonance  (A-type)  and  the  antiresonance  (B-type)  modes  for  a  rectangular-shape  har 
PZT  resonator  plotted  as  a  function  of  vibration  velocity .6)  The  sample  size  is  indicated  in 
the  figure  (43  mm  x  7  mm  x  2  mm).  Note  that  an  “effective”  vibration  velocity  v0  is  a 
material’s  constant  independent  of  the  sample  size,  and  is  defined  as  V2  7t  f  u  where  f  is 

max 

the  resonance  or  antiresonance  frequency  and  u  is  the  maximum  vibration  amplitude  of 
the  piezoelectric  device.8)  It  is  again  noteworthy  that  the  mechanical  quality  factor  decreases 


significantly  above  a  certain  critical  vibration  velocity  (0.1  m/s),  where  a  steep  tempreature 
rise  starts.  We  have  suggested  that  the  heat  generation  is  mainly  attributed  to  the  intrinsic 
dielectric  loss  rather  than  the  mechanical  loss.  Note  also  that  QB  is  higher  than  QA  over  the 
entire  vibration  velocity  range,  and  that  the  temperature  rise  of  the  sample  is  less  for  the  B- 
type  resonance  (antiresonance)  than  for  the  A-type  resonance  for  the  same  vibration  level. 
This  indicates  an  intriguing  idea  that  the  antiresonance  mode  should  be  superior  to  the 
conventional  resonance  mode,  particularly  for  high  power  applications  such  as  ultrasonic 
motors.  In  a  typical  piezoelectric  material  with  k3]  around  30  %,  the  plate  edge  is  not  a 

vibration  nodal  point  and  can  generate  a  large  vibration  velocity. 

CONCLUSIONS 

(1)  Various  techniques  for  measuring  the  electric,  mechanical  and  piezoelectric  coupling 
losses  separately  have  been  discussed: 

(a)  D  vs.  E,  x  vs.  X,  x  vs.  E  and  D  vs.  X  curves  for  dielectric,  mechanical  and 
piezoelectric  losses 

(b)  heat  generation  at  a  resonance  or  an  off-resonance  frequency  for  an  extrinsic 
mechanical  or  dielectric  loss 

(c)  resonance/antiresonance  technique  for  extrinsic  and  intrinsic  mechanical  losses, 
respectively 

By  combining  the  above  methods,  we  can  investigate  the  loss  mechanisms  in  practical 
piezoelectric  materials. 

(2)  The  piezoelectric  losses  tan0’  and  tan0  are  not  so  small  as  previously  believed,  but 
comparable  to  the  dielectric  and  elastic  losses  in  soft  PZTs.  Also  it  is  noteworthy  that 
the  intrinsic  dielectric  loss  tan8  increases  significantly  with  an  increase  of  the  intensive 
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parameter,  i.e.,  the  applied  electric  field,  while  the  intrinsic  elastic  loss  tan<J>  is  rather 
insensitive  to  the  intensive  parameter,  i.e.,  the  applied  compressive  stress. 

(3)  Heat  generation  is  caused  mainly  by  the  extrinsic  dielectric  loss  tan  8’  (i.e.,  P-E 
hysteresis  loss)  for  an  off-resonance  state  under  a  high  drive  electric  field,  and  by  the 
extrinsic  mechanical  loss  tan  4>’  for  a  resonance  state.  Both  situations  are  attributed  to 
the  large  “intrinsic  dielectric  loss”  enhanced  by  a  large  external  electric  field  or  stress. 
In  order  to  suppress  the  temperature  rise  practically,  a  transducer  design  with  larger 
surface  area  is  recommended  (for  example,  a  tube  rather  than  a  rod). 

(4)  A  significant  decrease  in  mechanical  with  an  increase  of  vibration  level  was  observed 

in  resonant  piezoelectric  ceramic  devices,  and  the  data  obtained  by  a  conventional 
impedance  analyzer  with  a  small  voltage/power  do  not  provide  data  relevant  to  high 
power  materials. 

(5)  Since  the  mechanical  quality  factor  QB  at  an  antiresonance  frequency  is  larger  than  Qa  at 
a  resonance  frequency,  the  antiresonance  mode  seems  to  be  superior  to  the  conventional 
usage  of  the  resonance  mode,  particularly  for  high  power  applications  such  as  ultrasonic 
motors. 
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Table  and  Figure  Captions 

Table  I  Loss  and  overall  heat  transfer  coefficient  for  PZT  multilayer  samples  (E  =  3  kV/mm, 
f  =  300  Hz). 

Fig.  1  (a)  D  vs.  E  (stress  free),  (b)  *  vs.  X  (short-circuit),  (c)  x  vs.  E  (stress  free)  and  (d)  D 

vs.  X  (open-circuit)  curves  with  a  slight  hysteresis  in  each  relation. 

Fig.2  Loss  tan  5’  as  a  function  of  sample  temperature  (3  kV/mm,  300  Hz). 

Fig.3  Loss  tan  8’  as  a  function  of  electric  field  (T  =  25°C,  f  =  300  Hz). 

Fig.4  Loss  tan  8’  as  a  funtion  of  frequency  ( T  =  25°C,  E  =  2  kV/mm). 

Fig.5  Conceptual  figure  for  explaining  the  relation  between  £X  and  e\  sE  and  s°. 

Fig.6  Dissipation  factors  determined  from  (a)  D  vs.  E  (stress  free),  (b)  x  vs.  X  (short- 

circuit),  (c)  x  vs.  E  (stress  free)  and  (d)  D  vs.  X  (open-circuit)  curves  for  a  PZT 

based  actuator. 

Fig.7  Extrinsic  piezoelectric  dissipation  factor  tan0’  as  a  function  of  electric  field  or 
compressive  stress,  measured  for  a  PZT  based  actuator. 

Fig.  8  Intrinsic  loss  factors,  tanS,  tan<|)  and  tan0  as  a  function  of  electric  field  or 
compressive  stress,  measured  for  a  PZT  based  actuator. 

Fig.9  Temperature  rise  for  various  actuators  while  driven  at  300  Hz  and  3  kV/mm. 

Fig.  10  Temperature  rise  versus  Ve/A  (3  kV/mm,  300  Hz),  where  Ve  is  the  effective  volume 
generating  the  heat  and  A  is  the  surface  area  dissipating  the  heat. 

Fig.  11  k(T)  as  a  function  of  applied  electric  field  (400  Hz,  data  from  the  actuator  with 
dimensions  of  7  mm  x  7  mm  x  2  mm). 
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Fig.  12  Overall  heat  transfer  coefficient  k(7)  as  a  function  of  frequency. 

Fig.  13  Longitudinal  vibration  through  the  transverse  piezoelectric  effect  (c^,)  in  a 
rectangular  plate. 

Fig.  14  Strain  generation  in  the  resonance  or  antiresonance  state.  The  strain  magnitude  is 
plotted  in  the  vertical  axis  as  a  function  of  the  x  coordinate. 

Fig.  15  Equivalent  circuit  of  a  piezoelectric  device  for  the  resonance  (a)  and  the 
antiresonance  (b). 

Fig.  16  Mechanical  Qm  versus  basic  composition  x  at  two  effective  vibration  velocities 
vo=0.05  m/s  and  0.5  m/s  for  PtyZrT^  x)03  doped  with  2.1  at.%  of  Fe. 

Fig.  17  Vibration  velocity  dependence  of  the  resistances  Rd  and  Rm  in  the  equivalent  electric 
circuit  for  a  PZT  sample. 

Fig.  18  Vibration  velocity  dependence  of  the  quality  factor  Q  (QA,QB)  and  temperature  rise 
for  both  A  (resonance)  and  B  (antiresonance)  type  resonances  of  a  longitudinally 
vibrating  PZT  ceramic  transducer  through  the  transverse  piezoelectric  effect  d3]. 
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Table  I  Loss  and  overall  heat  transfer  coefficient  for  PZT  multilayer  samples 
(E=3  kV/mm,  f  =  300  Hz). 


Actuator 

4.5x3 .5x2  mm 

7x7x2  mm 

17x3.5x1  mm 

Total  loss  (xl03J/m3) 

u-SVT) 

Qt  t->0 

19.2 

19.9 

19.7 

P-E  hysteresis  loss 

(xl03J/m3) 

18.5 

17.8 

17.4 

k(7)  (W/m2K) 

38.4 

39.2 

34.1 
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(d)  D  vs.  X  (open-circuit)  curves  with 
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Q  >"«««»  i  »  t  j  i  l  i  »  i  t  1  i  i  i  i  I  t  i  i  i  1  i  i  i.-i-l 

0  0.5  1  1.5  2  2.5  3 

Electric  Field  (kV/mm) 


Fig.3  Loss  tan  8’  as  a  function  of  electric  field  (T  =  25°C,  f  =  300  Hz). 
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Frequency  (kHz) 


Fig.4  Loss  tan  5’  as  a  funtion  of  frequency  ( T  =  25°C,  E  =  2  kV/mm). 
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clamped  Electrically 
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Itmff  -  tan'f' 


Electric  Field  (kV/mm)  Compressive  Stress  (MPa) 


Fig.  6  Dissipation  factors  determined  from  (a)  D  vs.  E  (stress  free),  (b)  x  vs.  X  (short- 
circuit),  (c)  x  vs.  E  (stress  free)  and  (d)  D  vs.  X  (open-circuit)  curves  for  a  PZT 
based  actuator. 
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Electric  Field  (kV/mm) 


Fig.7  Extrinsic  piezoelectric  dissipation  factor  tan0’  as  a  function  of  electric  field  or 
compressive  stress,  measured  for  a  PZT  based  actuator. 
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Electric  Field  (kV/mm) 


Compressive  Stress  (MPa) 


Fig.8  Intrinsic  loss  factors,  tan8,  tan<]>  and  tanS  as  a  function  of  electric  field  or 
compressive  stress,  measured  for  a  PZT  based  actuator. 
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Actuator  Temperature  (°C) 
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Fig.  10  Temperature 


rise  vi 


f  (kHz) 


Fig.  12  Overall  heat  transfer  coefficient  k(T)  as  a  function  of  frequency. 
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Fig.  13  Longitudinal  vibration  through  the  transverse  piezoelectric  effect  (d3])  in  a 
rectangular  plate. 
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Resonance 


Antiresonance 


m  =  l 


Low  coupling 


High  coupling 


m  =  2  m  =  2 


Fig.  14  Strain  generation  in  the  resonance  or  antiresonance  state.  The  strain  magnitude  is 
plotted  in  the  vertical  axis  as  a  function  of  the  x  coordinate. 


Mechanical  quality  fact  or  Q 


Pb(ZrxTii_x)03  +  2.1  at%  Fe 


0.48  0.50  0.52  0.54  0.56  0.58 


Mole  Fraction  of  Zr  (x) 


Fig.  16  Mechanical  Qm  versus  basic  composition  x  at  two  effective  vibration  velocities 
v0=0.05  m/s  and  0.5  m/s  for  Pb(ZrTi  x)03  doped  with  2.1  at.%  of  Fe. 
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Mechanical  quality  fact  or  Q 
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Fig.  17  Vibration  velocity  dependence  of  the  resistances  Rrf  and  Rm  in  the  equivalent  electric 
circuit  for  a  PZT  sample. 
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Fig.  18  Vibration  velocity  dependence  of  the  quality  factor  Q  (QA,QB)  and  temperature  rise 
for  both  A  (resonance)  and  B  (antiresonance)  type  resonances  of  a  longitudinally 
vibrating  PZT  ceramic  transducer  through  the  transverse  piezoelectric  effect  d3). 
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Temperature  Rise  (OC) 


